Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,045,268 T→C 12.6% intergenic (+38/+68) PP_0905 → / ← PP_0906 hypothetical protein/RND family multidrug transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,045,2680TC12.6% 218.9 / 19.8 87intergenic (+38/+68)PP_0905/PP_0906hypothetical protein/RND family multidrug transporter
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base C (5/6);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

CGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAA  >  NC_002947/1045122‑1045413
                                                                                                                                                  |                                                                                                                                                  
cgcgTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCg                                                                                                                                                 <  1:1960257/150‑1 (MQ=255)
 gcgTTGCCGCGGGATCCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGAGTCCAGGGTTGCGGTCGCCGAGGTGAACCGTGCGGTAGCTGTAATAGTCAGAGGCCTCGCTATTGCgg                                                                                                                                                <  2:1843668/150‑1 (MQ=255)
   gTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGTGATGAACGCTGGGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCgggg                                                                                                                                              >  2:2099296/1‑150 (MQ=255)
   gTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCgggg                                                                                                                                              >  1:587327/1‑150 (MQ=255)
   gTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCgggg                                                                                                                                              <  1:892883/150‑1 (MQ=255)
         gcgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCCCGCAATCGCGAGGCCTCTg                                                                                                                                        <  1:402245/150‑1 (MQ=255)
         gcgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGGCACAGGGATGAAGCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTg                                                                                                                                        <  2:269486/150‑1 (MQ=255)
           gggATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGt                                                                                                                                      >  2:1871779/1‑150 (MQ=255)
            ggATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtt                                                                                                                                     <  1:1528808/150‑1 (MQ=255)
             gATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGttt                                                                                                                                    <  1:315877/150‑1 (MQ=255)
              aTGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtttt                                                                                                                                   <  1:2216316/150‑1 (MQ=255)
                gCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtttttt                                                                                                                                 <  2:1711722/150‑1 (MQ=255)
                  aaaaCCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGc                                                                                                                               <  1:367256/150‑1 (MQ=255)
                  aaaaCCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTGCCGGTTTGGCCTGCGGCAAGGGTTGCGGTCACAGGGATCAACCGTGCGTTAGCTGGAAAGGCCAGAGGCCTCGCGCTTGCCGGGCCTTTGGTTTTTTGc                                                                                                                               >  2:966958/1‑150 (MQ=255)
                        gcagcagcaGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGc                                                                                                                         >  2:250377/1‑150 (MQ=255)
                              gcagGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGcc                                                                                                                   >  1:818354/1‑150 (MQ=255)
                                agGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGGTGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCgg                                                                                                                 >  1:1834917/1‑150 (MQ=255)
                                  gCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAAAGGCCCCGCAATCGCGAGGCCtctg                                                                                                                                        >  2:2093627/1‑122 (MQ=38)
                                  gCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAAAGGCCCCGCAATCGCGAGGCCTCTgg                                                                                                                                       >  1:35457/1‑126 (MQ=38)
                                  gCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAAAGGCCCCGCAATCGCGAGGCCTCTgg                                                                                                                                       <  2:35457/126‑1 (MQ=38)
                                  gCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAAAGGCCCCGCAATCGCGAGGCCTCTg                                                                                                                                        <  1:2093627/125‑1 (MQ=255)
                                         gcgtgcCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTcgc                                                                                                       <  1:1603226/150‑1 (MQ=255)
                                          cgagcCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTcgcg                                                                                                      <  2:800525/150‑1 (MQ=255)
                                                    cAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGTGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccgc                                                                                            >  1:64488/1‑150 (MQ=255)
                                                     aGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccgcc                                                                                           <  2:1981903/150‑1 (MQ=255)
                                                      ggACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccgccc                                                                                          <  2:818354/150‑1 (MQ=255)
                                                              tGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCg                                                                                  <  1:2099296/150‑1 (MQ=255)
                                                              tGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCg                                                                                  <  2:889332/150‑1 (MQ=255)
                                                                   ttCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCa                                                                             >  2:1964548/1‑150 (MQ=255)
                                                                   ttCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAAAGGCCCCGCAATCGCGAGGCCtc                                                                                                                                          <  2:846492/90‑2 (MQ=37)
                                                                   ttCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAAAGGCCCCGCAATCGCGAGGCCtc                                                                                                                                          >  1:846492/1‑89 (MQ=37)
                                                                    tCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGCGCAGCCCGCCCCAAGGTCGGTGCaa                                                                            >  2:994185/1‑150 (MQ=255)
                                                                     ccTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTg                                                                                                                                        <  1:787472/90‑1 (MQ=255)
                                                                     ccTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTg                                                                                                                                        >  2:787472/1‑90 (MQ=255)
                                                                         tttGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGAt                                                                       >  2:192906/1‑150 (MQ=255)
                                                                          ttGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGAtt                                                                      <  1:2054823/150‑1 (MQ=255)
                                                                             gCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCa                                                                                                  >  1:1940162/1‑119 (MQ=255)
                                                                             gCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCa                                                                                                  <  2:1940162/119‑1 (MQ=255)
                                                                              ccTGCGTCAAGGGTTGTGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCATAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAg                                                                  <  2:64488/150‑1 (MQ=255)
                                                                              ccTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTg                                                                     >  2:247195/1‑147 (MQ=255)
                                                                              ccTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTg                                                                     <  1:247195/147‑1 (MQ=255)
                                                                                  cGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGc                                                              <  2:1834917/150‑1 (MQ=255)
                                                                                   gTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAg                                                                                      <  2:2304358/125‑1 (MQ=255)
                                                                                   gTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAg                                                                                      >  1:2304358/1‑125 (MQ=255)
                                                                                      aaGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCggg                                                                                                                                               <  1:980679/66‑1 (MQ=255)
                                                                                      aaGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCggg                                                                                                                                               >  2:980679/1‑66 (MQ=255)
                                                                                      aaGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTg                                                          >  1:142439/1‑150 (MQ=255)
                                                                                          gTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGgcgc                                                      >  2:425103/1‑150 (MQ=255)
                                                                                                   acaGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGCCTTATCCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTtgt                                             >  2:575710/1‑150 (MQ=255)
                                                                                                     aGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGGCGGCCTGCTCGCGGGCAGCCCGCGCCAATGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTgc                                           >  1:88799/1‑150 (MQ=255)
                                                                                                     aGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTgc                                           >  2:396414/1‑150 (MQ=255)
                                                                                                      gggATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTgcg                                          >  1:2421974/1‑150 (MQ=255)
                                                                                                         aTGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCca                                       >  2:1146610/1‑150 (MQ=255)
                                                                                                            aaCCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCcagca                                    >  2:764113/1‑150 (MQ=255)
                                                                                                            aaCCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCcagca                                    <  1:1140418/150‑1 (MQ=255)
                                                                                                            aaCCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCcagca                                    >  1:332000/1‑150 (MQ=255)
                                                                                                              cccTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAgg                                  <  1:966958/150‑1 (MQ=255)
                                                                                                              cccTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAgg                                  >  1:631444/1‑150 (MQ=255)
                                                                                                              cccTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAgg                                  <  1:1146610/150‑1 (MQ=255)
                                                                                                              cccTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCAATTGCGGGGCCTTTGGTTTTTTTCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGACTGCAGGTGCCTTGGCGCTCGCCTTGTCCGCcagcaga                                  >  2:716022/1‑149 (MQ=255)
                                                                                                               ccTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAggg                                 <  1:192906/150‑1 (MQ=255)
                                                                                                                cTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGATCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGgt                                <  1:1871779/150‑1 (MQ=255)
                                                                                                                     ttAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCTTTGTGCGCCAGCAGGGTGTAGa                           <  1:1882521/150‑1 (MQ=255)
                                                                                                                      tAGCTGGAAAAGCCAGAGGCCTCGCGATTGCTGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGCCCTCGCTTGCGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGAt                          <  1:1860753/150‑1 (MQ=255)
                                                                                                                        gCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGt                                                                                    >  1:1093297/1‑90 (MQ=255)
                                                                                                                        gCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGt                                                                                    <  2:1093297/90‑1 (MQ=255)
                                                                                                                          tGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCAc                      >  2:2149232/1‑150 (MQ=255)
                                                                                                                          tGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCAc                      >  1:2052606/1‑150 (MQ=255)
                                                                                                                             aaaaGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGc                   <  2:631444/150‑1 (MQ=255)
                                                                                                                                aGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAac                <  1:425103/150‑1 (MQ=255)
                                                                                                                                aGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAac                >  2:1447329/1‑150 (MQ=255)
                                                                                                                                  ccaaagGCCCCGCAATCGCGAGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAgg                                  <  2:1374161/126‑1 (MQ=255)
                                                                                                                                  ccAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCgg                                                                                                     >  1:1836003/1‑63 (MQ=255)
                                                                                                                                  ccAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCgg                                                                                                     <  2:1836003/63‑1 (MQ=255)
                                                                                                                                  ccAAAGGCCCCGCAATCGCGAGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAgg                                  >  1:1374161/1‑130 (MQ=255)
                                                                                                                                   caaagGCCCCGCAATCGCGAGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGc                        <  2:1209495/136‑1 (MQ=255)
                                                                                                                                   caaagGCCCCGCAATCGCGAGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGc                        >  1:1209495/4‑139 (MQ=255)
                                                                                                                                    agagGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaa            >  1:60754/1‑150 (MQ=255)
                                                                                                                                     gagGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaaa           >  1:185186/1‑150 (MQ=255)
                                                                                                                                     gagGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaaa           >  1:807359/1‑150 (MQ=255)
                                                                                                                                        gCCTCGCGATTGCGGGGCCTTTGGTTTTTTGTGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTg                                                          <  2:1557432/100‑1 (MQ=255)
                                                                                                                                        gCCTCGCGATTGCGGGGCCTTTGGTTTTTTGTGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTg                                                          >  1:1557432/1‑100 (MQ=255)
                                                                                                                                        gCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaaacaa        <  2:2396560/150‑1 (MQ=255)
                                                                                                                                         ccTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAg       <  1:250377/150‑1 (MQ=255)
                                                                                                                                          cTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAgg      >  1:1497657/1‑150 (MQ=255)
                                                                                                                                          cTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAgg      >  1:1777709/1‑150 (MQ=255)
                                                                                                                                              cgATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTaaa  <  1:396414/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAA  >  NC_002947/1045122‑1045413

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