Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,688,242 C→G 8.6% intergenic (+60/‑171) recJ → / → PP_1478 ssDNA‑specific exonuclease RecJ/xenobiotic reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,688,2420CG 8.6% 881.3 / 55.7 268intergenic (+60/‑171)recJ/PP_1478ssDNA‑specific exonuclease RecJ/xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base C (128/117);  new base G (10/13);  total (138/130)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 294 total aligned reads displayed.

GGGTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCACCATTG  >  NC_002947/1688094‑1688391
                                                                                                                                                    |                                                                                                                                                     
gggTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCcgcg                                                                                                                                                      <  2:6237219/150‑1 (MQ=255)
 ggTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCcgcgc                                                                                                                                                     >  2:3068922/1‑150 (MQ=255)
 ggTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCcgcgc                                                                                                                                                     <  1:2651508/150‑1 (MQ=255)
  gTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCcgcgcg                                                                                                                                                    <  2:5775047/150‑1 (MQ=255)
  gTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCcgcgcg                                                                                                                                                    >  2:3430767/1‑150 (MQ=255)
    ggAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCgcggc                                                                                                                                                  <  2:3610073/150‑1 (MQ=255)
       gTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGGTGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                               >  2:6637414/1‑150 (MQ=255)
       gTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                               <  2:2287755/150‑1 (MQ=255)
          ggCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                            <  1:1187409/150‑1 (MQ=255)
           gccaaCAAGCAGGATGAGAACGAGTTGCTCGGTAATGAGAGTGACAACTTGATCATTGCGCATATGGATCAGACCAGAGTATAGCTTCATGAGCCTGCCATTGAACTTGCGCTGGGGCACAAATCGAGAGCCGCCCGCGCGGCGCTCGAt                                                                                                                                           >  2:1116623/5‑150 (MQ=255)
           gCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                                           >  1:2707307/1‑150 (MQ=255)
               aCAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGTCACAAAACGAGCGCCGCCCGCGCGGCGCTCGATCTca                                                                                                                                       >  1:4219381/1‑150 (MQ=255)
                cAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAGATCGAGCGCCGCGCGGGCGGCGCTCGAtttgtg                                                                                                                                      <  2:1043577/150‑6 (MQ=255)
                 aaGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTcaca                                                                                                                                     >  1:5016771/1‑150 (MQ=255)
                 aaGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTcaca                                                                                                                                     >  1:237734/1‑150 (MQ=255)
                  aGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAg                                                                                                                                    <  1:6090581/150‑1 (MQ=255)
                      ggATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGACCGCGCGGCGCTCGATCTCACAGgcgc                                                                                                                                >  2:5009938/1‑150 (MQ=255)
                         tgtgCACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGc                                                                                                                             <  2:1841508/150‑1 (MQ=255)
                          gtgAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCa                                                                                                                            >  1:4354450/1‑150 (MQ=255)
                            gAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCaaa                                                                                                                          >  1:3796014/1‑150 (MQ=255)
                             aaCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAAc                                                                                                                         >  2:1319070/1‑150 (MQ=255)
                             aaCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAAc                                                                                                                         >  2:4781558/1‑150 (MQ=255)
                             aaCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAAc                                                                                                                         >  2:199607/1‑150 (MQ=255)
                             aaCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAAc                                                                                                                         >  2:697219/1‑150 (MQ=255)
                                  gTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTc                                                                                                                    <  2:3206944/150‑1 (MQ=255)
                                  gTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTc                                                                                                                    >  2:3858689/1‑150 (MQ=255)
                                   tttCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCg                                                                                                                   >  2:4479908/1‑150 (MQ=255)
                                   tttCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCg                                                                                                                   >  1:6264814/1‑150 (MQ=255)
                                   tttCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCg                                                                                                                   >  1:6479202/1‑150 (MQ=255)
                                         ggTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAAc                                                                                                             >  2:2963345/1‑150 (MQ=255)
                                                agagTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGc                                                                                                      >  2:4052012/1‑150 (MQ=255)
                                                     gtgCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATcc                                                                                                 <  2:4336688/150‑1 (MQ=255)
                                                     gtgCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATcc                                                                                                 <  2:1684331/150‑1 (MQ=255)
                                                        cAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCcac                                                                                              >  1:2037001/1‑150 (MQ=255)
                                                         aGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCcaca                                                                                             <  2:4133977/150‑1 (MQ=255)
                                                         aGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCcaca                                                                                             >  1:2271338/1‑150 (MQ=255)
                                                           ttGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGa                                                                                           >  2:6801281/1‑150 (MQ=255)
                                                             gatgatTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAAc                                                                                         >  2:6772904/1‑150 (MQ=255)
                                                                       cATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCccac                                                                               <  2:2627877/150‑1 (MQ=255)
                                                                        atatGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCccacc                                                                              <  2:5212926/150‑1 (MQ=255)
                                                                           tGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAAt                                                                           >  2:5186306/1‑150 (MQ=255)
                                                                           tGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAAt                                                                           >  2:1724183/1‑150 (MQ=255)
                                                                            ggAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATc                                                                          <  1:6769133/150‑1 (MQ=255)
                                                                              aGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGGTGTAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTg                                                                        <  2:3783988/150‑1 (MQ=255)
                                                                               gCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTgg                                                                       <  1:4440724/150‑1 (MQ=255)
                                                                                 cgcgCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTc                                                                     <  2:2441919/150‑1 (MQ=255)
                                                                                  gcgcTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCt                                                                    >  1:3187318/1‑150 (MQ=255)
                                                                                   cgcTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTa                                                                   <  1:6596979/150‑1 (MQ=255)
                                                                                        aGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCACTCTGGTCTAGtctc                                                              >  1:1484709/1‑150 (MQ=255)
                                                                                          tttAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCAt                                                            >  1:1889119/1‑150 (MQ=255)
                                                                                          tttAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCAt                                                            >  1:3473378/1‑150 (MQ=255)
                                                                                             aGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTaa                                                         <  1:2543160/150‑1 (MQ=255)
                                                                                             aGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTaa                                                         <  1:3887091/150‑1 (MQ=255)
                                                                                                gTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGa                                                      <  1:4500088/150‑1 (MQ=255)
                                                                                                gTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGa                                                      >  2:5625737/1‑150 (MQ=255)
                                                                                                   aGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCg                                                   <  1:1013885/150‑1 (MQ=255)
                                                                                                    ggTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGTGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCgg                                                  >  2:3425170/1‑150 (MQ=255)
                                                                                                    ggTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCCATCTGGTCTAGTCTCATTAATGACCgg                                                  >  1:329949/1‑150 (MQ=255)
                                                                                                          tgcctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCgg                                            <  1:5806648/150‑1 (MQ=255)
                                                                                                           gcctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCggg                                           >  1:3868348/1‑150 (MQ=255)
                                                                                                           gcctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCggg                                           >  2:2690787/1‑150 (MQ=255)
                                                                                                             ctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGcc                                         >  1:5900165/1‑150 (MQ=255)
                                                                                                              ttGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCa                                        >  1:693092/1‑150 (MQ=255)
                                                                                                              ttGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCa                                        <  1:5625737/150‑1 (MQ=255)
                                                                                                                ggATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAgg                                      <  2:5067074/150‑1 (MQ=255)
                                                                                                                 gATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAggg                                     <  1:5152085/150‑1 (MQ=255)
                                                                                                                   tCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGAt                                   <  2:5272541/150‑1 (MQ=255)
                                                                                                                   tCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGAt                                   <  2:6264814/150‑1 (MQ=255)
                                                                                                                   tCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGAt                                   <  1:542771/150‑1 (MQ=255)
                                                                                                                     ttGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGAtgt                                 <  2:1609694/150‑1 (MQ=255)
                                                                                                                       gcgcTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGc                               <  1:6359303/150‑1 (MQ=255)
                                                                                                                        cgcTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCa                              <  2:3707757/150‑1 (MQ=255)
                                                                                                                           tGGGGCACAAATCGTGCGCCGCCCGCGCGGCGCTCGATCTCACAGGAGCTGCCATCCTTTCGGCGCACCCCTCGCCATCCCACAGAACATCCCCCCACCAATCTGGTCTAGTCTCATTAAAGACCGGCACCGGGCCCGGGATGGGCAcgt                           >  1:1232272/1‑148 (MQ=38)
                                                                                                                           tGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACtt                           <  2:6479202/150‑1 (MQ=255)
                                                                                                                             gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                         <  1:3021966/150‑1 (MQ=255)
                                                                                                                             gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                         <  2:5565912/150‑1 (MQ=255)
                                                                                                                             gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                         <  2:5078055/150‑1 (MQ=255)
                                                                                                                             gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                         <  1:1619573/150‑1 (MQ=255)
                                                                                                                             gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                         <  2:2100770/150‑1 (MQ=255)
                                                                                                                             gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                         <  2:2193288/150‑1 (MQ=255)
                                                                                                                               gCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgt                       <  2:2365762/150‑1 (MQ=255)
                                                                                                                               gCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgt                       <  1:3068922/150‑1 (MQ=255)
                                                                                                                               gCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgt                       <  2:6016075/150‑1 (MQ=255)
                                                                                                                                    gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCg                  >  1:5488793/2‑150 (MQ=255)
                                                                                                                                    aaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCg                  >  1:4854103/1‑150 (MQ=255)
                                                                                                                                      tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCggg                >  1:683892/1‑150 (MQ=255)
                                                                                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGcc              >  2:2742947/1‑150 (MQ=255)
                                                                                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGcc              >  2:798062/1‑150 (MQ=255)
                                                                                                                                         aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAATCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCg             >  2:6108968/1‑150 (MQ=255)
                                                                                                                                          gcgcCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCgg            >  2:4886641/1‑150 (MQ=255)
                                                                                                                                           cgccgcCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCCATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGa           >  2:895123/1‑150 (MQ=255)
                                                                                                                                           cgccgcCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGa           >  1:4093953/1‑150 (MQ=255)
                                                                                                                                           cgccgcCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAAGCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGa           <  2:4229939/150‑1 (MQ=255)
                                                                                                                                              cgcccgcGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCa        >  1:792515/1‑150 (MQ=255)
                                                                                                                                              cgcccgcGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCa        <  2:3868405/150‑1 (MQ=255)
                                                                                                                                                   gcgcgGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCACCAtt   <  1:2180838/150‑1 (MQ=255)
                                                                                                                                                    cgcgGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCACTCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCACCATTg  <  1:5205305/150‑1 (MQ=255)
                                                                                                                                                    cgcgGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCACCATTg  <  1:41820/150‑1 (MQ=255)
                                                                                                                                                    cgcgGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCACCATTg  <  1:3429498/150‑1 (MQ=255)
                                                                                                                                                    cgcgGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCACCATTg  <  1:5048128/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCACCATTG  >  NC_002947/1688094‑1688391

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