Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,108 A→C 12.0% intergenic (+91/‑107) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1080AC12.0% 169.6 / 2.9 92intergenic (+91/‑107)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (38/39);  new base C (8/5);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAA  >  NC_000913/3780968‑3781255
                                                                                                                                             |                                                                                                                                                   
aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATCCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACGCttttt                                                                                                                                             <  1:211543/149‑1 (MQ=255)
aaGAGGTGCCTGAGGCACTGGCCCCCATGGCGAAAGGGAATGCGGCATCGTCGTTTGCGCCCGCTATCCGTACCCGTCTCCCTCCGGGAGTGGGGACCGTTCGGCGCTGTATG‑TCCTCCCTCACGCTGAAACGACACCGCACTCttttt                                                                                                                                             <  1:851617/149‑1 (MQ=255)
   tgTTTCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATACGCGTTTGCGCCCCCTCAGCCTAACCCTCTCCCTCAGGGAGAGGGGAACGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtc                                                                                                                                          <  2:174668/148‑1 (MQ=255)
    ggtGCCGGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                                                                         <  2:564301/147‑1 (MQ=255)
     ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCACTCACTCTGAAACGACACCGCACTCTTTTTTtatc                                                                                                                                        >  2:138601/1‑149 (MQ=255)
     tggCCTGCGGCACTGGCTCGCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCACGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                                                                                        <  1:538678/147‑1 (MQ=255)
     gtGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTAGG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                                                                                        <  1:543089/148‑1 (MQ=255)
      tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGCGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                                                                                       <  2:797983/149‑1 (MQ=255)
      tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                                                                                       <  2:154058/149‑1 (MQ=255)
           gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCGCCCTCGGGGAGAGGGGCCCGTTCGGGGCGGGCTG‑TACACCCGCACTCTGCAACGGCACCGGCCGCTTTTGTTCGCCCgcgc                                                                                                                                  >  1:668999/1‑149 (MQ=255)
           gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                                                                                  <  2:91650/149‑1 (MQ=255)
           gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCACTCCCTCAGGGAGAGCGGAACGATCGGCGCTGTATG‑TACAGCCGCCCTCTGAAACGACACCCCCCTCTATTTTACACCCACGc                                                                                                                                  >  2:65287/1‑149 (MQ=255)
                aCTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAACCTATCCCTCAGGGAGAGGGGTCCGGTCGGCGCTGTATG‑TACTGCATCACTCAGCACCGCCCCCGCACTCGTTTTTTCTCCCTCGcccgac                                                                                                                             >  2:56917/1‑146 (MQ=255)
                    gcacccATGGCGATAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTCACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGATACGACACCGCACTCTTTTTATCTCCCTCGCCCCTCCggg                                                                                                                         <  1:233548/146‑1 (MQ=255)
                    gCTGCCTTGGCGAACGGGAATGGGGCATAGTCGTTTGCGCGCCCTCTCCCTAACCCTCTCCCTCATGTAGAGGGGAGCGGTCGGCGCTGTAAG‑GACTCCCTCACTCTGTAACGAAACGGCTCTCTTTTTTTCTCCCTCGCCCCTCCgtg                                                                                                                         <  1:818653/149‑3 (MQ=255)
                    gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGAGGGGGCCGTTGGGGGGTGTTTG‑TTCTGCCGCACAGGGATACGCCCCCGCACTCTTTTTTTGGCCCGCGCCCCTACggg                                                                                                                         >  2:217121/1‑149 (MQ=255)
                       cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCTCACCCTAAGGGTGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGAACCCGCACTATTTTTTTATGCCTCGCCCCACCGGGgag                                                                                                                      >  2:425448/1‑149 (MQ=255)
                         cATGGCGAAAGGGAATGCGGCATAGTCGATTGCGCCCCCTCACCCGAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                                    <  1:563304/149‑1 (MQ=255)
                          aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGGTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCGGAGGGGCGAgg                                                                                                                   <  1:483687/149‑1 (MQ=255)
                          aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCAAACCCAAACCCACTCCCACAGGGAGGGGGGCCCGGTGGGCGCGGTTTG‑GACTGCCAGACTCAGAAACGGCACCCGACGGGTTTTTTAGACCAAGCCCCTCCggggcgggg                                                                                                                   >  2:805433/1‑144 (MQ=255)
                           tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                                                                                                                  <  2:622511/149‑1 (MQ=255)
                           tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TAATCCCTCACTCTGCAACGACACCGCCCGCTTTATTTATCCCTCGCCCCGCCGGGGAGAggg                                                                                                                  >  1:1006231/1‑149 (MQ=255)
                            ggagAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCACCCTAAAGGAGAGGGGCCCGTTCGGCCCTGTATG‑TACTCCCTCAAACTGAAACGACCCCGCACGATATTTTTCTCCCTCCCCACACCGGGGAGAGGGc                                                                                                                 >  2:621909/4‑149 (MQ=255)
                            ggCGAACGGGACGGCGGCAGAGGCGTTTGCGCCCCCTCACCCTAACCCACGCCCTCAGGGAGAGGGGACCGTTCGGCGCTGGATG‑TACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGc                                                                                                                 <  1:722676/149‑1 (MQ=255)
                              cGAAAGGGAATGCGGCTTAGTCGGTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCGTCACTCTGAAACGACACCGCACTGTTTTTATCGCCCACGCCCCTCCGGGGAGAGGGCCg                                                                                                               <  1:904046/149‑1 (MQ=255)
                                aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCGGGATG‑TACTCCCTCACGGTGAAACGGCACCGCACGCGTTTTTTCTCCCTCGCCCCGCCGGGGGGAGGGCCggg                                                                                                             >  1:13489/1‑149 (MQ=255)
                                aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                                             <  1:150385/149‑1 (MQ=255)
                                aaaGGCAAGGCGGCATAGACGTTGGCGCCCCCTCACCCACCCCCTCTCCCACAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                                             <  1:570438/149‑1 (MQ=255)
                                 aaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCAAACCCACTCCCACAGGGTGAGGGGGGCGGTGGGCGGGGTTTG‑GACGGCCTCACACATCAACGACACCGCCCACGTTTTTTATCCCTACACCCTCCCGGGGGcggggcgggg                                                                                                            >  2:827190/1‑149 (MQ=255)
                                    ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCCCGCTTTTTTTCTCCCACGCCCCGCCGGGGAGAGGGCCGGGGTgc                                                                                                         >  2:629606/1‑148 (MQ=255)
                                     gAATGCGGCATTGGCGTTTGCGCCCCCTCACCCTACCCCACTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                                                        <  1:942344/149‑1 (MQ=255)
                                      aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCACCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                       <  1:928844/149‑1 (MQ=255)
                                      aaTGCGGCATAGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCTGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                       <  2:621534/149‑1 (MQ=255)
                                           ggCATAGTCGTTTGCGCCCCCTCCCCCGAACCCTCTCCCTCAGGGAGAGGGGCCCGTTCGGCGCTGGATG‑TACTCCCTCACTCTGAAACGACACCGCAGTCGTTTTTTCTCCCTCGCCCCTCCGGGGGGAGGGCCCGGGGGAGGGGaaa                                                                                                  >  2:625183/1‑149 (MQ=255)
                                            gCATAGTCGTTTGCGCCCCCTCACCCTAACCCGCTCCCTCAGGGAGAGGGGACCGTGCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                                 <  2:345335/149‑1 (MQ=255)
                                              aTAGTCGTTGGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCGCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                                                                               <  2:388489/149‑1 (MQ=255)
                                                  tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCGCAGGGAGAGGGGACCGTTCGGGGCGGTATG‑TACTCCCTCCCGCTGAAACGCCGCCGCACTCGTTTTTTCACCCTCGCCCCGCCGGGGAGCGGGCCGGGGTGGGGGGTAAAGGCCGc                                                                                           >  2:252869/1‑149 (MQ=255)
                                                          gcCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCa                                                                                   >  2:26158/1‑149 (MQ=255)
                                                          ccccccTCCCGCTCAGCCACCCCGTCGGGGCGAGGGGACGGTTCGTCGCTGTATT‑TACTCCATCACTCTGAAAGCACACCGCATTCTTTTTTTCTCCCTCGCCCCTCCGGGGTTAGGGCCGGGGTTCGGGGAAAAGGCCGCACTGCTCa                                                                                   <  2:199192/148‑1 (MQ=255)
                                                           cccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGGACTGCTCAc                                                                                  >  1:24109/1‑149 (MQ=255)
                                                           cccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAc                                                                                  >  1:20439/1‑149 (MQ=255)
                                                                tCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtt                                                                             >  1:964262/1‑149 (MQ=255)
                                                                  aCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                                           >  2:305218/1‑149 (MQ=255)
                                                                  aCCCTAACCCTCACCCTCCGGGAGTGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGTAAAAGGCCGCACTGCTCACCAAtttt                                                                           <  2:964262/149‑1 (MQ=255)
                                                                   cccTCACCCTCTCCCCCCGGGAGAGGGGTCCGTGCGGCGCGGTATG‑TACACCCACCCTCTCAACCCACACCCCCCTCTTTTTTTCTCCCTCCCCCCGCCGGGGAGCGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                                          <  2:74884/149‑1 (MQ=255)
                                                                    ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                         >  1:963818/1‑149 (MQ=255)
                                                                    ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCGCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGCGGGCCGGGGTGAGGGGAAAACGCCGCACTGCTCACCCATTTTCt                                                                         >  1:871616/1‑149 (MQ=255)
                                                                    ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCCCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                         >  1:178209/1‑149 (MQ=255)
                                                                    ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGATGTATG‑TACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTATCTCCCTCGCCCATCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCCCGGCGCACCAATTTTCt                                                                         >  2:212064/1‑149 (MQ=255)
                                                                     cacaCCCGCCCCCTCCGGGAGCGGGGACCGGTCGGCGCTGTATG‑TACACCCTCACTCGGACACCCCCCCGCCCTCTTTTTTTCCCCCACCCCCCTCGGGGGAGAGGGCCGGGGGGCGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                        <  1:998289/146‑1 (MQ=255)
                                                                     caaccccTCTCCCTCCGGGTTTGGGGACCGTTCGGGGCTTTTTG‑TACTCCCTCCCACTGAAACGCCACCTCACTTTTTTTTTCTCCCTCGCCCCTCCGGGGAGTGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                        <  2:57967/145‑1 (MQ=255)
                                                                     caaaCCCTCCCCCACAGGCTGAGGGGACCGTCCGCCGCTGTTTG‑ACCTCCCACACTCTGAACCGACCCCCCACTCTTTTTTTCACCCTCGCCCCGCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                        <  2:949679/147‑1 (MQ=255)
                                                                      tAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGTAAAGGCCGCACTGCTCCCCAATTTTCTGc                                                                       >  2:237421/1‑149 (MQ=255)
                                                                         cccTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATC‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTAt                                                                    <  1:731480/149‑1 (MQ=255)
                                                                          ccTCTCCCTCAGGGTGTGGGGACCGTTCGGCGCTGTATG‑AACTCCCTCACTCTTAAACTACACCGAACTCTTTTTTTATCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGTTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                   <  2:20439/149‑1 (MQ=255)
                                                                          ccTCTCCCTCAGGGTGAGGGGTCGTTTCGGGGCTGTATG‑TACTCCCTAAATCTTAAACGACACCGAAATCTTTTTTTCTCCCTCGCCCCTCCGGGGAGTGGGCCGGGGTGTGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                   <  2:24109/149‑1 (MQ=255)
                                                                          ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACCCCGCCCTCTGTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGCGGATCGGCCGCACTGGTCACCAAGTTTCTGCTATc                                                                   >  1:409990/1‑149 (MQ=255)
                                                                           ctctCGCTCAGGGAGAGGGGACCGTTCGGCGCTGTAGG‑TACTCCCTCACTCTGAAACGCGCCCGCCCTGTTTTTTTCTACCTCGCCCCTCCGGGGCGAGGGCCGGGGTGAGGGGAAAGGGCCGCCCTGCTCACCAAGTTTCTGCTATcc                                                                  >  1:114076/1‑149 (MQ=255)
                                                                             cacccTCCGGGCGTGGGGACCGTTGGGGGCTGTATT‑TACTCCCTCACTTTTAAACCACACCGCACTATTTTTTTCTCCCACGCCCCTCCGGGGAGTGGGCCGGGGTGCGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                                <  2:453361/147‑1 (MQ=255)
                                                                                    ctggtgagGGGACCGGTCGGCGCTG‑ATGCTCCTCCCTCTCCCTTACACGCCCCCGCACTGTTTTTAGCTCCCGCGCCCCCCCGGGGTGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                                         <  1:953428/144‑1 (MQ=255)
                                                                                      ggAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGCGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                                       >  1:471016/1‑149 (MQ=255)
                                                                                            gtgACCGTTCGGCGCTGTATG‑AACTCCATAACTTTCAAACTACACCGCACATTTTTTTACTCCCTCGCCCCCCCGGGGAGTGGGCCGGGGTTAGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                                 <  1:138601/147‑1 (MQ=255)
                                                                                            gggACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGCCACCGCACTCTTATTTTCTCCCTCGCCCCTCCGGGGATAGGGCCGTGGGGAGGGGTAAAGGCCGCAATGCTCACCAATTTTGTGCTATCCTGCACCCGCACAAAggg                                                 >  2:668836/1‑149 (MQ=255)
                                                                                            gggACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                                 >  1:217526/1‑149 (MQ=255)
                                                                                            gggACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                                 >  1:210227/1‑149 (MQ=255)
                                                                                             ggACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGAAACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGTGAAAACGCCGCACTCATCGCCAATTTTCTGATATCCTGCCCCACCACTAAgggg                                                >  1:283494/1‑149 (MQ=255)
                                                                                             ggACCGTGCGGCGCTGTATG‑AACTCCCTCCCTCTGAAACGACCCCGCACCCTTTTTTTCTCCCGCGCCCCGCCGGGGCGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                                <  1:631354/149‑1 (MQ=255)
                                                                                              gACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTATggggg                                               >  2:69355/1‑149 (MQ=255)
                                                                                              gACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                               >  2:870328/1‑149 (MQ=255)
                                                                                               aCCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                              >  2:657219/1‑149 (MQ=255)
                                                                                                ccGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTTCTCACAAATTTTCTGCTATCCTGCCCCCGCACGAAGGGGGCa                                             >  2:957499/1‑149 (MQ=255)
                                                                                                 cGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                            >  2:740267/1‑149 (MQ=255)
                                                                                                     cggcggGGTATG‑TCCTCCCTCCCACTCAAACCCCCCCGCACACTTTTTTTCTCCCGCCCCCCGCCGGGGAGCGGGCCGGGGTGAGGGGACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATg                                        <  2:994136/149‑1 (MQ=255)
                                                                                                       gcgcTGTATG‑TACTCCCTCACTCTGAAACGACACCGCCCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                                      <  2:178209/149‑1 (MQ=255)
                                                                                                         gcTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCACTCCGGGGAGAGGGCCGGGTTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTATGCTATCCTGCCCACGCAGTAAGGGGGCAGTATGCTaa                                    >  2:838419/1‑149 (MQ=255)
                                                                                                         gcTGTAAG‑TACTCCCTCACGCTGAAACGACACCGCACTCGTTATTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGCGTGAGGGGAACCGGCCGCCCTCCTCCCCAATTTTCTGCTATCCTGCCCCCGCTCTTAGGGGGCAGGATGCTaa                                    >  2:808152/1‑149 (MQ=255)
                                                                                                          cTGTATG‑TACTCCCTCACTCTCAACCCCCCCCCCACTCTTTTTTTCCCCCTCCCCCCCCGGGGGAGTGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                                   <  1:981647/149‑1 (MQ=255)
                                                                                                                  tccccccTCACTCTGAAACGACACCCCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa                            <  1:442951/145‑1 (MQ=255)
                                                                                                                    ccccccCACTCTGAACGGCCACCGCACTCTTTTTTGCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGTTGTGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt                          <  2:975393/147‑1 (MQ=255)
                                                                                                                     tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACtt                         >  1:620128/1‑149 (MQ=255)
                                                                                                                     tCCCTCACTCTGAAACGACACCGCACTCGTTTTTTCTCACTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTACCCAATTTTCTGATATCATGCCCACGCACTAAGGTGGCAGTATGCTATACATCGGACtt                         >  2:923228/1‑149 (MQ=255)
                                                                                                                       cccccctctGAAAGTCCCCCGCACTCTTTTTTTCTCCCGCGCCCCCCGGGGTTGTGGGCCGGGGGGCGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa                       <  2:529792/144‑1 (MQ=255)
                                                                                                                       ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCTGGGAGAGGCCCGGGGTGAGGGGAAAAGGCCGCACTCCTCCCCAAGTTTCTGCTATCCTGCCCCCGCACGAAGGGGGCAGTATGCTAAACATCGTACttaa                       >  1:120294/1‑147 (MQ=255)
                                                                                                                        cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGTAAAGGCCGCACTGCTGACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAGGCTAAACCTGGTACTTTAc                      >  2:232239/1‑149 (MQ=255)
                                                                                                                         agaCTCTGAAACGCCACCGCCCTGTTTTTTGCACCCTCGCCCCTCCGGGGCGAGGGCGGGGGTTCGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg                     <  1:957499/147‑1 (MQ=255)
                                                                                                                            ctctAAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAc                  >  2:397947/1‑149 (MQ=255)
                                                                                                                               tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGTGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg               >  1:305853/1‑149 (MQ=255)
                                                                                                                                  aaCGACACCCCACTCTTTTTTTCTCCCTCGCCCCTCAGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACATTACGAACCAGaat            >  2:119099/1‑148 (MQ=255)
                                                                                                                                    cGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAtt          >  1:880562/1‑149 (MQ=255)
                                                                                                                                      aaacCGCCCTATTTTTAACTCCCACCCCCCGCCGGGGTGAGGGCGGGGGTCCGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTcc        <  1:421915/147‑1 (MQ=255)
                                                                                                                                           gCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGTGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCaa   >  1:428426/1‑149 (MQ=255)
                                                                                                                                            cACTCTTTTTATCTCCCACGCCCCGCCGGGGAGATGACCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAATTTCAGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCaaa  <  1:102590/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
AAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAA  >  NC_000913/3780968‑3781255

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