Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,863,124 A→G 8.4% intergenic (‑108/+129) 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,1240AG 8.4% 448.7 / 19.9 158intergenic (‑108/+129)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base A (91/45);  new base G (4/9);  total (97/61)
Fisher's exact test for biased strand distribution p-value = 1.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

GAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGC  >  NC_002947/4862975‑4863271
                                                                                                                                                      |                                                                                                                                                    
gAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACaa                                                                                                                                                      >  1:1638807/1‑150 (MQ=255)
  gCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAgg                                                                                                                                                    >  1:5018170/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)
          cGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCACCCACAAGGGCggccca                                                                                                                                            >  1:3221542/1‑145 (MQ=255)
                tCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCt                                                                                                                                      <  2:2635605/150‑1 (MQ=255)
                  gAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAg                                                                                                                                    <  1:1765905/150‑1 (MQ=255)
                   aGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAgg                                                                                                                                   >  2:453055/1‑150 (MQ=255)
                   aGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAgg                                                                                                                                   <  2:1537724/150‑1 (MQ=255)
                   aGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAgg                                                                                                                                   >  2:3502077/1‑150 (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)
                          gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGcatgtga                                                                                                                           >  1:222503/1‑144 (MQ=21)
                          gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                           <  1:3015326/150‑1 (MQ=255)
                          gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                           <  1:1882365/150‑1 (MQ=255)
                           ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                          <  1:5105596/150‑1 (MQ=255)
                           ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                          <  2:1325867/150‑1 (MQ=255)
                           ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                          <  2:1325698/150‑1 (MQ=255)
                           ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                          <  2:2228670/150‑1 (MQ=255)
                                  gcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTg                                                                                                                   <  1:4262312/150‑1 (MQ=255)
                                     gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                                                <  1:5170034/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:3848958/150‑1 (MQ=255)
                                            gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                         <  2:4856083/150‑1 (MQ=255)
                                            gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCGCCCACAACGGCCTCACAAGCCCAAGGGCGGGGATGTACCGGTAGGAGCggg                                                                                                         >  1:3294550/1‑150 (MQ=255)
                                             aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                        <  1:3403243/150‑1 (MQ=14)
                                             aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                        <  2:1807067/150‑1 (MQ=14)
                                             aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                        <  1:138900/150‑1 (MQ=14)
                                              tttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGtt                                                                                                       <  1:2986704/150‑1 (MQ=14)
                                               ttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGttt                                                                                                      <  2:1132315/150‑1 (MQ=255)
                                                   cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAccc                                                                                                  <  2:1879393/150‑1 (MQ=14)
                                                   cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGGAGGAGCCGGTTTAccc                                                                                                  >  2:938547/1‑150 (MQ=255)
                                                    ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACACGCGCACCCCACGGCGGCACAAGCCCACGGGCCGTGATTTAACCGTTGGATCCGGTTTAaccc                                                                                                 >  1:1105550/1‑149 (MQ=255)
                                                     tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCACCCACACGGGCGTCACATGCCCGAGGGCGGTGATGTAACAGTAGGAGCGGGGTTAccccc                                                                                                >  1:4596654/1‑148 (MQ=255)
                                                       tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGTATACCCCGCCCACAAGGGCGGGACAAGCCCGAGGGCGGGGGTGTGCCGGGTGGTGCGGGTGTTACCGCGa                                                                                              >  1:636720/1‑150 (MQ=255)
                                                       tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                              >  2:4818295/1‑150 (MQ=255)
                                                       tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                              >  2:4421977/1‑150 (MQ=255)
                                                       tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGCTGTACCGGTAGGAGCGGGTTTACCCGCGa                                                                                              >  2:5231210/1‑150 (MQ=255)
                                                        ctTGCAAAGCAAAAGAGCAGCTGGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAGGCCCTAGGGCGGTGTTGTACCTGTAGGAGCGGGTTTACCCGCGaa                                                                                             <  2:3961019/150‑1 (MQ=255)
                                                         ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAag                                                                                            <  2:3242638/150‑1 (MQ=255)
                                                           gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCCCGAagag                                                                                          >  1:4482330/1‑150 (MQ=255)
                                                             aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAACCCCCGCCCACAACGGCGGCACAAGCCCAAAGGCGGGGATGTACCGGTAGGAGCGGGGTTAACCCCGAAGAGGc                                                                                        >  1:2834589/1‑150 (MQ=255)
                                                                gcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCgg                                                                                     >  2:809346/1‑150 (MQ=255)
                                                                 caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTATGGCGGTGATGTACCGGTAGGAGCGGGTTTAACCGCGAAGAGGCCGGc                                                                                    >  2:2924829/1‑150 (MQ=255)
                                                                  naaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCCCACGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGca                                                                                   >  1:5380945/2‑150 (MQ=255)
                                                                  aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTGAACCCCCCCCCCCCAGGGGGGCACAAGCCCCCGGGCGGGGGTGGACCCGTGGGGGCGGGGGTGACCGCCAAGAGggcggcc                                                                                   >  1:656113/1‑149 (MQ=255)
                                                                  aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcc                                                                                   >  2:4278694/1‑149 (MQ=255)
                                                                  aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGTGATGTACCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGca                                                                                   >  2:4278432/1‑150 (MQ=255)
                                                                   aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcac                                                                                  >  1:4432688/1‑150 (MQ=255)
                                                                   aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCACAAGGGCGTGCCCAGCCCGTGGGGGGTGGTGTACCTGTAGGAGGGGGTTTTCCCCCGAAGAGGCCGCcac                                                                                  >  1:4753578/1‑150 (MQ=255)
                                                                   aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCACAAGGGCGTCACAAGCCCAAGGGCGGTGCTGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcac                                                                                  >  1:3946110/1‑150 (MQ=38)
                                                                    aaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCTGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcaca                                                                                 >  2:5013574/1‑150 (MQ=255)
                                                                               ggggTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAcc                                                                                                   <  1:4939437/121‑1 (MQ=255)
                                                                               ggggTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAcc                                                                                                   >  2:4939437/1‑121 (MQ=255)
                                                                                  gtggtgCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATaa                                                                   >  2:1112186/1‑149 (MQ=255)
                                                                                               ttCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAg                                                      >  2:380411/1‑150 (MQ=255)
                                                                                               ttCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAg                                                      >  2:380257/1‑150 (MQ=255)
                                                                                                     gtgtCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATg                                                                                                                         >  2:2308434/1‑77 (MQ=255)
                                                                                                     gtgtCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATg                                                                                                                         <  1:2308434/77‑1 (MQ=255)
                                                                                                         gtGTTCACTGGGCCCGCCTCTGCGGGGTTAAACCCGCCCCCACAGGGGCTACACATGCCCTCGGGGGGTTTTGTCCGTGTTGGCGGGGTTTTCCCCGCCAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCa                                            <  2:3789027/149‑1 (MQ=255)
                                                                                                          tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCCGCACGCCCTCAg                                           >  2:3586521/1‑150 (MQ=255)
                                                                                                          tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCGCCCACAAGGGCGGCACATGCCCGAGGGCGGGGGTGTACCGGGAGGAGCGGGGTTACCCGCGCAGAGGCCGGCCCAGGCCACGCCGCAACCCGAACCGGCCAGCACGCCCGCAg                                           >  2:1098691/1‑150 (MQ=25)
                                                                                                          tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCACCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGCACAGGCCCGCCCCGGCAACGACGATACCCGAACCGGCCAGCACGCCCACAg                                           >  2:4707107/1‑150 (MQ=38)
                                                                                                            ttGACTGTCCCCCCCTCTTCGCGGTTAAACCCGCACCCAAACGGGGGAAACTTCCCCTTGGGGGTTTTTTTACCTGTTGGAGGGTTTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGgc                                         <  2:2041170/150‑1 (MQ=255)
                                                                                                                     cccgCCTTTTCGCGTTTAACCCCGCTCCCCAAGGGGCGTCCCTGGCCCTCGGGCGTGGTTGTACCTGTCGGAGGGGTTTTACCGGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                                <  1:4532465/147‑1 (MQ=255)
                                                                                                                     ccAGCCTCTTCGCGGGTAAACCCCCGCCCACAACGGCCTCACATGCCCCAGGGCGGGGATGTACCGGGAGGCGCCGGGGTACCCGCGAAGAGGCCCGCCCCCGCAACCCCCAACACCGCACCGGCCCGCCCCCCCACAggcggcgcgtgc                                >  2:5138286/1‑144 (MQ=255)
                                                                                                                         ccTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGGCACATGCCCGAGGGCGGGGATGTACCGGTAGGAGCGGGGGTACCCGCGCAGAGGCCGGCACCGGCCACGACGATACCCGAACCGGCCAGCACGCCCGCAGGCGGCGCTTTCGGCt                            >  2:5345196/1‑150 (MQ=255)
                                                                                                                           tctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTgg                          >  2:3631964/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)
                                                                                                                                 cgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGt                    >  1:2589715/1‑150 (MQ=255)
                                                                                                                                    ggaacccccGCTCCCACCCGGGCGCCACGGGCCCGGGGGGGGTGTTGCCCGTGTCGGCGGGGTTTCCCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                  <  2:656113/144‑1 (MQ=255)
                                                                                                                                    ggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                  >  2:34345/1‑150 (MQ=255)
                                                                                                                                    ggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                  <  2:1770828/150‑1 (MQ=255)
                                                                                                                                     gTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAg                 >  2:1846280/1‑150 (MQ=255)
                                                                                                                                      tAAACCCGCTCCCACCCGGGCGTCACCTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGAGACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgc                <  2:5361721/150‑1 (MQ=255)
                                                                                                                                         aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt             >  2:1119734/1‑150 (MQ=14)
                                                                                                                                         aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTGTCGTCTGGATCTGTCAGCGCt             >  1:52675/1‑150 (MQ=14)
                                                                                                                                          cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt            >  2:3098937/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc           >  1:1210814/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc           >  2:787639/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc           >  1:1099513/1‑150 (MQ=255)
                                                                                                                                             gCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTg         >  1:1256361/1‑150 (MQ=255)
                                                                                                                                               tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg       >  2:5280309/1‑150 (MQ=255)
                                                                                                                                               tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg       >  2:2160229/1‑150 (MQ=255)
                                                                                                                                                cccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGt      >  1:2263428/1‑150 (MQ=255)
                                                                                                                                                   acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg   >  1:509860/1‑150 (MQ=255)
                                                                                                                                                    caAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgc  <  1:3098937/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGC  >  NC_002947/4862975‑4863271

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