Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,101 A→C 18.6% intergenic (+84/‑114) 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,1010AC18.6% 127.2 / 2.5 90intergenic (+84/‑114)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base C (8/10);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATT  >  NC_000913/3780953‑3781247
                                                                                                                                                     |                                                                                                                                                  
tGCAGGGGCCGTGTAAAGAGTTGCCTGCGGCAATGGCTCCCATGGCGAAAGGGAATGCGGAATATTCTTTTGCGCCCCCTCAACTTAACCCTCTCCGTCAGGGAGAGGTGACCGTTCGGAGCTGTATG‑TACTCCCTCACTCTGAAACGa                                                                                                                                                    <  2:647611/149‑1 (MQ=255)
ggCAGGGGCTGGGTAAAGAGTGGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGa                                                                                                                                                    <  2:996758/148‑1 (MQ=255)
 gCAGGGGCTGGGTAAAGAGTTGCCTGCGGAACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCCTCGCTGTATG‑TACTCCCTCACTCTGAACCGac                                                                                                                                                   >  2:602701/1‑149 (MQ=255)
  cAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGaca                                                                                                                                                  <  1:827034/149‑1 (MQ=255)
   aGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGacac                                                                                                                                                 <  1:775272/149‑1 (MQ=255)
    ggggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTGGGCGCTGTATG‑TACTCCCTCACTCTGAAACGCCAcc                                                                                                                                                >  2:442951/1‑149 (MQ=255)
               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:210227/1‑149 (MQ=255)
                                                                                                           gggACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                         >  1:217526/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)
                                                                                                                                                     |                                                                                                                                                  
TGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATG‑TACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATT  >  NC_000913/3780953‑3781247

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: ‑