Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 69,705 C→T 5.1% V48M (GTG→ATG)  gmhB ← D‑glycero‑beta‑D‑manno‑heptose‑1,7‑bisphosphate 7‑phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294769,7050CT 5.1% 388.3 / 6.8 128V48M (GTG→ATG) gmhBD‑glycero‑beta‑D‑manno‑heptose‑1,7‑bisphosphate 7‑phosphatase
Reads supporting (aligned to +/- strand):  ref base C (59/61);  new base T (3/4);  total (62/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 149 total aligned reads displayed.

GCACCCTTCGTCCGGCCCATGCGGGCAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTAC  >  NC_002947/69556‑69853
                                                                                                                                                     |                                                                                                                                                     
 cACCCTTCGTCCGGCCCATGCGGGCAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTggccacg                                                                                                                                                      <  1:282027/150‑1 (MQ=255)
    ccTTCGTCCGGCCCATGCGGGCAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGTGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTggccacggcc                                                                                                                                                   >  2:1737247/1‑150 (MQ=255)
    ccTTCGTCCGGCCCATGCGGGCAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTggccacggcc                                                                                                                                                   >  2:2097123/1‑150 (MQ=255)
    ccTTCGTCCGGCCCATGCGGGCAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTggccacggcc                                                                                                                                                   >  1:1831440/1‑150 (MQ=255)
      ttCGTCCGGCCCATGCGGGCAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAAGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTggccacggccac                                                                                                                                                 >  2:2130018/1‑150 (MQ=255)
       tggaccGGACCATTTGGGCAATACACGATAAGGCCCACCACGCCTACCTGCTCGTCCACCAGCCCGCGTAGGCACGCATGTCTGGCCTTAAGTGTCTCCAGCCAATAGTAGCCACGAGCAATGCCGGACTGGTTGGTGGCCACGGCCAcc                                                                                                                                                <  2:303418/146‑1 (MQ=255)
         gTCCGGCCCATGCGGGCAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCGTGGCCAccaa                                                                                                                                              <  2:1319724/150‑3 (MQ=255)
          tCCGGCCCATGCGGGCAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGGCTGGACGGTGGCCGTGGCCAccaac                                                                                                                                             <  2:1615663/150‑4 (MQ=255)
              gCCCATGCGGGCAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGc                                                                                                                                         >  1:2036159/1‑150 (MQ=255)
                   tGCGGGCAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGcc                                                                                                                                    <  2:1765755/150‑1 (MQ=255)
                     cGGGTAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCtt                                                                                                                                  <  2:2096605/150‑1 (MQ=255)
                     cGGGCAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCtt                                                                                                                                  <  2:2036159/150‑1 (MQ=255)
                        tcAATATACGATAGGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCt                                                                                                                               <  2:624087/149‑1 (MQ=255)
                           aTACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCaa                                                                                                                            <  1:2230704/150‑1 (MQ=255)
                            tACACGATCTGGCCCACCTCACCGCCCGGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAAc                                                                                                                           <  2:785780/150‑1 (MQ=255)
                            tACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTggccacggcca                                                                                                                                                  >  1:1587190/1‑127 (MQ=255)
                             acacGATATCGCCCACCTCACCGCCATGCTCGGCGACCAGAGCGGGCAGACGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATATTAGCCACGGGCAATGCCGGAATGGTTGGTGGCCACGGCCACCGGCCAGCACGCCTTGCGCAACt                                                                                                                          >  2:322451/1‑150 (MQ=255)
                               acGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCACCTgc                                                                                                                        >  2:464148/1‑150 (MQ=255)
                                  atatGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTCCTCAACcgcgcg                                                                                                                    >  2:161350/1‑150 (MQ=255)
                                   tatGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTggccacggcca                                                                                                                                                  >  1:297662/1‑120 (MQ=255)
                                    atGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGAt                                                                                                                  <  2:26429/150‑1 (MQ=255)
                                     tGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATc                                                                                                                 <  1:666676/150‑1 (MQ=255)
                                      ggCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCg                                                                                                                <  2:1831440/150‑1 (MQ=255)
                                       gCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGc                                                                                                               >  2:532001/1‑150 (MQ=255)
                                        cccACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCa                                                                                                              <  2:1152888/150‑1 (MQ=255)
                                        cccACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCa                                                                                                              >  1:1724069/1‑150 (MQ=255)
                                            ccTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGa                                                                                                          <  1:1287030/150‑1 (MQ=255)
                                             cTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGAc                                                                                                         >  2:1247706/1‑150 (MQ=255)
                                              tCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGATTCGCATGATCGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATcagat                                                                                                        >  2:1270098/1‑146 (MQ=255)
                                                     ccTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCgg                                                                                                 <  1:555942/150‑1 (MQ=255)
                                                       tGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGc                                                                                                                                <  1:917391/118‑1 (MQ=255)
                                                       tGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGc                                                                                                                                >  2:917391/1‑118 (MQ=255)
                                                          tCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTCGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATc                                                                                            <  1:1952395/150‑1 (MQ=255)
                                                            ggCCACCAGCCCGCCCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATACTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCCCGGCCCCCCTCCAGCCCGCCCTCCCCAACTGCGCCACCGCCTCGACCGCACCGCGGaaccc                                                                                          >  2:2172053/1‑146 (MQ=255)
                                                               caccaGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTTGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCAc                                                                                       <  1:900882/150‑1 (MQ=255)
                                                                accaGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCAcc                                                                                      <  1:1153540/150‑1 (MQ=255)
                                                                 ccaGCGCTCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCAccc                                                                                     >  1:885454/1‑150 (MQ=255)
                                                                  caGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCa                                                                                    <  2:285777/150‑1 (MQ=255)
                                                                  caGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGACCGGCACCCa                                                                                    >  2:621051/1‑150 (MQ=255)
                                                                   catgcgcgcCTGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATc                                                                                            <  2:1042400/138‑1 (MQ=255)
                                                                   catgcgcgcCTGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATc                                                                                            >  1:1042400/4‑141 (MQ=255)
                                                                    gcgcacgcAGGAGCCCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCAtt                                                                                  >  2:92802/1‑150 (MQ=255)
                                                                     cgcgcgcAGGCGCGCATGCATGGCCTCCAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCAttc                                                                                 <  2:306280/150‑1 (MQ=255)
                                                                          gcAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCa                                                                            <  1:2052801/150‑1 (MQ=255)
                                                                            aGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGc                                                                          >  1:1218418/1‑150 (MQ=255)
                                                                            aGGCGCGCATGCATGACCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTgcgc                                                                                                                      <  2:1963154/107‑1 (MQ=255)
                                                                            aGGCGCGCATGCATGACCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTgcgc                                                                                                                      >  1:1963154/1‑107 (MQ=255)
                                                                                gcgcATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCt                                                                      <  1:522782/150‑1 (MQ=255)
                                                                                 cgcATGCATGGCCTCAAGGGTTGCTAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGTCACCGCCCCCGTCCAGCCCGCCTTGCTCAACTGCGCGCTCCCATCGACCGAACCGCGCCTCGGCACCCATTCTTCCAGCGTcgg                                                                     <  2:702904/150‑3 (MQ=255)
                                                                                 cgcATGCATGGCCTCAAGGGTTGCCAGCGGATAGTATCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCtt                                                                     <  1:1737247/150‑1 (MQ=255)
                                                                                 cgcATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCtt                                                                     <  1:2097123/150‑1 (MQ=255)
                                                                                       catgGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCCTTCTTCAAGCGTCTTCATTTa                                                               >  2:1714906/1‑150 (MQ=255)
                                                                                       catgGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCAt                                                                                   >  1:2165640/1‑130 (MQ=255)
                                                                                       catgGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCAt                                                                                   <  2:2165640/130‑1 (MQ=255)
                                                                                        atgGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCa                                                                                                                                           >  2:297074/1‑74 (MQ=255)
                                                                                        atgGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCa                                                                                                                                           <  1:297074/74‑1 (MQ=255)
                                                                                        atgGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTAc                                                              <  1:1392436/150‑1 (MQ=255)
                                                                                        atgGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTAc                                                              <  2:1724069/150‑1 (MQ=255)
                                                                                             cgcAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCACTCTTCCAGCGTCTTGATGTACGCATc                                                         >  1:2270143/3‑150 (MQ=255)
                                                                                                 aGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCGACCGTACAGCCCGACTTGCTAAACACGCATGCCACAGCTGCCAAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAg                                                     >  1:550187/1‑150 (MQ=255)
                                                                                                  gggTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGt                                                    <  1:1714906/150‑1 (MQ=255)
                                                                                                   ggTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTc                                                   <  2:885454/150‑1 (MQ=255)
                                                                                                      tGCCAGCGGATAGTAGCCACGGGACCCCCCTGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTa                                                <  1:943/150‑1 (MQ=255)
                                                                                                         cAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATc                                                                                            <  1:2357668/103‑1 (MQ=255)
                                                                                                         cAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATc                                                                                            >  2:2357668/1‑103 (MQ=255)
                                                                                                         cAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGtt                                             >  1:27185/1‑150 (MQ=255)
                                                                                                          aGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGGCCAGCCCGCCTTGCTCAACTGCGCGAGCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTg                                            >  2:2350161/1‑150 (MQ=255)
                                                                                                           gCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGa                                           >  2:1010900/1‑150 (MQ=255)
                                                                                                            cGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGAt                                          <  2:254784/150‑1 (MQ=255)
                                                                                                            cGGATAGTAGCCAAGGGCAATGACGGACAGGATGGTGGCCACGGCCACCTTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGAt                                          >  2:1481611/1‑150 (MQ=255)
                                                                                                                   tagCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGAt                                                                                             <  1:2402170/92‑1 (MQ=255)
                                                                                                                     gCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTc                                 >  1:2112491/1‑150 (MQ=255)
                                                                                                                      ccACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGc                                                                                         <  2:943287/93‑1 (MQ=255)
                                                                                                                      ccACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGc                                                                                         >  1:943287/1‑93 (MQ=255)
                                                                                                                      ccACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCa                                >  1:392559/1‑150 (MQ=255)
                                                                                                                       cACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCAc                               >  1:348270/1‑150 (MQ=255)
                                                                                                                        aCGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCAc                                       <  2:416273/141‑1 (MQ=255)
                                                                                                                        aCGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACg                              >  1:1441489/1‑150 (MQ=255)
                                                                                                                             cAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCga                         <  1:1197752/150‑1 (MQ=255)
                                                                                                                              aaTGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCgag                        >  2:866014/1‑150 (MQ=255)
                                                                                                                                   cGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCAtt                                                                                  <  2:1538969/87‑1 (MQ=255)
                                                                                                                                   cGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCAtt                                                                                  >  1:1538969/1‑87 (MQ=255)
                                                                                                                                     gACTGGTTGGTGGCCACGGCCAGCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGtgg                                            <  1:148725/123‑3 (MQ=255)
                                                                                                                                     gACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGtgg                                            >  2:148725/1‑121 (MQ=255)
                                                                                                                                      ccTGGTTGGTGGCCAAGGCCCCCGTCCAGCCCGCCTTGCTCAACCGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATcagca                <  1:2188515/149‑1 (MQ=255)
                                                                                                                                        tggttggtGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGt              >  2:1554040/1‑150 (MQ=255)
                                                                                                                                        tggttggtGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGt              >  2:1386478/1‑150 (MQ=255)
                                                                                                                                        tggttggtGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCCTCATGAGGGACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGt              >  2:2131775/1‑150 (MQ=255)
                                                                                                                                          gttggtGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGttt            >  2:321394/1‑150 (MQ=255)
                                                                                                                                          gacggtggCCGTGGCCACCAACCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGttt            >  2:2034529/4‑150 (MQ=255)
                                                                                                                                           ttggtGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCgg                                                                                                 <  2:969573/64‑1 (MQ=255)
                                                                                                                                           ttggtGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCAc                                                                                       <  2:544433/74‑1 (MQ=255)
                                                                                                                                           ttggtGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCAc                                                                                       >  1:544433/1‑74 (MQ=255)
                                                                                                                                            cggtggCCGTGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCTAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATcagca                <  2:770853/143‑1 (MQ=255)
                                                                                                                                            cggtggCCGTGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCTAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATcagca                >  1:770853/2‑144 (MQ=255)
                                                                                                                                              gtggCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAg        <  2:1218418/150‑1 (MQ=255)
                                                                                                                                                tgccacggccacCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAg                                               <  1:1140941/108‑1 (MQ=255)
                                                                                                                                                ggccacggccacCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAg                                               >  2:1140941/1‑109 (MQ=255)
                                                                                                                                                   cacggccacCGTCCAGCCCGCCTTGCTCAACTGCGCGATCACATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCaa         >  2:789817/1‑144 (MQ=255)
                                                                                                                                                   cacggccacCGTCCAGCCCGCCTTGCTCAACTGCGCGATCACATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCaa         <  1:789817/144‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCACCCTTCGTCCGGCCCATGCGGGCAATACACGATATGGCCCACCTCACCGCCCTGCTCGGCCACCAGCGCGCGCAGGCGCGCATGCATGGCCTCAAGGGTTGCCAGCGGATAGTAGCCACGGGCAATGCCGGACTGGTTGGTGGCCACGGCCACCGTCCAGCCCGCCTTGCTCAACTGCGCGATCGCATCGACCGAACCGGGGATCGGCACCCATTCTTCCAGCGTCTTGATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTAC  >  NC_002947/69556‑69853

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