Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,045,272 G→A 13.8% intergenic (+42/+64) 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,2720GA13.8% 186.3 / 20.6 80intergenic (+42/+64)PP_0905/PP_0906hypothetical protein/RND family multidrug transporter
Reads supporting (aligned to +/- strand):  ref base G (35/34);  new base A (5/6);  total (40/40)
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.71e-01

GCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAACAACGTG  >  NC_002947/1045123‑1045420
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                                                 >  1:846492/1‑89 (MQ=37)
                                                                  ttCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAAAGGCCCCGCAATCGCGAGGCCtc                                                                                                                                                 <  2:846492/90‑2 (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)
                                                                                                                                                    ggggCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAACAACGTg  <  2:587327/150‑1 (MQ=255)
                                                                                                                                                    ggggCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAACAACGTg  >  2:2413085/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAACAACGTG  >  NC_002947/1045123‑1045420

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: