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

AAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAA  >  NC_000913/3780968‑3781255
                                                                                                                                             |                                                                                                                                                   
aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATCCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACGCttttt                                                                                                                                             <  1:211543/149‑1 (MQ=255)
aaGAGGTGCCTGAGGCACTGGCCCCCATGGCGAAAGGGAATGCGGCATCGTCGTTTGCGCCCGCTATCCGTACCCGTCTCCCTCCGGGAGTGGGGACCGTTCGGCGCTGTATGTCCTCCCTCACGCTGAAACGACACCGCACTCttttt                                                                                                                                             <  1:851617/149‑1 (MQ=255)
   tgTTTCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATACGCGTTTGCGCCCCCTCAGCCTAACCCTCTCCCTCAGGGAGAGGGGAACGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtc                                                                                                                                          <  2:174668/148‑1 (MQ=255)
    ggtGCCGGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                                                                         <  2:564301/147‑1 (MQ=255)
     ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCACTCACTCTGAAACGACACCGCACTCTTTTTTtatc                                                                                                                                        >  2:138601/1‑149 (MQ=255)
     tggCCTGCGGCACTGGCTCGCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCACGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                                                                                        <  1:538678/147‑1 (MQ=255)
     gtGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTAGGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                                                                                        <  1:543089/148‑1 (MQ=255)
      tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGCGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                                                                                       <  2:797983/149‑1 (MQ=255)
      tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                                                                                       <  2:154058/149‑1 (MQ=255)
           gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCGCCCTCGGGGAGAGGGGCCCGTTCGGGGCGGGCTGTACACCCGCACTCTGCAACGGCACCGGCCGCTTTTGTTCGCCCgcgc                                                                                                                                  >  1:668999/1‑149 (MQ=255)
           gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                                                                                  <  2:91650/149‑1 (MQ=255)
           gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCACTCCCTCAGGGAGAGCGGAACGATCGGCGCTGTATGTACAGCCGCCCTCTGAAACGACACCCCCCTCTATTTTACACCCACGc                                                                                                                                  >  2:65287/1‑149 (MQ=255)
                aCTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAACCTATCCCTCAGGGAGAGGGGTCCGGTCGGCGCTGTATGTACTGCATCACTCAGCACCGCCCCCGCACTCGTTTTTTCTCCCTCGcccgac                                                                                                                             >  2:56917/1‑146 (MQ=255)
                    gcacccATGGCGATAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTCACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGATACGACACCGCACTCTTTTTATCTCCCTCGCCCCTCCggg                                                                                                                         <  1:233548/146‑1 (MQ=255)
                    gCTGCCTTGGCGAACGGGAATGGGGCATAGTCGTTTGCGCGCCCTCTCCCTAACCCTCTCCCTCATGTAGAGGGGAGCGGTCGGCGCTGTAAGGACTCCCTCACTCTGTAACGAAACGGCTCTCTTTTTTTCTCCCTCGCCCCTCCgtg                                                                                                                         <  1:818653/149‑3 (MQ=255)
                    gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGAGGGGGCCGTTGGGGGGTGTTTGTTCTGCCGCACAGGGATACGCCCCCGCACTCTTTTTTTGGCCCGCGCCCCTACggg                                                                                                                         >  2:217121/1‑149 (MQ=255)
                       cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCTCACCCTAAGGGTGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGAACCCGCACTATTTTTTTATGCCTCGCCCCACCGGGgag                                                                                                                      >  2:425448/1‑149 (MQ=255)
                         cATGGCGAAAGGGAATGCGGCATAGTCGATTGCGCCCCCTCACCCGAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                                    <  1:563304/149‑1 (MQ=255)
                          aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGGTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCGGAGGGGCGAgg                                                                                                                   <  1:483687/149‑1 (MQ=255)
                          aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCAAACCCAAACCCACTCCCACAGGGAGGGGGGCCCGGTGGGCGCGGTTTGGACTGCCAGACTCAGAAACGGCACCCGACGGGTTTTTTAGACCAAGCCCCTCCggggcgggg                                                                                                                   >  2:805433/1‑144 (MQ=255)
                           tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                                                                                                                  <  2:622511/149‑1 (MQ=255)
                           tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTAATCCCTCACTCTGCAACGACACCGCCCGCTTTATTTATCCCTCGCCCCGCCGGGGAGAggg                                                                                                                  >  1:1006231/1‑149 (MQ=255)
                            ggagAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCACCCTAAAGGAGAGGGGCCCGTTCGGCCCTGTATGTACTCCCTCAAACTGAAACGACCCCGCACGATATTTTTCTCCCTCCCCACACCGGGGAGAGGGc                                                                                                                 >  2:621909/4‑149 (MQ=255)
                            ggCGAACGGGACGGCGGCAGAGGCGTTTGCGCCCCCTCACCCTAACCCACGCCCTCAGGGAGAGGGGACCGTTCGGCGCTGGATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGc                                                                                                                 <  1:722676/149‑1 (MQ=255)
                              cGAAAGGGAATGCGGCTTAGTCGGTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCGTCACTCTGAAACGACACCGCACTGTTTTTATCGCCCACGCCCCTCCGGGGAGAGGGCCg                                                                                                               <  1:904046/149‑1 (MQ=255)
                                aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCGGGATGTACTCCCTCACGGTGAAACGGCACCGCACGCGTTTTTTCTCCCTCGCCCCGCCGGGGGGAGGGCCggg                                                                                                             >  1:13489/1‑149 (MQ=255)
                                aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                                             <  1:150385/149‑1 (MQ=255)
                                aaaGGCAAGGCGGCATAGACGTTGGCGCCCCCTCACCCACCCCCTCTCCCACAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                                             <  1:570438/149‑1 (MQ=255)
                                 aaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCAAACCCACTCCCACAGGGTGAGGGGGGCGGTGGGCGGGGTTTGGACGGCCTCACACATCAACGACACCGCCCACGTTTTTTATCCCTACACCCTCCCGGGGGcggggcgggg                                                                                                            >  2:827190/1‑149 (MQ=255)
                                    ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCCCGCTTTTTTTCTCCCACGCCCCGCCGGGGAGAGGGCCGGGGTgc                                                                                                         >  2:629606/1‑148 (MQ=255)
                                     gAATGCGGCATTGGCGTTTGCGCCCCCTCACCCTACCCCACTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                                                        <  1:942344/149‑1 (MQ=255)
                                      aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCACCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                       <  1:928844/149‑1 (MQ=255)
                                      aaTGCGGCATAGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCTGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                       <  2:621534/149‑1 (MQ=255)
                                           ggCATAGTCGTTTGCGCCCCCTCCCCCGAACCCTCTCCCTCAGGGAGAGGGGCCCGTTCGGCGCTGGATGTACTCCCTCACTCTGAAACGACACCGCAGTCGTTTTTTCTCCCTCGCCCCTCCGGGGGGAGGGCCCGGGGGAGGGGaaa                                                                                                  >  2:625183/1‑149 (MQ=255)
                                            gCATAGTCGTTTGCGCCCCCTCACCCTAACCCGCTCCCTCAGGGAGAGGGGACCGTGCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                                 <  2:345335/149‑1 (MQ=255)
                                              aTAGTCGTTGGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCGCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                                                                               <  2:388489/149‑1 (MQ=255)
                                                  tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCGCAGGGAGAGGGGACCGTTCGGGGCGGTATGTACTCCCTCCCGCTGAAACGCCGCCGCACTCGTTTTTTCACCCTCGCCCCGCCGGGGAGCGGGCCGGGGTGGGGGGTAAAGGCCGc                                                                                           >  2:252869/1‑149 (MQ=255)
                                                          gcCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCa                                                                                   >  2:26158/1‑149 (MQ=255)
                                                          ccccccTCCCGCTCAGCCACCCCGTCGGGGCGAGGGGACGGTTCGTCGCTGTATTTACTCCATCACTCTGAAAGCACACCGCATTCTTTTTTTCTCCCTCGCCCCTCCGGGGTTAGGGCCGGGGTTCGGGGAAAAGGCCGCACTGCTCa                                                                                   <  2:199192/148‑1 (MQ=255)
                                                           cccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGGACTGCTCAc                                                                                  >  1:24109/1‑149 (MQ=255)
                                                           cccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAc                                                                                  >  1:20439/1‑149 (MQ=255)
                                                                tCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtt                                                                             >  1:964262/1‑149 (MQ=255)
                                                                  aCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                                           >  2:305218/1‑149 (MQ=255)
                                                                  aCCCTAACCCTCACCCTCCGGGAGTGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGTAAAAGGCCGCACTGCTCACCAAtttt                                                                           <  2:964262/149‑1 (MQ=255)
                                                                   cccTCACCCTCTCCCCCCGGGAGAGGGGTCCGTGCGGCGCGGTATGTACACCCACCCTCTCAACCCACACCCCCCTCTTTTTTTCTCCCTCCCCCCGCCGGGGAGCGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                                          <  2:74884/149‑1 (MQ=255)
                                                                    ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                         >  1:963818/1‑149 (MQ=255)
                                                                    ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCGCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGCGGGCCGGGGTGAGGGGAAAACGCCGCACTGCTCACCCATTTTCt                                                                         >  1:871616/1‑149 (MQ=255)
                                                                    ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCCCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                         >  1:178209/1‑149 (MQ=255)
                                                                    ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGATGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTATCTCCCTCGCCCATCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCCCGGCGCACCAATTTTCt                                                                         >  2:212064/1‑149 (MQ=255)
                                                                     cacaCCCGCCCCCTCCGGGAGCGGGGACCGGTCGGCGCTGTATGTACACCCTCACTCGGACACCCCCCCGCCCTCTTTTTTTCCCCCACCCCCCTCGGGGGAGAGGGCCGGGGGGCGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                        <  1:998289/146‑1 (MQ=255)
                                                                     caaccccTCTCCCTCCGGGTTTGGGGACCGTTCGGGGCTTTTTGTACTCCCTCCCACTGAAACGCCACCTCACTTTTTTTTTCTCCCTCGCCCCTCCGGGGAGTGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                        <  2:57967/145‑1 (MQ=255)
                                                                     caaaCCCTCCCCCACAGGCTGAGGGGACCGTCCGCCGCTGTTTGACCTCCCACACTCTGAACCGACCCCCCACTCTTTTTTTCACCCTCGCCCCGCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                        <  2:949679/147‑1 (MQ=255)
                                                                      tAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGTAAAGGCCGCACTGCTCCCCAATTTTCTGc                                                                       >  2:237421/1‑149 (MQ=255)
                                                                         cccTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATCTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTAt                                                                    <  1:731480/149‑1 (MQ=255)
                                                                          ccTCTCCCTCAGGGTGTGGGGACCGTTCGGCGCTGTATGAACTCCCTCACTCTTAAACTACACCGAACTCTTTTTTTATCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGTTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                   <  2:20439/149‑1 (MQ=255)
                                                                          ccTCTCCCTCAGGGTGAGGGGTCGTTTCGGGGCTGTATGTACTCCCTAAATCTTAAACGACACCGAAATCTTTTTTTCTCCCTCGCCCCTCCGGGGAGTGGGCCGGGGTGTGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                   <  2:24109/149‑1 (MQ=255)
                                                                          ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACCCCGCCCTCTGTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGCGGATCGGCCGCACTGGTCACCAAGTTTCTGCTATc                                                                   >  1:409990/1‑149 (MQ=255)
                                                                           ctctCGCTCAGGGAGAGGGGACCGTTCGGCGCTGTAGGTACTCCCTCACTCTGAAACGCGCCCGCCCTGTTTTTTTCTACCTCGCCCCTCCGGGGCGAGGGCCGGGGTGAGGGGAAAGGGCCGCCCTGCTCACCAAGTTTCTGCTATcc                                                                  >  1:114076/1‑149 (MQ=255)
                                                                             cacccTCCGGGCGTGGGGACCGTTGGGGGCTGTATTTACTCCCTCACTTTTAAACCACACCGCACTATTTTTTTCTCCCACGCCCCTCCGGGGAGTGGGCCGGGGTGCGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                                <  2:453361/147‑1 (MQ=255)
                                                                                    ctggtgagGGGACCGGTCGGCGCTGATGCTCCTCCCTCTCCCTTACACGCCCCCGCACTGTTTTTAGCTCCCGCGCCCCCCCGGGGTGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                                         <  1:953428/144‑1 (MQ=255)
                                                                                      ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGCGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                                       >  1:471016/1‑149 (MQ=255)
                                                                                            gtgACCGTTCGGCGCTGTATGAACTCCATAACTTTCAAACTACACCGCACATTTTTTTACTCCCTCGCCCCCCCGGGGAGTGGGCCGGGGTTAGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                                 <  1:138601/147‑1 (MQ=255)
                                                                                            gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTATTTTCTCCCTCGCCCCTCCGGGGATAGGGCCGTGGGGAGGGGTAAAGGCCGCAATGCTCACCAATTTTGTGCTATCCTGCACCCGCACAAAggg                                                 >  2:668836/1‑149 (MQ=255)
                                                                                            gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                                 >  1:217526/1‑149 (MQ=255)
                                                                                            gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                                 >  1:210227/1‑149 (MQ=255)
                                                                                             ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGAAACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGTGAAAACGCCGCACTCATCGCCAATTTTCTGATATCCTGCCCCACCACTAAgggg                                                >  1:283494/1‑149 (MQ=255)
                                                                                             ggACCGTGCGGCGCTGTATGAACTCCCTCCCTCTGAAACGACCCCGCACCCTTTTTTTCTCCCGCGCCCCGCCGGGGCGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                                <  1:631354/149‑1 (MQ=255)
                                                                                              gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTATggggg                                               >  2:69355/1‑149 (MQ=255)
                                                                                              gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                               >  2:870328/1‑149 (MQ=255)
                                                                                               aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                              >  2:657219/1‑149 (MQ=255)
                                                                                                ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTTCTCACAAATTTTCTGCTATCCTGCCCCCGCACGAAGGGGGCa                                             >  2:957499/1‑149 (MQ=255)
                                                                                                 cGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                            >  2:740267/1‑149 (MQ=255)
                                                                                                     cggcggGGTATGTCCTCCCTCCCACTCAAACCCCCCCGCACACTTTTTTTCTCCCGCCCCCCGCCGGGGAGCGGGCCGGGGTGAGGGGACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATg                                        <  2:994136/149‑1 (MQ=255)
                                                                                                       gcgcTGTATGTACTCCCTCACTCTGAAACGACACCGCCCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                                      <  2:178209/149‑1 (MQ=255)
                                                                                                         gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCACTCCGGGGAGAGGGCCGGGTTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTATGCTATCCTGCCCACGCAGTAAGGGGGCAGTATGCTaa                                    >  2:838419/1‑149 (MQ=255)
                                                                                                         gcTGTAAGTACTCCCTCACGCTGAAACGACACCGCACTCGTTATTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGCGTGAGGGGAACCGGCCGCCCTCCTCCCCAATTTTCTGCTATCCTGCCCCCGCTCTTAGGGGGCAGGATGCTaa                                    >  2:808152/1‑149 (MQ=255)
                                                                                                          cTGTATGTACTCCCTCACTCTCAACCCCCCCCCCACTCTTTTTTTCCCCCTCCCCCCCCGGGGGAGTGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                                   <  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)
                                                                                                                                             |                                                                                                                                                   
AAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAA  >  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: