Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,145,781 G→A 7.8% R839Q (CGG→CAG)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,145,7810GA 7.8% 510.3 / 21.6 153R839Q (CGG→CAG) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (74/67);  new base A (6/6);  total (80/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

GACTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAAT  >  NC_000913/2145633‑2145929
                                                                                                                                                    |                                                                                                                                                    
taCTATCACTTTATATGGGTAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGTCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTggccggggccgg                                                                                                                                                     <  1:390758/149‑1 (MQ=255)
 aCTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGgccggggccggg                                                                                                                                                    >  2:612381/1‑150 (MQ=255)
 aCTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGgccggggccggg                                                                                                                                                    >  1:1651661/1‑150 (MQ=255)
    aTCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCAGGTGa                                                                                                                                                 >  2:815081/1‑150 (MQ=255)
    aTCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCAGGTGa                                                                                                                                                 >  1:162508/1‑150 (MQ=255)
    aTCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCAGGTGa                                                                                                                                                 >  1:684977/1‑150 (MQ=255)
      cACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACg                                                                                                                                               >  1:2126637/1‑150 (MQ=255)
      cACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACg                                                                                                                                               >  1:1785338/1‑150 (MQ=255)
          ttATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGttt                                                                                                                                           <  2:768365/150‑1 (MQ=255)
          ttATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGttt                                                                                                                                           <  2:248154/150‑1 (MQ=255)
                gggAAGGCCGGGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGAGTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAAc                                                                                                                                     <  2:2048197/150‑1 (MQ=255)
                 ggAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAAcc                                                                                                                                    >  1:2257046/1‑150 (MQ=255)
                   aaGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCgc                                                                                                                                  >  1:1431397/1‑150 (MQ=255)
                        cGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCaa                                                                                                                             >  1:897519/1‑150 (MQ=255)
                        cGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCaa                                                                                                                             >  1:293199/1‑150 (MQ=255)
                                   ggTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGcaca                                                                                                                  <  2:670398/150‑1 (MQ=255)
                                      aGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAg                                                                                                               <  1:255913/150‑1 (MQ=255)
                                      aGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAg                                                                                                               <  1:842344/150‑1 (MQ=255)
                                      aGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAg                                                                                                               >  2:1701435/1‑150 (MQ=255)
                                         tGCCAGTGCCGGGAGTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGa                                                                                                            <  2:1550309/150‑1 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           >  1:2322039/1‑150 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           >  2:633336/1‑150 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           >  2:1294245/1‑150 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           >  1:338951/1‑150 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           >  2:1231472/1‑150 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           >  2:1231469/1‑150 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           >  2:1141285/1‑150 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           >  2:1826036/1‑150 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           >  2:2114470/1‑150 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           >  1:1686669/1‑150 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           >  2:1825951/1‑150 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           >  2:1738725/1‑150 (MQ=255)
                                          gCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGaa                                                                                                           <  2:676456/150‑1 (MQ=255)
                                               tGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGggc                                                                                                      <  1:1826036/150‑1 (MQ=255)
                                                gCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGggcg                                                                                                     <  1:1079835/150‑1 (MQ=255)
                                                gCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGggcg                                                                                                     <  1:729692/150‑1 (MQ=255)
                                                   gggATTACCTTTATTGATGACAACAATCATCAGGATGCATATGTGATGTCGCAGGCTGTTATCGACTGTTATGCCTCCAAATATGGTGGCCGGGGCCGGGTGTCGGTTTACGATCAGCAGCAAGCTGCCGCACATAGCGTACGGgcggcg                                                                                                  <  2:2464624/150‑1 (MQ=255)
                                                   gggATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGgcggcg                                                                                                  <  2:1432094/150‑1 (MQ=255)
                                                   gggATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGgcggcg                                                                                                  <  2:1432377/150‑1 (MQ=255)
                                                          ccTTGATTGATGACAACACTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGc                                                                                           <  2:2608319/150‑1 (MQ=255)
                                                          ccTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGc                                                                                           <  2:1868/150‑1 (MQ=255)
                                                          ccTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGc                                                                                           <  2:2289162/150‑1 (MQ=255)
                                                          ccTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGc                                                                                           <  2:915580/150‑1 (MQ=255)
                                                          ccTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGc                                                                                           <  2:612526/150‑1 (MQ=255)
                                                          ccTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGc                                                                                           >  1:924986/1‑150 (MQ=255)
                                                          ccTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGc                                                                                           >  1:1042978/1‑150 (MQ=255)
                                                               attgatGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTtgat                                                                                      <  1:438243/150‑1 (MQ=255)
                                                                  gatgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGaa                                                                                   <  1:269749/150‑1 (MQ=255)
                                                                   atgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAAc                                                                                  >  1:1007394/1‑150 (MQ=255)
                                                                   atgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAAc                                                                                  >  1:1455811/1‑150 (MQ=255)
                                                                   atgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAAc                                                                                  >  1:1866659/1‑150 (MQ=255)
                                                                    tgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACa                                                                                 <  1:2063637/150‑1 (MQ=255)
                                                                     gaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAg                                                                                >  1:1207322/1‑150 (MQ=255)
                                                                      acaacaATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGt                                                                               <  1:1628017/150‑1 (MQ=255)
                                                                        aacaaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTgg                                                                             >  2:1107510/1‑150 (MQ=255)
                                                                        aacaaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTgg                                                                             >  1:1916686/1‑150 (MQ=255)
                                                                         acacTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCAGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGc                                                                            <  1:2481351/150‑1 (MQ=255)
                                                                         acaaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGc                                                                            <  1:567652/150‑1 (MQ=255)
                                                                         acaaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGc                                                                            <  2:227771/150‑1 (MQ=255)
                                                                         acaaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGc                                                                            <  2:1954800/150‑1 (MQ=255)
                                                                         acaaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCAGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGc                                                                            <  1:656159/150‑1 (MQ=255)
                                                                         acaaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCAGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGc                                                                            <  2:684914/150‑1 (MQ=255)
                                                                         acaaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCAGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGc                                                                            <  2:162508/150‑1 (MQ=255)
                                                                          cactcatcaGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTCGGCCGCCACAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACAGAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCg                                                                           <  1:961412/147‑1 (MQ=255)
                                                                          caaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAACTGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCg                                                                           <  2:1538813/150‑1 (MQ=255)
                                                                          caaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCg                                                                           <  1:2358208/150‑1 (MQ=255)
                                                                          caaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCg                                                                           <  2:36222/150‑1 (MQ=255)
                                                                          caaTCATCAGGCGGCAGAAGTGATGTCACAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCg                                                                           <  1:2395998/150‑1 (MQ=255)
                                                                                    gcggcAGAAGTGATGTCGCAGACTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTaaa                                                                 >  2:1639188/1‑150 (MQ=255)
                                                                                     cggcAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAag                                                                <  2:459405/150‑1 (MQ=255)
                                                                                     cggcAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAag                                                                <  2:1204618/150‑1 (MQ=255)
                                                                                      ggcAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAaga                                                               >  2:1205868/1‑150 (MQ=255)
                                                                                          gAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAAt                                                           <  2:525447/150‑1 (MQ=255)
                                                                                          gAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAAt                                                           >  2:978379/1‑150 (MQ=255)
                                                                                          gAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAAt                                                           <  1:2454386/150‑1 (MQ=255)
                                                                                          gAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATCATTAAAGAGAAt                                                           >  2:2120168/1‑150 (MQ=255)
                                                                                           aaGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATc                                                          >  2:645249/1‑150 (MQ=255)
                                                                                            aGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCa                                                         >  2:1098207/1‑150 (MQ=255)
                                                                                            aGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCa                                                         >  1:501181/1‑150 (MQ=255)
                                                                                              tGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGt                                                       <  2:2257046/150‑1 (MQ=255)
                                                                                                 tGTCGCAGGCTGATATCGCCTATTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTtga                                                    <  2:712863/150‑1 (MQ=255)
                                                                                                  gTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTtgat                                                   >  1:2633179/1‑150 (MQ=255)
                                                                                                   tCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTtgatg                                                  <  2:877462/150‑1 (MQ=255)
                                                                                                   tCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTtgatg                                                  >  1:1951765/1‑150 (MQ=255)
                                                                                                   tCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTtgatg                                                  <  2:235417/150‑1 (MQ=255)
                                                                                                   tCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTtgatg                                                  >  1:2232590/1‑150 (MQ=255)
                                                                                                        ggCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCt                                             >  2:2626579/1‑150 (MQ=255)
                                                                                                           tGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGc                                          >  1:1196306/1‑150 (MQ=255)
                                                                                                           tGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGc                                          >  1:1196769/1‑150 (MQ=255)
                                                                                                            gATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGcc                                         >  1:1231761/1‑150 (MQ=255)
                                                                                                                         aTGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCtt                            >  1:2611477/1‑150 (MQ=255)
                                                                                                                              tCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCa                       <  2:337974/150‑1 (MQ=255)
                                                                                                                              tCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCa                       >  2:2483848/1‑150 (MQ=255)
                                                                                                                              tCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCa                       <  1:1005866/150‑1 (MQ=255)
                                                                                                                              tCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCa                       >  2:1933749/1‑150 (MQ=255)
                                                                                                                                   aaaTGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGAt                  >  1:2621293/1‑150 (MQ=255)
                                                                                                                                   aaaTGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGAt                  >  1:1589009/1‑150 (MQ=255)
                                                                                                                                   aaaTGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGAt                  >  2:2567701/1‑150 (MQ=255)
                                                                                                                                          ggccggggccggGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGaa           <  2:2181880/150‑1 (MQ=255)
                                                                                                                                                   cgggTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAAt  <  2:897519/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAAT  >  NC_000913/2145633‑2145929

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