Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,459,121 A→G 6.3% S316G (AGC→GGC)  sthA → pyridine nucleotide transhydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,459,1210AG 6.3% 284.2 / 6.2 95S316G (AGC→GGC) sthApyridine nucleotide transhydrogenase
Reads supporting (aligned to +/- strand):  ref base A (45/44);  new base G (4/2);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

GTGCAACGGCCGTACCGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAGGAGCTGAC  >  NC_002947/2458973‑2459269
                                                                                                                                                    |                                                                                                                                                    
gtgCAACGGCCGTACCGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAg                                                                                                                                                     >  1:1592428/1‑150 (MQ=255)
 tgCAACGGCCGTACCGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGc                                                                                                                                                    <  1:2002599/150‑1 (MQ=255)
 tgCAACGGCCGTACCGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGc                                                                                                                                                    <  1:2382342/150‑1 (MQ=255)
 tgCAACGGCCGTACCGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGc                                                                                                                                                    <  1:1480504/150‑1 (MQ=255)
  gCAACGGCCGTACCGGCAACACTGAAAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGcc                                                                                                                                                   >  2:2121572/1‑150 (MQ=255)
    aaCGGCCGTACCGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTg                                                                                                                                                 <  1:595799/150‑1 (MQ=255)
    aaCGGCCGTACCGGCAACACTGACAAGCTGGGCCGGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTg                                                                                                                                                 <  1:621002/150‑1 (MQ=255)
      cGGCCGTACCGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACAGTACGACCGTGCCGAACAACAAGGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGc                                                                                                                                               >  2:572641/1‑150 (MQ=255)
         ccGTACCGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAg                                                                                                                                            >  2:812554/1‑150 (MQ=255)
           gTACCGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGtg                                                                                                                                          >  1:278499/1‑150 (MQ=255)
           gTACCGGCAACACTGACAAGCGGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAAATCGAGGTCGACGAGGCCTACCGTACCAGCGTGGCGAACATCTACGGTGCACGCGACGTGATCGGCTTGCCAAGCCTCGCCatcg                                                                                                                                          >  2:72244/1‑147 (MQ=255)
             ctcGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGANCNTCTACGNTGCCGGCGACGTNATCGGCTGGCCAGGCTTGGCCAgcnga                                                                                                                                        <  1:1515776/148‑5 (MQ=255)
               cGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGtgctgc                                                                                                                                      >  2:1764089/1‑150 (MQ=255)
                     cacTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGa                                                                                                                                >  1:1464989/1‑150 (MQ=255)
                     cacTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGa                                                                                                                                <  1:1144929/150‑1 (MQ=255)
                        tGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCa                                                                                                                             >  2:209925/1‑150 (MQ=255)
                        tGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCa                                                                                                                             >  2:867816/1‑150 (MQ=255)
                        tGAAAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCa                                                                                                                             <  1:2121572/150‑1 (MQ=255)
                           cAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAggg                                                                                                                          <  1:458921/150‑1 (MQ=255)
                           cAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAggg                                                                                                                          <  1:812554/150‑1 (MQ=255)
                               cTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCgg                                                                                                                      <  2:1464989/150‑1 (MQ=255)
                                    ccTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGc                                                                                                                 >  1:733830/1‑150 (MQ=255)
                                    ccTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGc                                                                                                                 >  2:168494/1‑150 (MQ=255)
                                     cTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGcc                                                                                                                >  2:779826/1‑150 (MQ=255)
                                       ggAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGc                                                                                                              <  2:733830/150‑1 (MQ=255)
                                        gAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCCTGCCGAACTCCTACGGTGCCTTCGACTTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCt                                                                                                             <  1:2098409/150‑1 (MQ=255)
                                         aaaaCATCGGCATCAAGGTGAACAGCCATGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGAc                                                                                                                               <  1:629657/131‑1 (MQ=255)
                                         aaaaCATCGGCATCAAGGTGAACAGCCATGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGAc                                                                                                                               >  2:629657/1‑131 (MQ=255)
                                              aTCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGgcagc                                                                                                       >  1:1952965/1‑150 (MQ=255)
                                               tCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGgcagca                                                                                                      <  1:168494/150‑1 (MQ=255)
                                               tCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGgcagca                                                                                                      >  1:1486295/1‑150 (MQ=255)
                                                 ggCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCAtt                                                                                                    <  1:209925/150‑1 (MQ=255)
                                                   cATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGt                                                                                                  <  2:1486295/150‑1 (MQ=255)
                                                    aTCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTc                                                                                                 <  1:779826/150‑1 (MQ=255)
                                                    aTCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTc                                                                                                 <  1:867816/150‑1 (MQ=255)
                                                      cAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGa                                                                                               <  1:72244/150‑1 (MQ=255)
                                                       aaGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGAt                                                                                              >  2:1486789/1‑150 (MQ=255)
                                                            gAACAGCCGTGGGCAGATCGAGGGCGACGAGGCCTGCCGTACCACCGTGCCGAACATCTACGGTGCCGGGGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACgg                                                                                         <  2:403890/150‑1 (MQ=255)
                                                                 gCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCt                                                                                    >  1:859314/1‑150 (MQ=255)
                                                                  ccGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTgg                                                                                                                                                >  1:1110139/1‑89 (MQ=255)
                                                                  ccGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTgg                                                                                                                                                <  2:1110139/89‑1 (MQ=255)
                                                                  ccGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGc                                                                                                                                               <  1:2344515/90‑1 (MQ=255)
                                                                  ccGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGc                                                                                                                                               >  2:2344515/1‑90 (MQ=255)
                                                                  ccGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTACCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTg                                                                                   >  2:887942/1‑150 (MQ=255)
                                                                     tGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGgcg                                                                                >  1:2365300/1‑150 (MQ=255)
                                                                         cAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCTCTGGCAGCATTGTCGATAACGGCATCTGGCGCTTc                                                                            >  2:540198/1‑150 (MQ=255)
                                                                         cAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCt                                                                                                             >  1:1299693/1‑117 (MQ=255)
                                                                         cAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCt                                                                                                             <  2:1299693/117‑1 (MQ=255)
                                                                         cAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTc                                                                            <  2:2136405/150‑1 (MQ=255)
                                                                            aTCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTc                                                                         >  2:118096/1‑150 (MQ=255)
                                                                                            ccTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACgg                                                         <  1:1486789/150‑1 (MQ=255)
                                                                                             cTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCCGGCGCTTCGTCAACGACGTGCCGACggg                                                        >  1:1437100/1‑150 (MQ=255)
                                                                                              tACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCCATCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTTGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGt                                                       >  2:660772/1‑150 (MQ=255)
                                                                                              tACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGgcag                                                                                                        >  2:392000/1‑101 (MQ=255)
                                                                                              tACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGgcag                                                                                                        <  1:392000/101‑1 (MQ=255)
                                                                                              tACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAg                                                                                      >  1:1813794/1‑119 (MQ=255)
                                                                                              tACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAg                                                                                      <  2:1813794/119‑1 (MQ=255)
                                                                                               aCCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTa                                                      >  2:1522021/1‑150 (MQ=255)
                                                                                                ccGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTAt                                                     <  1:1799042/150‑1 (MQ=255)
                                                                                                   tACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTa                                                  <  1:1522021/150‑1 (MQ=255)
                                                                                                     ccaccGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTaca                                                >  2:1255668/1‑150 (MQ=255)
                                                                                                      caccGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGTTTCGTCAACAACGCGCCGACGGGTATCTacac                                               <  2:344194/150‑1 (MQ=255)
                                                                                                      caccGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTacac                                               <  2:278499/150‑1 (MQ=255)
                                                                                                      caccGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTacac                                               <  1:1255668/150‑1 (MQ=255)
                                                                                                          gTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCAtt                                           >  1:51082/1‑150 (MQ=255)
                                                                                                          gTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCAtt                                           <  2:1592428/150‑1 (MQ=255)
                                                                                                            gCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGCGTATCTGTACAATTcc                                         >  2:516430/1‑150 (MQ=255)
                                                                                                            gCCGAACATCTACAGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCATCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTcc                                         >  2:1955498/1‑150 (MQ=255)
                                                                                                              cGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCg                                       <  2:859314/150‑1 (MQ=255)
                                                                                                               gAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCga                                      >  1:1343551/1‑150 (MQ=255)
                                                                                                               gAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCga                                      <  2:1437100/150‑1 (MQ=255)
                                                                                                                 aCATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCgaga                                    <  1:1123059/150‑1 (MQ=255)
                                                                                                                  cATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGAt                                   >  1:516131/1‑150 (MQ=255)
                                                                                                                   caatACGATGCCGGCGCGGTGATCTGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATc                                  >  2:2331917/4‑150 (MQ=255)
                                                                                                                     cTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAg                                >  2:2315511/1‑150 (MQ=255)
                                                                                                                             ccGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCg                        >  2:230862/1‑150 (MQ=255)
                                                                                                                                gCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATCCCGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCa                     >  1:1852361/1‑150 (MQ=255)
                                                                                                                                 cGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCaa                    <  2:2171659/150‑1 (MQ=255)
                                                                                                                                    cGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGaa                 <  2:1952965/150‑1 (MQ=255)
                                                                                                                                      tGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAgg                                                                                                                           >  1:1451046/1‑42 (MQ=255)
                                                                                                                                      tGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAgg                                                                                                                           <  2:1451046/42‑1 (MQ=255)
                                                                                                                                           ggCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACgg                                                         <  2:268676/103‑1 (MQ=255)
                                                                                                                                           ggCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACgg                                                         >  1:268676/1‑103 (MQ=255)
                                                                                                                                           ggCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAg          <  1:643654/150‑1 (MQ=255)
                                                                                                                                           cacTGGCCAGGCTTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGACCGGCAAGAACCAGCac          >  2:425404/3‑149 (MQ=255)
                                                                                                                                            acTGGCCAGGCTTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGa                  <  2:1668572/140‑1 (MQ=255)
                                                                                                                                            acTGGCCAGGCTTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGa                  >  1:1668572/2‑141 (MQ=255)
                                                                                                                                            acTGGCCAGGCTTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGGATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAgg         >  2:272247/2‑150 (MQ=255)
                                                                                                                                             gtGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAGGa        <  2:2254284/149‑1 (MQ=255)
                                                                                                                                             cTGGCCAGGCTTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAGGa        >  2:199738/1‑150 (MQ=255)
                                                                                                                                             cTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAGGa        <  1:1764089/150‑1 (MQ=255)
                                                                                                                                              tGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAGGAg       <  1:887942/150‑1 (MQ=255)
                                                                                                                                                gCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAATAGGGAGCAGGAgtt     >  2:1267993/1‑148 (MQ=255)
                                                                                                                                                  cAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAGGAGCTGa   >  1:1821172/1‑150 (MQ=255)
                                                                                                                                                  cAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACAACATTCCTAATGTCAGCTCGATCGGCAAGAACGAGCAGGAGCTGa   >  2:1770391/1‑150 (MQ=255)
                                                                                                                                                   aaGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTaaat                                                                      >  1:654341/1‑81 (MQ=255)
                                                                                                                                                   aaGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAac                                                                      <  2:654341/82‑1 (MQ=255)
                                                                                                                                                   aaGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAGGAg       >  2:646956/1‑145 (MQ=255)
                                                                                                                                                   aaGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAGGAGCTGAc  >  2:220461/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGCAACGGCCGTACCGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAGGAGCTGAC  >  NC_002947/2458973‑2459269

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: