Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,863,120 C→T 6.6% intergenic (‑104/+133) ligA ← / ← PP_4275 DNA ligase/cell division protein ZipA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,863,1200CT 6.6% 474.6 / 14.8 154intergenic (‑104/+133)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base C (90/52);  new base T (4/6);  total (95/59)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 163 total aligned reads displayed.

GATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGT  >  NC_002947/4862971‑4863267
                                                                                                                                                      |                                                                                                                                                     
gATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCACTccc                                                                                                                                                       >  1:2182762/1‑150 (MQ=255)
gATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTccc                                                                                                                                                       >  2:2589291/1‑150 (MQ=255)
 aTCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCca                                                                                                                                                      <  2:5224626/150‑1 (MQ=255)
 aTCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGGTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCca                                                                                                                                                      <  1:2359406/150‑1 (MQ=255)
  tCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCcac                                                                                                                                                     >  2:351243/1‑150 (MQ=255)
      gCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAgg                                                                                                                                                 <  2:4575720/150‑1 (MQ=255)
          ggCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGt                                                                                                                                             >  1:447695/1‑150 (MQ=255)
           gCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTc                                                                                                                                            >  2:3966157/1‑150 (MQ=255)
            cACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAACGGCGTcc                                                                                                                                           >  2:4573518/1‑149 (MQ=255)
              cGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCACCCACAAGGGCggccca                                                                                                                                         >  1:3221542/1‑145 (MQ=255)
                      gAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAg                                                                                                                                 <  1:1765905/150‑1 (MQ=255)
                          aTACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCg                                                                                                                            >  1:1132315/1‑150 (MQ=255)
                          aTACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCg                                                                                                                            <  2:3838747/150‑1 (MQ=255)
                          aTACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACCGCCCTAGGGcat                                                                                                                            <  2:2188840/150‑3 (MQ=255)
                           tACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCgg                                                                                                                           <  2:5018170/150‑1 (MQ=255)
                            aCGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGt                                                                                                                          <  2:1052618/150‑1 (MQ=255)
                              gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                        <  1:3015326/150‑1 (MQ=255)
                               ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                       <  2:1325867/150‑1 (MQ=255)
                               ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                       <  1:2025180/150‑1 (MQ=255)
                               ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                       <  1:5105596/150‑1 (MQ=255)
                                       cggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGt                                                                                                               <  1:5245636/150‑1 (MQ=255)
                                         gTCATGGGATTTGTGTTCTCTTGCATTGCAATAGAGCAGCTTGGGGTGGTGGGTGGGGTTGCAGGTGTCTGTTGACTGTGCCAGCCTCTTGGCGGGTTTAGCCGCTGCGAGAAGGGGGTCAGTTTGCCTTGGGCGGTGTTGTTGCTGttg                                                                                                             >  2:2960619/1‑148 (MQ=255)
                                         gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                                             <  1:5170034/150‑1 (MQ=255)
                                         gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                                             <  1:5023431/150‑1 (MQ=255)
                                         gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAg                                                                                                             >  2:2960647/1‑150 (MQ=255)
                                            aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTACCACAACGGCGTCACATTCCCAAGGGCGGTGATGTACCAGTAGGAg                                                                                                          >  1:3517425/1‑150 (MQ=255)
                                             tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                         >  1:2515065/1‑150 (MQ=255)
                                             tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                         <  2:3610661/150‑1 (MQ=255)
                                              gggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCg                                                                                                       <  1:1472658/150‑1 (MQ=255)
                                               ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                      <  2:2515065/150‑1 (MQ=255)
                                                gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCggg                                                                                                     <  1:3916244/150‑1 (MQ=14)
                                                gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                     <  2:4856083/150‑1 (MQ=255)
                                                gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                     <  2:3848958/150‑1 (MQ=255)
                                                gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCGCCCACAACGGCCTCACAAGCCCAAGGGCGGGGATGTACCGGTAGGAGCggg                                                                                                     >  1:3294550/1‑150 (MQ=255)
                                                 aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                    <  1:138900/150‑1 (MQ=14)
                                                 aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                    <  2:1807067/150‑1 (MQ=14)
                                                 aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGt                                                                                                    <  2:3221918/150‑1 (MQ=255)
                                                  tttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGtt                                                                                                   <  1:2986704/150‑1 (MQ=14)
                                                       cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAccc                                                                                              <  2:1879393/150‑1 (MQ=14)
                                                       cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGGAGGAGCCGGTTTAccc                                                                                              >  2:938547/1‑150 (MQ=255)
                                                        ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACACGCGCACCCCACGGCGGCACAAGCCCACGGGCCGTGATTTAACCGTTGGATCCGGTTTAaccc                                                                                             >  1:1105550/1‑149 (MQ=255)
                                                         tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcgc                                                                                            >  2:4950666/1‑150 (MQ=255)
                                                         tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCGAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCcgc                                                                                            >  2:1843698/1‑150 (MQ=255)
                                                           tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGTATACCCCGCCCACAAGGGCGGGACAAGCCCGAGGGCGGGGGTGTGCCGGGTGGTGCGGGTGTTACCGCGa                                                                                          >  1:636720/1‑150 (MQ=255)
                                                           tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                          >  2:4421977/1‑150 (MQ=255)
                                                            ctTGCAAAGCAAAAGAGCAGCTGGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAGGCCCTAGGGCGGTGTTGTACCTGTAGGAGCGGGTTTACCCGCGaa                                                                                         <  2:3961019/150‑1 (MQ=255)
                                                             ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATTCCCTAGGGCGGTGAAGTACCTGTAGGAGCGGGTTTACCCGCGAag                                                                                        >  2:1119677/1‑150 (MQ=255)
                                                               gctttgcaagAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAagag                                                                                      >  2:2776503/6‑150 (MQ=255)
                                                               gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCCCGAagag                                                                                      >  1:4482330/1‑150 (MQ=255)
                                                                caaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCTGTCGGAGCGGGTTTACCCGCGAAGCgg                                                                                     >  2:3728179/1‑150 (MQ=255)
                                                                 aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAACCCCCGCCCACAACGGCGGCACAAGCCCAAAGGCGGGGATGTACCGGTAGGAGCGGGGTTAACCCCGAAGAGGc                                                                                    >  1:2834589/1‑150 (MQ=255)
                                                                    gcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCTGGTGTCTGTTGACTGGGCCAGCCTCTTCGCGGTTAAACCCGCTCCCACAAGGGCGTCACATGCCCTCGGGCGGTGATGTACCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCgg                                                                                 <  2:4826177/150‑1 (MQ=38)
                                                                    gcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCgg                                                                                 >  2:809346/1‑150 (MQ=255)
                                                                     caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAACCCCCGCCCCCAACGGCGTCAAAAGCCCAAGGGCGGGGATGTACCTGGAGGAGCGGGGTTTAACCCGAAGAGGCCCGc                                                                                >  1:4688467/1‑150 (MQ=255)
                                                                      naaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCCCACGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGca                                                                               >  1:5380945/2‑150 (MQ=255)
                                                                      aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTGAACCCCCCCCCCCCAGGGGGGCACAAGCCCCCGGGCGGGGGTGGACCCGTGGGGGCGGGGGTGACCGCCAAGAGggcggcc                                                                               >  1:656113/1‑149 (MQ=255)
                                                                      aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGca                                                                               >  2:2071994/1‑150 (MQ=255)
                                                                      aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGTGATGTACCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGca                                                                               >  2:4278432/1‑150 (MQ=255)
                                                                       aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcac                                                                              >  1:4432688/1‑150 (MQ=255)
                                                                       aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCGCCCACAAGGGCGTCACAAGCCCTAGGGCGGGGCTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcac                                                                              >  1:4432581/1‑150 (MQ=255)
                                                                       aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCACAAGGGCGTGCCCAGCCCGTGGGGGGTGGTGTACCTGTAGGAGGGGGTTTTCCCCCGAAGAGGCCGCcac                                                                              >  1:4753578/1‑150 (MQ=255)
                                                                       aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCACAAGGGCGTCACAAGCCCAAGGGCGGTGCTGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcac                                                                              >  1:3946110/1‑150 (MQ=38)
                                                                                gtgggggTGGTGGGGGGGTTGCCTGGTGTCTGTTGAGTGGGCCCGCCTCTTCGGGGTTAAACCCGCCCCAACCGGGGCGACACTGGCCCTGGGGGGGTTTTGTACGTTTGGGCGGGGTTTTCCCCGCGAAGAGGCCGGCACAGGCAacga                                                                     <  2:1105550/147‑1 (MQ=255)
                                                                                   ggggTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAcc                                                                                               <  1:4939437/121‑1 (MQ=255)
                                                                                      gtggtgCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATaa                                                               >  2:1112186/1‑149 (MQ=255)
                                                                                                  gTTCCAGGTGTCTGTTGACTGGGCCAGCCTCTGCGCGGGTAACCCCGCTCCCACACGGGCGTCACAGGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCa                                                   <  2:5353530/150‑1 (MQ=38)
                                                                                                   ttCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAg                                                  >  2:380411/1‑150 (MQ=255)
                                                                                                   ttCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAg                                                  >  2:380257/1‑150 (MQ=255)
                                                                                                   ttCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCGCCCACAAGGGGGTGACAAGCCCTAGGGCGGTGATGTACCAGTAGGAGCGGGTTTACACGCGAAGAGGCCGGCACCGGCCACGACGATACCCGAACCCGCCAg                                                  >  2:5015977/1‑150 (MQ=38)
                                                                                                      cAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCACCCACAAGGGCCTCACAAGCCCTAGGGCCGTGATGTACCGGTAGGAGCGGGTTTACCCCCGAAGAGGCCGGCACAGGCAACCACGATACCCGAACccgcccgccc                                               >  2:3408642/1‑148 (MQ=37)
                                                                                                         gtgtCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATg                                                                                                                      <  1:2308434/77‑1 (MQ=255)
                                                                                                         gtgtCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATg                                                                                                                      >  2:2308434/1‑77 (MQ=255)
                                                                                                             gtGTTGAGTGGGCCAGCCTCTTCGGGGTTAAACCCGCGCCCACCCGGGCGCCACTGGCCCTCGGGGGTTGTTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCa                                        <  2:5182292/149‑1 (MQ=25)
                                                                                                              tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCCGCACGCCCTCAg                                       >  2:3586521/1‑150 (MQ=255)
                                                                                                              tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAg                                       >  2:1262574/1‑150 (MQ=255)
                                                                                                              tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCACCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGCACAGGCCCGCCCCGGCAACGACGATACCCGAACCGGCCAGCACGCCCACAg                                       >  2:4707107/1‑150 (MQ=38)
                                                                                                                ttGACTGTCCCCCCCTCTTCGCGGTTAAACCCGCACCCAAACGGGGGAAACTTCCCCTTGGGGGTTTTTTTACCTGTTGGAGGGTTTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGgc                                     <  2:2041170/150‑1 (MQ=255)
                                                                                                                         cccgCCTTTTCGCGTTTAACCCCGCTCCCCAAGGGGCGTCCCTGGCCCTCGGGCGTGGTTGTACCTGTCGGAGGGGTTTTACCGGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                            <  1:4532465/147‑1 (MQ=255)
                                                                                                                             ccTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGGCACATGCCCGAGGGCGGGGATGTACCGGTAGGAGCGGGGGTACCCGCGCAGAGGCCGGCACCGGCCACGACGATACCCGAACCGGCCAGCACGCCCGCAGGCGGCGCTTTCGGCt                        >  2:5345196/1‑150 (MQ=255)
                                                                                                                               tctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATTCCCGAGGGCGGTGTTGTACCTGTAGGAGCGGGTTTACCCCCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCACAGGCTGCTCGTTCGTCTgg                      >  2:1492684/1‑150 (MQ=38)
                                                                                                                               tctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTgg                      >  2:3631964/1‑150 (MQ=255)
                                                                                                                               tctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTgg                      >  2:2213076/1‑150 (MQ=255)
                                                                                                                               tctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACCATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTgg                      >  1:2229900/1‑150 (MQ=255)
                                                                                                                               tctGCGCGGGAACACCCGCTCCCACACGGGCGTCACATCCCCTTGGGGGGTGTGGTACGTGTTGGCGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTgg                      <  2:2834589/150‑1 (MQ=37)
                                                                                                                                 ttCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTATGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTAGTCGGGAt                    >  2:2490477/1‑150 (MQ=255)
                                                                                                                                    gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                 >  1:1140630/1‑150 (MQ=255)
                                                                                                                                    gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCGGGATCTg                 >  1:4183843/1‑150 (MQ=255)
                                                                                                                                       gggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTc               >  1:1372937/1‑150 (MQ=255)
                                                                                                                                       gggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTc               >  2:1668817/1‑150 (MQ=255)
                                                                                                                                        ggaacccccGCTCCCACCCGGGCGCCACGGGCCCGGGGGGGGTGTTGCCCGTGTCGGCGGGGTTTCCCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa              <  2:656113/144‑1 (MQ=255)
                                                                                                                                        ggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa              >  2:34345/1‑150 (MQ=255)
                                                                                                                                         gTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGAAACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAg             >  1:1758815/1‑150 (MQ=255)
                                                                                                                                             aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt         >  2:1119734/1‑150 (MQ=14)
                                                                                                                                               ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc       >  2:787639/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc       >  1:1099513/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc       >  1:1210814/1‑150 (MQ=255)
                                                                                                                                               ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc       >  1:3284044/1‑150 (MQ=255)
                                                                                                                                                  cgcccACCAGGGCGTCACCTGCCCTTGGGGGGTGATGTACGTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTgg    <  2:4308116/148‑1 (MQ=38)
                                                                                                                                                   tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg   >  2:5280309/1‑150 (MQ=255)
                                                                                                                                                   tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg   >  1:2696610/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGT  >  NC_002947/4862971‑4863267

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