Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,103 A→C 13.9% intergenic (+86/‑112) 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,1030AC13.9% 167.8 / 2.7 90intergenic (+86/‑112)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (45/30);  new base C (5/9);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCG  >  NC_000913/3780955‑3781250
                                                                                                                                                    |                                                                                                                                                   
cAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGacac                                                                                                                                                    <  2:1141089/150‑1 (MQ=255)
cAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCGAACCCACTACCTCAAGGAGAGGGGACCGTTCGGCGCGGTATGTACAACCTAACACTGAAACCacac                                                                                                                                                    >  2:1185415/1‑150 (MQ=255)
   ggtcTGTGTAATGCGTTGCCTGCGTCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCGTCACCCTATCCCGCTCCCTCTGCGTGAGGGGACCGTTCGGCGCTGTATGTACTGCCTCACTCTTAAACGACACCGc                                                                                                                                                 <  2:1458766/147‑1 (MQ=255)
   gggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCGCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGTGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGc                                                                                                                                                 <  1:433968/150‑1 (MQ=255)
   gggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGGGAGGGGACCGTTAGGGGCTGTAGGGACTCCCTCACTCTGAAACGCCACCGc                                                                                                                                                 >  1:1022936/1‑150 (MQ=255)
   gggCTGGGTAAAGAGTGGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGc                                                                                                                                                 <  2:1458762/150‑1 (MQ=255)
    ggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCAAACCCAAACCCTCACCCACAGGGAGAGGGGACCGTTAGGCGCAGTATGTACTCCCTCACTCTTAAACGACACCCCa                                                                                                                                                >  2:327993/1‑150 (MQ=255)
          gTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                                                                                          <  2:1341128/150‑1 (MQ=255)
           tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttt                                                                                                                                         <  2:933690/150‑1 (MQ=255)
               gagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtc                                                                                                                                     <  1:1462079/150‑1 (MQ=255)
                 gTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                                                                                   <  1:57664/150‑1 (MQ=255)
                   tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATTCGGCATAGTCGTATGCGCCCACTAACAATAACCCTCACCCTCAGTGAGAGCGGACCGTTCGGCGCTGTAAGTAATCCCACCCTATGCAACGACCCCGCACTCTTTTTTTATccc                                                                                                                                 >  1:258020/1‑150 (MQ=255)
                   tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATTGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGTGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTccc                                                                                                                                 <  1:839356/150‑1 (MQ=255)
                     ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                                                                                               <  1:538591/150‑1 (MQ=255)
                        gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcc                                                                                                                            <  2:1654354/150‑1 (MQ=255)
                              cTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCg                                                                                                                      >  1:1096415/1‑150 (MQ=255)
                                 gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgggg                                                                                                                   <  1:1829184/150‑1 (MQ=255)
                                 gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgggg                                                                                                                   <  1:1339704/150‑1 (MQ=255)
                                 gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCTAACCCAAACCCGGAGGGCGGGGGGGCCGTTAGGGGCGGTTTGTAATCCCTCACGCGGAAACCGCCCCGCCCGCTTTGTTTCGGCCTCGCCCCACCcggg                                                                                                                   >  1:737914/1‑150 (MQ=255)
                                  cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTCACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTAATCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                                                                                  <  1:84504/150‑1 (MQ=255)
                                     ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGGCCGTTCGGCGCTGGATGTACACCCTCACTCTGAAACGGAACCGCACTCATTTTTACTACATCGCCCCTCCGGGGagag                                                                                                               >  2:76380/1‑150 (MQ=255)
                                          gCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCTCTCCCGCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACCCTGAACCGACACCCCCCTCTTGTTTTCTCCCGCGCCCCTCCGGGgggcggggcg                                                                                                          >  2:291537/1‑150 (MQ=255)
                                           cGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCTAACCCACAACCGCAGGGAGAGGGGACCCTTGGGCGCTGTATGTACAACCTCACTCTGGAACGACACCGCCCACATTTTTTATACCACCCCCCGCCGGGGGGAGGGGCgg                                                                                                         >  2:454764/1‑150 (MQ=255)
                                           cGAAAGGGAATGCGGCATAGTCGCTTGCGCCCCCACACCCTAACCCTCTCCCTCCGGGAGAGGGGACCGTCCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                         <  1:156963/150‑1 (MQ=255)
                                            gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                                        <  1:589719/150‑1 (MQ=255)
                                               aGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCCCCGCACACTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                                                                     >  2:61054/1‑150 (MQ=255)
                                                 ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGTGAGGGGACCGGTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCCCTCTTTTTTTCTCCCTCGCCCCGCCGGGGAGAGGGGGCGGGTGAg                                                                                                   >  2:700657/1‑150 (MQ=255)
                                                 ggAATGCGGCATAGTCGTTTGCGCCCCCTAACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATTTACTCCCTCACTCAGAAAAGCCACCGGAAGCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                                                   >  2:365513/1‑150 (MQ=255)
                                                   aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACGCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                                                                 <  1:543416/150‑1 (MQ=255)
                                                   aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCCCAGGGAGAGGGGACCGTTCGGCGCTGTATGGACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                                                                 <  1:1651083/150‑1 (MQ=255)
                                                   aaTGCGGCATAGGCGTTTGCGCCCCCTCAGCCACACCCTCACCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                                                                 <  1:1459183/150‑1 (MQ=255)
                                                       cGGCATAGTCGTTTGCGCCCCCACACCCAAACCCACACCCAAAGGGAGCGGGGGCCGTGCGGCGCAGGATGTACTACCACACACGGAAACCCCACCGCCAGCTTTTTTTATCCCTAGCCCCTCCGGGGAGAGGGCCGGGGGGGGGGGaaa                                                                                             >  2:1754516/1‑150 (MQ=255)
                                                        ggCATAGTCGTTTGCGCCCCCGCCCCCAAACCCACTCCCTCAGGGAGCGGGGGCCGGTCGGCGCGGTTTGTTCTGCCACCCACTGCAACGCCCCCCCCCTCGTTTTTTCTCCCGCGCCCCGCCGGGGCGCGGGCCGGGGGGAGGGGgaaa                                                                                            >  2:698656/1‑150 (MQ=255)
                                                                cGTTTGCGCCCCCACACCCTAACCCTCACCCTCAGGGAGAGGGGACCGTTAGGCGCTGTATGTACTCCCGCCGTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCGCCGGGCAGACGGCCGGGGTGAGGGGAGAAGGCCGCAc                                                                                    >  2:1283872/1‑150 (MQ=255)
                                                                    tGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGCACCGACCCCGCACTCTTTTTTTCTCCCTCACCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCt                                                                                >  1:1218001/1‑150 (MQ=255)
                                                                      cgcCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCGCACTCTGAAACGACACCGCACTCTTTTTTTCTACCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCGCTGCTcc                                                                              >  1:1612300/1‑149 (MQ=255)
                                                                       gcCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCCCTGCGCAc                                                                             >  2:894160/1‑150 (MQ=255)
                                                                        cccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGGCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAcc                                                                            >  1:216143/1‑150 (MQ=255)
                                                                          cccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCCCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCaa                                                                          >  1:889725/1‑150 (MQ=255)
                                                                          cccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAATGGCCGCCGTGCGCACCac                                                                          >  2:95328/1‑149 (MQ=255)
                                                                              cACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                                      >  1:787517/1‑150 (MQ=255)
                                                                               ccccaaccccTCTCCCTCCGTGAGTGGGGACGGGTCGGCGCTGTTTGTACTCCCTCCCTCTCAAACGACACGGCAGTCTTTTTTTCCCCCTGCCCCCTCGGGGGCGTGGGCGGGGGTGCGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                                     <  2:431163/143‑1 (MQ=255)
                                                                               ccccTACCCCTCTCCCTCAGGGCTTGGGGACCGTTCGGCGCTGAATGTACTCCCTCCCTCTGACCCGACCCCCCCCTCTTTTTTTCTCCCTCGCCCCTCCGGGGTGTGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                                     <  2:1051762/149‑1 (MQ=255)
                                                                               aCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                                     >  2:1174144/1‑150 (MQ=255)
                                                                                cccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTACCTCACTCTGAAACGACGCCGCAATCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCCCCAATGTTcg                                                                    >  2:283577/1‑149 (MQ=255)
                                                                                cccAAACCCTCCCCCATCGGGTGAGGGGACCGACCGTCGCTGTATGTACTCCAACACTCAGAACCCACCCCCCAATTTTTTTTTCTCCCCCGCCCCTCCGGGGAGGGGGCCGGGGTGTGGGTAAAAGGCCGCACTGCTCACCAATTTTCt                                                                    <  2:1661518/150‑1 (MQ=255)
                                                                                 ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                   <  2:1288187/150‑1 (MQ=255)
                                                                                  aaaaCCCTCTCCCTCTGGGAGAGGGGACGTTTCGGCGCTGTATGTACTCCCTCACGCTGAAACGACACCGCACACTTTTTTTCGCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                                  <  2:121538/148‑1 (MQ=255)
                                                                                       ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCGCACCAATTTTCTGCTATcc                                                             >  1:834676/1‑150 (MQ=255)
                                                                                         tctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                           >  2:1284341/1‑150 (MQ=255)
                                                                                           tcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                                         >  1:590415/1‑150 (MQ=255)
                                                                                           tcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                                         <  1:1081016/150‑1 (MQ=255)
                                                                                              cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAACCGACACCGCCCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                                      >  1:1097/1‑150 (MQ=255)
                                                                                              cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                                      >  2:5372/1‑150 (MQ=255)
                                                                                               tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                                     >  1:373649/1‑150 (MQ=255)
                                                                                                  gggtgcggggACCGTCCGGCGGTGTATGACCACCCTCCCTCTCAACCCACCCCGCACGCTTTTTTTCCCCCTCGCCCCTCCGGGGCGCGGGCGGGGGGGAGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                                  <  1:132168/144‑1 (MQ=255)
                                                                                                   ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCACCGCGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTCCTATCCTGCCCCCGCACt                                                 >  2:733972/1‑150 (MQ=255)
                                                                                                    gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCGCACTaa                                               >  2:870125/1‑150 (MQ=255)
                                                                                                       tggggACCGTGCGGCGCTTTATTTACTCCCTCACTCTGAACCGACACCGCACTATTTTTTTCTCCCTCGCCCCTCGGGGGTGAGGGCCGGGGTGTGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                             <  1:1890706/149‑1 (MQ=255)
                                                                                                        ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCGCTAAggg                                            >  2:475305/1‑150 (MQ=255)
                                                                                                        ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                            >  1:1752339/1‑150 (MQ=255)
                                                                                                        ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                            >  1:9216/1‑150 (MQ=255)
                                                                                                        ggggACCGTTCGGCGCTGTATGTACACCCTCACTCTGAAACGACACCGCACCCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTTAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                            <  2:64068/150‑1 (MQ=255)
                                                                                                           gcccGTGCGGCGCTGTTTTTACCCCCTCACTCTCAACCCCCCCCCCACTTTTTTTTCCTCCCCCCCCCCGCGGGGGTGTGGGCGGGGGTGCGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                         <  2:414710/148‑1 (MQ=255)
                                                                                                           gACGGTTCGGCGCTGTATGTACTCCCTCACTCTGAACCGACCCCGCACTTTTTTTTTCTCCGTCCCCCCTCCGGGGTGTGGGCCGGGGTGCGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                         <  1:327993/150‑1 (MQ=255)
                                                                                                           gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                         <  2:235492/150‑1 (MQ=255)
                                                                                                           gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGATATCCTGCCCCCGCACTAAGGGGGc                                         >  2:1792215/1‑150 (MQ=255)
                                                                                                             ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAACCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCCATTTTCTGCTATCCTGCCCCCGCCCGAAGGGGGCAg                                       >  1:963736/1‑150 (MQ=255)
                                                                                                             ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                       >  1:1610236/1‑150 (MQ=255)
                                                                                                             ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAGGGGGCAGt                                      >  2:1095010/1‑150 (MQ=255)
                                                                                                              cGTGCGGCGGTTTTTGAACCCCCTCACTCACAAACCCCCCCGCACTTTTTTTTTCTCCCCCGCCCCGCCGGGGCGGGGGCGGGGGTGCGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                                      <  2:483566/150‑1 (MQ=255)
                                                                                                                ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                    >  1:771762/1‑150 (MQ=255)
                                                                                                                  gggCGCTGTATGTACTCCCGCACTCGCAACCGACACCGCCCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                                  <  1:77872/149‑1 (MQ=255)
                                                                                                                  cGGCGCTGTTTGTACTCCCTCACTCTCAACCGACCCCCCACTCTTTTTTTCTCCCGCGCCCCTCCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCTTAGTGCGGGGGCAGGAtag                                  <  1:61054/150‑3 (MQ=255)
                                                                                                                    ccgcTGTATGTACTCCCTCCGTGTGAACCGCCACCGCACTCTTTTTTTCACCCTCGCCACTCCGTGGATAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                                <  2:1366423/149‑1 (MQ=255)
                                                                                                                      gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                              >  2:1813297/1‑150 (MQ=255)
                                                                                                                             taacaccccCACTCTCAAACCCCACCTCCCTTTTTTTACCCCCCGCCCCCCGCGGGGGGGTGGGCGGGGGGGCGGGCAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa                       <  2:1474806/145‑1 (MQ=255)
                                                                                                                               accccctccctCTTACCCCACCCCCCACTCTTTTTTTCCCCCCCGCCCCGCGGGGGTGTGGGGGGGGGGGTGGGCAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt                     <  2:886750/147‑1 (MQ=255)
                                                                                                                                cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACtt                    >  1:1852388/1‑150 (MQ=255)
                                                                                                                                 ccccTCACTCTGACACGACACCGCACTATTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt                   <  2:910138/149‑1 (MQ=255)
                                                                                                                                  cctacACTCTAAAACCACCCCGCATTATTTTTTTCTCCATCCCCCCTCCGGGTTGTGGGCGGGGGTTTGGGGACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa                  <  2:1665850/146‑1 (MQ=255)
                                                                                                                                  cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCATCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTATGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa                  >  2:1610751/1‑150 (MQ=255)
                                                                                                                                  cccTCACTCTGAAACGACACAGCACTCATATATTCTCCCTCGCCACTCCGGGGAGAGGTCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTATCTGCTATCCTGCCCCCGCTCTTAGGGGGCAGTATGCTAAACATCCTACTTTa                  >  1:940579/1‑150 (MQ=255)
                                                                                                                                   ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc                 >  2:337292/1‑150 (MQ=255)
                                                                                                                                   ccTCACGCTGAAACGACACCCCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc                 <  2:306382/150‑1 (MQ=255)
                                                                                                                                        ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc            >  2:651999/1‑150 (MQ=255)
                                                                                                                                        ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc            >  2:615474/1‑150 (MQ=255)
                                                                                                                                          ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcaat          >  2:1232544/1‑147 (MQ=255)
                                                                                                                                          ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg          >  1:1781330/1‑150 (MQ=255)
                                                                                                                                          ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg          >  1:465232/1‑150 (MQ=255)
                                                                                                                                              aaCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAt      >  1:644403/1‑150 (MQ=255)
                                                                                                                                               aCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTATACATCGTACTTTACGAACCAGAAAtt     >  1:1468596/1‑150 (MQ=255)
                                                                                                                                                  acacCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCg  <  2:1096415/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCG  >  NC_000913/3780955‑3781250

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