Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,683,290 C→G 6.9% intergenic (+57/‑174) 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,683,2900CG 6.9% 673.4 / 29.8 203intergenic (+57/‑174)recJ/PP_1478ssDNA‑specific exonuclease RecJ/xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base C (86/103);  new base G (4/10);  total (90/113)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 207 total aligned reads displayed.

TGGGTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCAC  >  NC_002947/1683144‑1683437
                                                                                                                                                  |                                                                                                                                                   
 gggTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCcgcg                                                                                                                                                 <  2:4466566/150‑1 (MQ=255)
 gggTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCcgcg                                                                                                                                                 <  1:3025361/150‑1 (MQ=255)
   gTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCcgcgcg                                                                                                                                               <  2:4226071/150‑1 (MQ=255)
    tGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCgcgg                                                                                                                                              <  2:3929420/150‑1 (MQ=255)
    tGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCgcgg                                                                                                                                              <  1:7339182/150‑1 (MQ=255)
      gAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCgcggcg                                                                                                                                            <  2:3022401/150‑1 (MQ=255)
        gTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                          >  1:1014218/1‑150 (MQ=255)
         ttGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                         <  1:654967/150‑1 (MQ=255)
            gCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                      <  1:5375854/150‑1 (MQ=255)
               tACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGAtctc                                                                                                                                   >  2:5955963/1‑150 (MQ=255)
                aCAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAGATCGAGCGCCGCGCGGGCGGCGCTCGAtttgt                                                                                                                                  <  1:6716866/150‑5 (MQ=255)
                aCAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTca                                                                                                                                  <  2:8123468/150‑1 (MQ=255)
                   aGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAg                                                                                                                               >  2:7172369/1‑150 (MQ=255)
                   aGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAg                                                                                                                               >  2:4990768/1‑150 (MQ=255)
                    gCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAgg                                                                                                                              >  1:3226851/1‑150 (MQ=255)
                       ggATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGgcgc                                                                                                                           <  2:8158386/150‑1 (MQ=255)
                         aTGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTg                                                                                                                         <  1:5103255/150‑1 (MQ=255)
                          tgtgAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGc                                                                                                                        <  2:4870524/150‑1 (MQ=255)
                           gtgAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCCTATGGAGCCGCGCTGTGTTTAGCGTCCGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGTGCGCCGCCCGCGCGGCGCTCGCTCTCACAGGCGCGGCa                                                                                                                       >  2:5232231/1‑150 (MQ=255)
                             gAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCaaa                                                                                                                     <  1:5842142/150‑1 (MQ=255)
                             gAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCaaa                                                                                                                     <  1:8489727/150‑1 (MQ=255)
                              aaCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGTTCGATCTCACAGGCGCTGCAAAc                                                                                                                    >  1:261021/1‑150 (MQ=255)
                              aaCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAAc                                                                                                                    <  2:6052297/150‑1 (MQ=255)
                               aCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAAcc                                                                                                                   >  1:8917652/1‑150 (MQ=255)
                                cGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCt                                                                                                                  <  1:1094257/150‑1 (MQ=255)
                                    tttCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCg                                                                                                              >  1:7008482/1‑150 (MQ=255)
                                     ttCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCgg                                                                                                             <  2:2624245/150‑1 (MQ=255)
                                                gagaGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCg                                                                                                  >  1:8959022/1‑150 (MQ=255)
                                                gagaGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCg                                                                                                  >  2:4133026/1‑150 (MQ=255)
                                                gagaGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCg                                                                                                  <  1:3311050/150‑1 (MQ=255)
                                                gagaGTGTGAAGTTGATGATAGCGCAAATGGAGAAGCGCTGAGTTTAGCGTCAGGAGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCg                                                                                                  >  2:132965/1‑150 (MQ=38)
                                                      gtgCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATcc                                                                                            <  2:236267/150‑1 (MQ=255)
                                                        gCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCca                                                                                          <  1:2337387/150‑1 (MQ=255)
                                                        gCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCca                                                                                          >  1:2178583/1‑150 (MQ=255)
                                                          aGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCcaca                                                                                        <  2:1881695/150‑1 (MQ=255)
                                                          aGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCcaca                                                                                        <  2:1882318/150‑1 (MQ=255)
                                                            ttGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAGCCCCTCGCCATCCCACAGa                                                                                      <  2:5798446/150‑1 (MQ=255)
                                                             tgatgaTTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGaa                                                                                     <  1:8274017/150‑1 (MQ=255)
                                                                 gatTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTc                                                                                 >  1:5413378/1‑150 (MQ=255)
                                                                   ttGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTccc                                                                               <  1:1147998/150‑1 (MQ=255)
                                                                      cgcATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTcccccc                                                                            >  1:2729251/1‑150 (MQ=255)
                                                                       gcATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCcca                                                                           >  2:4944474/1‑150 (MQ=255)
                                                                        cATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCccac                                                                          <  2:3650897/150‑1 (MQ=255)
                                                                             ggAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAAta                                                                     >  2:3812141/1‑149 (MQ=255)
                                                                             ggAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATc                                                                     >  2:7699001/1‑150 (MQ=255)
                                                                             ggAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATc                                                                     >  2:3887604/1‑150 (MQ=255)
                                                                             ggAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATc                                                                     >  1:3971013/1‑150 (MQ=255)
                                                                               aGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTg                                                                   <  1:5545677/150‑1 (MQ=255)
                                                                               aGCCGCGCTGAGTTTAGCGTCACGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTg                                                                   <  2:7109433/150‑1 (MQ=255)
                                                                                  cgcgCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTc                                                                <  2:5455473/150‑1 (MQ=255)
                                                                                    cgcTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTa                                                              <  2:8829934/150‑1 (MQ=255)
                                                                                       tGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCcca                                                                           >  1:96489/1‑134 (MQ=255)
                                                                                       tGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGtc                                                           >  1:3485056/1‑150 (MQ=255)
                                                                                       tGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGtc                                                           >  2:4738103/1‑150 (MQ=255)
                                                                                        gAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGtct                                                          <  1:3455958/150‑1 (MQ=255)
                                                                                          gTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCa                                                        <  2:4641043/150‑1 (MQ=255)
                                                                                            ttAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCAtt                                                      <  1:1964102/150‑1 (MQ=255)
                                                                                            ttAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCAtt                                                      <  1:4196810/150‑1 (MQ=255)
                                                                                               gCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAAt                                                   >  2:6561654/1‑150 (MQ=255)
                                                                                                 gTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCCATCTGGTCTAGTCTCATTAATGa                                                 >  2:6226203/1‑150 (MQ=255)
                                                                                                  tCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGAc                                                <  2:613711/150‑1 (MQ=255)
                                                                                                  tCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGAc                                                <  2:5369579/150‑1 (MQ=255)
                                                                                                    aGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCg                                              <  2:969723/150‑1 (MQ=255)
                                                                                                    aGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCg                                              >  1:3212165/1‑150 (MQ=255)
                                                                                                    aGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCg                                              >  1:7827882/1‑150 (MQ=255)
                                                                                                    aGGTGCCTGCCTTGGAACTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCCCCGCGCTCGATCTCACAGGCGCTGCAAAACTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCCATCTCGTCAAGTCTCATTAATGACCg                                              >  2:8523431/1‑150 (MQ=255)
                                                                                                     ggTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCACCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCgg                                             <  2:163979/150‑1 (MQ=255)
                                                                                                      gTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCGCACAGGCGCTGCACACCTTGCGGCGAACCCCTCGCCATCCCACCGAACCTCCCCCCACCACGCGGGTCTCGGCTCCGTCATGACCGGc                                            >  2:7203284/1‑150 (MQ=37)
                                                                                                       tgcctgccTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCa                                           >  1:4272802/1‑150 (MQ=255)
                                                                                                        gcctgcctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTa                                                              <  1:2776370/130‑1 (MQ=255)
                                                                                                        gcctgcctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTa                                                              >  2:2776370/1‑130 (MQ=255)
                                                                                                          ctgcctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCg                                        <  1:7699001/150‑1 (MQ=255)
                                                                                                           tgcctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCgg                                       >  2:7933419/1‑150 (MQ=255)
                                                                                                           tgcctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCgg                                       >  1:6342366/1‑150 (MQ=255)
                                                                                                              ctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGcc                                    >  2:795494/1‑150 (MQ=255)
                                                                                                               ttGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCa                                   >  1:1540833/1‑150 (MQ=255)
                                                                                                               ttGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCa                                   >  1:991683/1‑150 (MQ=255)
                                                                                                                 ggATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCCATCTGGTCTCGTCTCATTAATGACCGGCACCGGGCCAgg                                 >  2:6357591/1‑150 (MQ=255)
                                                                                                                 ggATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAgg                                 <  1:471420/150‑1 (MQ=255)
                                                                                                                   aTCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGa                               <  1:6051255/150‑1 (MQ=255)
                                                                                                                   aTCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGa                               <  1:4395920/150‑1 (MQ=255)
                                                                                                                        gcgcTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGc                          >  1:1691934/1‑150 (MQ=255)
                                                                                                                              gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                    <  1:6178384/150‑1 (MQ=255)
                                                                                                                              gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                    <  1:3721518/150‑1 (MQ=255)
                                                                                                                              gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                    <  1:5433105/150‑1 (MQ=255)
                                                                                                                              gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                    <  1:5955963/150‑1 (MQ=255)
                                                                                                                                tgtgagaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGGACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgt                  >  1:1841976/7‑150 (MQ=38)
                                                                                                                                gCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgt                  >  2:3248696/1‑150 (MQ=255)
                                                                                                                                gCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgt                  <  1:449562/150‑1 (MQ=255)
                                                                                                                                gCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgt                  <  2:8717366/150‑1 (MQ=255)
                                                                                                                                  acaAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAgg                                 >  1:8239253/1‑133 (MQ=255)
                                                                                                                                    aaaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTcc              >  1:508291/1‑150 (MQ=255)
                                                                                                                                      aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCgg            >  1:7031165/1‑150 (MQ=255)
                                                                                                                                        cGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGc          >  2:828744/1‑150 (MQ=255)
                                                                                                                                         gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGcc         >  2:3717938/1‑150 (MQ=255)
                                                                                                                                         gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGcc         <  1:1451108/150‑1 (MQ=255)
                                                                                                                                         gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCCTTCGGCGAACCCCTCGCCATCCAACAGAACCTCCCCCCACCACTCTGCTCTATTCTCATTAATGACCGGCACCGGGCCAGGGATGTACACTTGAGTGTCCGGGcc         <  2:5669664/150‑1 (MQ=39)
                                                                                                                                               cgcccgcGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCa   >  1:7176217/1‑150 (MQ=255)
                                                                                                                                               cgcccgcGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCa   >  2:1253589/1‑150 (MQ=255)
                                                                                                                                                gcccgcGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCAc  <  2:7557142/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGGGTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCAC  >  NC_002947/1683144‑1683437

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