Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,683,290 C→G 13.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,2900CG13.9% 213.5 / 26.0 79intergenic (+57/‑174)recJ/PP_1478ssDNA‑specific exonuclease RecJ/xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base C (40/28);  new base G (6/5);  total (46/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

TGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCAC  >  NC_002947/1683148‑1683437
                                                                                                                                              |                                                                                                                                                   
tGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCgcgg                                                                                                                                              <  2:949637/150‑1 (MQ=255)
    gTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                          <  2:2267503/150‑1 (MQ=255)
    gTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                          >  2:939734/1‑150 (MQ=255)
        gCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCTCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTTGATCTTGCGCTGGGGCACAAAGCGAGCGCCGCGCGGGCGGCGCTCTAt                                                                                                                                      <  2:1696793/150‑1 (MQ=255)
        gCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                                      <  1:2170318/150‑1 (MQ=255)
              aaGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTcaca                                                                                                                                >  2:1134439/1‑150 (MQ=255)
                gCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAAGCGAGCTCCGCCCGCGCGGCGCTCGATCTCACAgg                                                                                                                              >  2:2310751/1‑150 (MQ=255)
                        tgAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCaa                                                                                                                      >  2:2349625/1‑150 (MQ=255)
                          aaCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCgcggc                                                                                                                                             >  2:1918075/1‑125 (MQ=255)
                          aaCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCgcggc                                                                                                                                             <  1:1918075/125‑1 (MQ=255)
                          aaCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAAc                                                                                                                    <  2:1710172/150‑1 (MQ=255)
                                      ggTAATGAGAGTGTGCAGTTGATGATTGCGCATACGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAAc                                                                                                        >  2:1807427/1‑150 (MQ=255)
                                      gatcgggagaGTGGGGAGGTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCGGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAAc                                                                                                        <  1:2306732/144‑1 (MQ=38)
                                          aTGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                       <  2:1529023/115‑1 (MQ=255)
                                          aTGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                       >  1:1529023/1‑115 (MQ=255)
                                            gagaGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTGGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCg                                                                                                  >  1:1207529/1‑150 (MQ=255)
                                             agagTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGc                                                                                                 >  1:2164027/1‑150 (MQ=255)
                                             agagTGTGCAGTTGATGATTGCGCATATAGAGCAGCCGTGAGCTGAATGTCAAACTGGTGCCTTGGATCTATGACTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATATCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGc                                                                                                 >  1:421468/1‑150 (MQ=21)
                                                      aGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCcaca                                                                                        <  2:1559285/150‑1 (MQ=255)
                                                        ttGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                         >  1:415001/1‑99 (MQ=255)
                                                        ttGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                         <  2:415001/99‑1 (MQ=255)
                                                        ttGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTCCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                           <  2:431026/97‑1 (MQ=255)
                                                        ttGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTCCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                           >  1:431026/1‑97 (MQ=255)
                                                          gatgatTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTca                                                                                                                                  >  2:1516787/1‑104 (MQ=255)
                                                          gatgatTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTca                                                                                                                                  <  1:1516787/104‑1 (MQ=255)
                                                          gatgatTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAAc                                                                                    >  1:73993/1‑150 (MQ=255)
                                                              atTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTcc                                                                                >  1:2083083/1‑150 (MQ=255)
                                                                 gcgcATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTGGCGCTGGGGCACAAAGGGAGCGCCGCCCGCGCGGCGCTCGGTCTCACGGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGGACCTccccc                                                                             >  2:210481/1‑150 (MQ=255)
                                                                   gcATATGGAGCCGCGCTGAGTCTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCGCAACTCGAGCGCCGCCGGCGCGGCGCTCGTTTTGACAGGCGCTGCAAAGCCTTCGGCGGACCCCTAGACATCCCTCCGCACCTCGTCcca                                                                           >  2:415613/1‑150 (MQ=255)
                                                                     atatGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAAc                                                                                                                    <  2:1760941/107‑1 (MQ=255)
                                                                     atatGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAAc                                                                                                                    >  1:1760941/1‑107 (MQ=255)
                                                                      tatGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTGGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCccacca                                                                        >  2:1394578/1‑150 (MQ=255)
                                                                           aGCCGCGCTGAGGTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACATCCCCCCACCAATCTg                                                                   <  2:1822410/150‑1 (MQ=255)
                                                                             ccGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGt                                                                 <  2:17161/150‑1 (MQ=255)
                                                                             ccGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGt                                                                 <  1:1134439/150‑1 (MQ=255)
                                                                                cgcTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGGCGCCGGCCTGGCGCTCGATGTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTa                                                              <  1:2349625/150‑1 (MQ=255)
                                                                                   tGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGtc                                                           >  2:1658186/1‑150 (MQ=255)
                                                                                       tttAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAAc                                                                                                        >  1:2208277/1‑101 (MQ=255)
                                                                                       tttAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAAc                                                                                                        <  2:2208277/101‑1 (MQ=255)
                                                                                       tttAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGAAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCAt                                                       >  2:868813/1‑150 (MQ=255)
                                                                                           gCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAAt                                                   <  2:642747/150‑1 (MQ=255)
                                                                                             gTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGt                                                            >  1:577760/1‑139 (MQ=255)
                                                                                             gTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGt                                                            <  2:577760/139‑1 (MQ=255)
                                                                                              tCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                         >  2:111671/1‑61 (MQ=21)
                                                                                              tCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                         <  1:111671/61‑1 (MQ=21)
                                                                                              tCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGAGGCCGTCTTATGGCGGACCCCTCGCCATCCCACAGCAGTCACTCACTGAAATCTGGTCTGGTCTCATTAATGAc                                                >  1:655833/1‑150 (MQ=19)
                                                                                              tCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGc                                                                                                                        >  1:222238/1‑78 (MQ=255)
                                                                                              tCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGc                                                                                                                        <  2:222238/78‑1 (MQ=255)
                                                                                                     cctgcctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCAcc                                         >  1:405131/1‑150 (MQ=255)
                                                                                                       tgcctTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCgg                                       >  2:963003/1‑150 (MQ=255)
                                                                                                           ttGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTc                                                                                 >  2:1143019/1‑104 (MQ=255)
                                                                                                           ttGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTc                                                                                 <  1:1143019/104‑1 (MQ=255)
                                                                                                            tGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAg                                  <  1:415613/150‑1 (MQ=255)
                                                                                                             ggATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGTCACCGGGCCAgg                                 <  2:2083083/150‑1 (MQ=255)
                                                                                                              gATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAAcccc                                                                                                     <  2:809701/81‑1 (MQ=255)
                                                                                                              gATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAAcccc                                                                                                     >  1:809701/1‑81 (MQ=255)
                                                                                                              gATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAggg                                >  1:321048/1‑150 (MQ=255)
                                                                                                               aTCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCGCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGa                               >  2:1617018/1‑150 (MQ=255)
                                                                                                                  ttGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGAtgt                            <  1:958952/150‑1 (MQ=255)
                                                                                                                  ttGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGAtgt                            >  2:2090578/1‑150 (MQ=255)
                                                                                                                          gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATc                                                                                             >  2:1568742/1‑77 (MQ=255)
                                                                                                                          gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATc                                                                                             <  1:1568742/77‑1 (MQ=255)
                                                                                                                          gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                    <  1:1658186/150‑1 (MQ=255)
                                                                                                                          gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                    <  1:1174049/150‑1 (MQ=255)
                                                                                                                          gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTCTCGGCGAACCCCTCGCCATc                                                                                             <  1:1568650/77‑1 (MQ=255)
                                                                                                                          gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTCTCGGCGAACCCCTCGCCATc                                                                                             >  2:1568650/1‑77 (MQ=255)
                                                                                                                            tgtgagaTCGAGCGCCGCGCGGGCGGCGCTCGATTTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgt                  >  1:578994/7‑150 (MQ=39)
                                                                                                                            gCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgt                  >  2:112161/1‑150 (MQ=255)
                                                                                                                             gtgagaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgtg                 >  1:1388493/6‑150 (MQ=255)
                                                                                                                              acaAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATcc                                                                                            <  2:2390463/74‑1 (MQ=255)
                                                                                                                              acaAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATcc                                                                                            >  1:2390463/1‑74 (MQ=255)
                                                                                                                              acaAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgtgt                >  1:1588813/1‑150 (MQ=255)
                                                                                                                                aaaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTTCACTTGAGTGTcc              >  2:913911/1‑150 (MQ=255)
                                                                                                                                 aaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACGTGAGTGTCCg             >  2:1620808/1‑150 (MQ=255)
                                                                                                                                  aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCgg            >  1:1518226/1‑150 (MQ=255)
                                                                                                                                   tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCggg           >  2:297737/1‑150 (MQ=255)
                                                                                                                                      aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAAt                                                                      <  2:494780/88‑1 (MQ=255)
                                                                                                                                      aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAAt                                                                      >  1:494780/1‑88 (MQ=255)
                                                                                                                                      aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCg        >  1:1179251/1‑150 (MQ=255)
                                                                                                                                            gcccccGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCAc  <  1:2310751/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCAC  >  NC_002947/1683148‑1683437

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