Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,459,124 C→T 6.4% L317L (CTG→TTG)  sthA → pyridine nucleotide transhydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,459,1240CT 6.4% 276.4 / 6.8 94L317L (CTG→TTG) sthApyridine nucleotide transhydrogenase
Reads supporting (aligned to +/- strand):  ref base C (46/42);  new base T (4/2);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

GCAACGGCCGTACCGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAGGAGCTGACCCAG  >  NC_002947/2458975‑2459273
                                                                                                                                                     |                                                                                                                                                     
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:1144929/150‑1 (MQ=255)
                   cacTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGa                                                                                                                                    >  1:1464989/1‑150 (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:812554/150‑1 (MQ=255)
                         cAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAggg                                                                                                                              <  1:458921/150‑1 (MQ=255)
                             cTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCgg                                                                                                                          <  2:1464989/150‑1 (MQ=255)
                                  ccTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGc                                                                                                                     >  2:168494/1‑150 (MQ=255)
                                  ccTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGc                                                                                                                     >  1:733830/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                                                                                                                                   >  2:629657/1‑131 (MQ=255)
                                       aaaaCATCGGCATCAAGGTGAACAGCCATGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGAc                                                                                                                                   <  1:629657/131‑1 (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                                                                                                                                                    <  2:1110139/89‑1 (MQ=255)
                                                                ccGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTgg                                                                                                                                                    >  1:1110139/1‑89 (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                                                                                                                 <  2:1299693/117‑1 (MQ=255)
                                                                       cAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCt                                                                                                                 >  1:1299693/1‑117 (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                                               <  2:1592428/150‑1 (MQ=255)
                                                                                                        gTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCAtt                                               >  1:51082/1‑150 (MQ=255)
                                                                                                          gCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGCGTATCTGTACAATTcc                                             >  2:516430/1‑150 (MQ=255)
                                                                                                          gCCGAACATCTACAGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCATCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTcc                                             >  2:1955498/1‑150 (MQ=255)
                                                                                                            cGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCg                                           <  2:859314/150‑1 (MQ=255)
                                                                                                             gAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCga                                          <  2:1437100/150‑1 (MQ=255)
                                                                                                             gAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCga                                          >  1:1343551/1‑150 (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                                                             >  1:268676/1‑103 (MQ=255)
                                                                                                                                         ggCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACgg                                                             <  2:268676/103‑1 (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)
                                                                                                                                                 aatccTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAGGAg           <  1:646956/142‑1 (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)
                                                                                                                                                     cTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAac                                                                          <  2:917038/78‑1 (MQ=255)
                                                                                                                                                     cTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAac                                                                          >  1:917038/1‑78 (MQ=255)
                                                                                                                                                     cTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCACAGGCGTGCCGACGGGTATCTACACCATTCCCGCGATCAGCTCGATCGGCAAGAACGAGCAGGAGCTGACCCAg  >  1:1621738/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAACGGCCGTACCGGCAACACTGACAAGCTGGGCCTGGAAAACATCGGCATCAAGGTGAACAGCCGTGGCCAGATCGAGGTCGACGAGGCCTACCGTACCACCGTGCCGAACATCTACGGTGCCGGCGACGTGATCGGCTGGCCAAGCCTGGCCAGTGCTGCCCATGACCAGGGCCGGTCGGCCGCTGGCAGCATTGTCGATAACGGCAGCTGGCGCTTCGTCAACGACGTGCCGACGGGTATCTACACCATTCCCGAGATCAGCTCGATCGGCAAGAACGAGCAGGAGCTGACCCAG  >  NC_002947/2458975‑2459273

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: