Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,906,608 G→C 5.1% intergenic (+25/+35) ygcG → / ← eno TPM domain protein, putative phosphatase/enolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,906,6080GC 5.1% 220.1 / 4.1 78intergenic (+25/+35)ygcG/enoTPM domain protein, putative phosphatase/enolase
Reads supporting (aligned to +/- strand):  ref base G (27/47);  new base C (2/2);  total (29/49)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

GAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAG  >  NC_000913/2906460‑2906755
                                                                                                                                                    |                                                                                                                                                   
gAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCg                                                                                                                                                     >  2:482294/1‑149 (MQ=255)
gAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTAACCCATAAAAAATGCCATCCCg                                                                                                                                                     >  2:499427/1‑149 (MQ=255)
 aaaaGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTATCTCAACCATAAAAAATGCCAGCCCgg                                                                                                                                                    >  2:52341/1‑149 (MQ=255)
 aaaaGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCAGGAGGCTCAGATAGCGGCGTATTTAGCGGCGGGGGCGGATCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCACCCATAAAAAATGCCAGACCgg                                                                                                                                                    >  2:549243/1‑149 (MQ=255)
  aaaGAGCATCTTGGAGAAGTTATCATTCAGACAACAGCGGTTCCGCTTGAGGCTCAGATAGCGGCGGCGTTAGAGGCGGGGGCGGTTCTTCTGGAGGGGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGa                                                                                                                                                   <  1:761605/149‑1 (MQ=255)
           cTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAt                                                                                                                                          <  1:991751/149‑1 (MQ=255)
             ggtactaGTTATCATTCTGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCAttt                                                                                                                                        <  1:233227/143‑1 (MQ=255)
             ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAttt                                                                                                                                        >  2:435231/1‑149 (MQ=255)
             ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCGGGAGGCTCAGATAGCGGCGGATTTCGCGGCGGGGGCGGTTCTTCTGGAGGGGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGGGGCTGGCAttt                                                                                                                                        <  2:874579/149‑1 (MQ=255)
              gTAGAAGTTATCATTCAGACAACAGCGGGTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAtttt                                                                                                                                       <  2:986376/149‑1 (MQ=255)
               tcgAAGGTCTCATTCAGACAACAGCCGTTCCGATGGAGGCGCAGATCGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGGGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCAttttt                                                                                                                                      <  2:872679/147‑1 (MQ=255)
                 gAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGCGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTACCACAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTaa                                                                                                                                    >  1:348333/1‑149 (MQ=255)
                  aaGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTaaa                                                                                                                                   <  2:725915/149‑1 (MQ=255)
                   aggtATCATGCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGCTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAAt                                                                                                                                  <  1:358867/146‑1 (MQ=255)
                       aTCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGTGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCAGGAGGCAGGCATTTTTAAATCAGa                                                                                                                              >  1:1012955/1‑149 (MQ=255)
                       aTCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGa                                                                                                                              >  2:154585/1‑149 (MQ=255)
                          aTTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaa                                                                                                                           >  1:313174/1‑149 (MQ=255)
                          aTTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCCGCGGGGGCGGTTCTTCAGGAGGTGGCGGCGCATCCCGCCGCTGGTAACTCAACCATAAAAAATGCCCGCCCGGAGGCGGGCATTGTTAAATCAGATaa                                                                                                                           >  1:707941/1‑149 (MQ=255)
                           ttCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGTCATTGTTAAATCAGATaaa                                                                                                                          >  1:954353/1‑149 (MQ=255)
                           ttCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaaa                                                                                                                          >  1:530543/1‑149 (MQ=255)
                                  tccaGCGGAACCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTTACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAagtcagtc                                                                                                                   >  1:770283/3‑149 (MQ=255)
                                  aacaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAagtcagtc                                                                                                                   >  1:626114/1‑149 (MQ=255)
                                   acaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCt                                                                                                                  <  1:149767/149‑1 (MQ=255)
                                     aGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCGTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATATAGTCAGTCttt                                                                                                                >  1:760298/1‑148 (MQ=255)
                                        ggTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGc                                                                                                             >  1:38211/1‑149 (MQ=255)
                                         gTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGcc                                                                                                            >  1:99285/1‑149 (MQ=255)
                                         gTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGcc                                                                                                            >  1:703208/1‑149 (MQ=255)
                                               cTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGAGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTAAGCCTGGCCt                                                                                                      <  2:745993/149‑1 (MQ=255)
                                                 ggAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCttt                                                                                                    <  2:718662/149‑1 (MQ=255)
                                                  gAGGCGCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTg                                                                                                   <  2:262056/149‑1 (MQ=255)
                                                    ggCTCAGTTAGCGGCGGATTTAGCGGCGGCGGCGGGTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGTCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTCGCCTTTGAt                                                                                                 <  2:855141/149‑1 (MQ=255)
                                                    ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAt                                                                                                 <  1:579908/149‑1 (MQ=255)
                                                     gCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtc                                                                                                <  1:370642/149‑1 (MQ=255)
                                                     gCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtc                                                                                                <  1:47746/149‑1 (MQ=255)
                                                      gtctgttAGCGGCGTTTTTAGCGGCGGGGGGGGTTCTTCTGGAGGTGGCGGCGCATGCGGGCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTATATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGTtct                                                                                               <  1:554840/143‑1 (MQ=255)
                                                        cAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATctct                                                                                             <  2:516593/149‑1 (MQ=255)
                                                        cAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATctct                                                                                             <  2:598809/149‑1 (MQ=255)
                                                        cAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATctct                                                                                             <  1:127109/149‑1 (MQ=255)
                                                         aGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCtt                                                                                            <  1:376917/149‑1 (MQ=255)
                                                          gATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCttt                                                                                           <  2:703208/149‑1 (MQ=255)
                                                          gATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCttt                                                                                           <  2:99285/149‑1 (MQ=255)
                                                           aTAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTa                                                                                          <  1:999419/149‑1 (MQ=255)
                                                           aTAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCATGGCCTTTGATCTCTTTa                                                                                          >  2:832476/1‑149 (MQ=255)
                                                           aTAGCGGCGGATTTAGAGGCGTGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACGATAAAAAATGCCAGCCCGTAGGCTGGCATTTTTAAATCAGATAAAGTCTGTCTTATGCCTGGCCTTTGATCTCTTTa                                                                                          <  2:141542/149‑1 (MQ=255)
                                                               cggcggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGAc                                                                                      <  2:791779/149‑1 (MQ=255)
                                                                   ggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGtt                                                                                  <  2:776594/149‑1 (MQ=255)
                                                                   ggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGtt                                                                                  <  1:652779/149‑1 (MQ=255)
                                                                         aGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAGTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACgg                                                                            >  1:606522/1‑149 (MQ=255)
                                                                            ggcggGGGCGGATCTTATGAAGGTGGCGGCGCATCCGGCCACTGGTAACTCAACCATAAAAAATGCCAGCCCGTAGTCTGGCATTTTTAAAACAGATAAAGTCATTCTTATGCCTGTCCTTTGATATCTTTACGTCCTTTGTACTGTGc                                                                         >  2:618990/1‑149 (MQ=255)
                                                                             gcggGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGTCCGCTGGTAACTAAACCATAAAAAATTCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTAATGCCTGGCCTTTGATCTCTTTACGACCGTAGTACGGTGCt                                                                        >  2:683991/1‑149 (MQ=255)
                                                                              cggGGGCGGTTCTTCTGGAGGTGGCGGGGCATCCGGCGGCTGGTAAGTCAACGATAAAATATGGCAGCCCGGAGGCTGGCATTTTTAAATCATATAAAGGCAGTGTTATGCCTGTCCTTTGATCTCTTTACGACCGTTGTACGGTGCtt                                                                       <  1:62187/149‑1 (MQ=255)
                                                                                 gggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCttttt                                                                    >  2:289065/1‑149 (MQ=255)
                                                                                 gggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCttttt                                                                    <  1:435231/149‑1 (MQ=255)
                                                                                  ggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTc                                                                   >  2:31809/1‑149 (MQ=255)
                                                                                   gCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCTTTGTACGGTGCTTTTTcc                                                                  >  2:37069/1‑148 (MQ=255)
                                                                                    cGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGc                                                                 <  1:29281/149‑1 (MQ=255)
                                                                                     ggTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGcc                                                                <  2:489698/149‑1 (MQ=255)
                                                                                          ttctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCagag                                                           >  1:104480/1‑149 (MQ=255)
                                                                                           tctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATTAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGc                                                          >  2:585315/1‑149 (MQ=255)
                                                                                               gAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttct                                                      <  1:832355/149‑1 (MQ=255)
                                                                                                aGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttcct                                                     <  2:512862/149‑3 (MQ=255)
                                                                                                 ggTGGCGGCGCATCCGTCCGCTGGTAACTCAACCATTAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATATAGTCAGTCTTATGCCTGGCCTTTTATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCttcttc                                                    <  1:549243/149‑1 (MQ=255)
                                                                                                 ggTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGGCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCttcttc                                                    <  2:370363/149‑1 (MQ=255)
                                                                                                 ggTGGCGGCGCATCCGGACGCTGGTAACTAAACCATAAAAAATGCCAGCCCGGATGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCttcttc                                                    <  2:856102/149‑1 (MQ=255)
                                                                                                    ggcggcGCATCCGGCCGCTGGTAACACAACCATAAAAAATGCCAGCACGGAGGCTGGCATTTTTAAATCAGATAAAGTCACTCTTATGACTGGCCTATGATCTCTTTACGACCGTTGTACGGTGCATTATAGCACATAGCTACTTCGaa                                                 >  2:310290/1‑148 (MQ=255)
                                                                                                     gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATa                                                <  2:305787/149‑1 (MQ=255)
                                                                                                     gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATa                                                <  1:482294/149‑1 (MQ=255)
                                                                                                     gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGACAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATa                                                <  2:114509/149‑1 (MQ=255)
                                                                                                      cggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATAc                                               <  2:626114/149‑1 (MQ=255)
                                                                                                          gcATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAAt                                           <  1:98457/149‑1 (MQ=255)
                                                                                                              ccGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGc                                       <  1:499427/149‑1 (MQ=255)
                                                                                                                   cGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGtt                                  <  1:787629/149‑1 (MQ=255)
                                                                                                                             tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAGTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTGTGATCTCTTTACGACCGTTGGACGGTGCTTTTTCGCCCAGCGCTTCTTCGCTACGAATCAGCGGGATGTATTTAGCa                        >  2:451209/1‑149 (MQ=255)
                                                                                                                                  cATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACg                   <  2:532017/149‑1 (MQ=255)
                                                                                                                                   tttaaaaaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGATTCAGCTGGTTGTATTTAGCAACACgg                  >  2:260547/4‑149 (MQ=255)
                                                                                                                                     aaaaaaTGCCTGGCCGGTGGCGGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTTATCTCTTTACGACCGATGTACGGTGCTTTTTCGCCCAGAGCTTATTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTc                <  2:815370/149‑1 (MQ=255)
                                                                                                                                     aaaaaaTGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGGCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTACGATACGAATCAGCTGGTTGTATTTAGCAACACGGTc                <  1:620225/149‑1 (MQ=255)
                                                                                                                                      aaaaaTGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGGCAGTCTTATGCCTGGCCTTTGAGCTCTTTACGATCGTTGTCCGGGGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACCCGGTCa               >  2:100888/1‑149 (MQ=255)
                                                                                                                                       aaaaTGCCAGCCCGGCGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACAGTTGTACGGTGCTTTTTCGCCCTGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTGAg              >  1:103127/1‑149 (MQ=255)
                                                                                                                                        aaaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGa             >  2:105425/1‑149 (MQ=255)
                                                                                                                                             ccAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGc        <  2:131749/149‑1 (MQ=255)
                                                                                                                                              acgCCCGGAGGCTGGCATTTTTAAATCAGAATCAGTCAGGCTTATGCCTGGCCTTTGATCTCTTTCCGACCGTTGTACGGTGCGTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCt       <  1:504500/147‑1 (MQ=255)
                                                                                                                                                  caggAGGCTGGCATTTTTAATTTAGATAAAGTCAGTCTTATGCCTGGCCTTTGTTCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGGAACACGGTCAGAACGGCTCATa   <  2:707941/147‑1 (MQ=255)
                                                                                                                                                   cggcggCTGGCATTATTAAATAAGTTAAAGTCAGTATTATTCCTGGCCTTTGATATCTTTAAGACCGTTGAACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAg  <  2:954353/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAG  >  NC_000913/2906460‑2906755

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑