Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,863,129 G→A 7.3% intergenic (‑113/+124) 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,1290GA 7.3% 705.4 / 29.7 238intergenic (‑113/+124)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base G (127/89);  new base A (2/15);  total (129/109)
Fisher's exact test for biased strand distribution p-value = 1.84e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

Only 100 of 240 total aligned reads displayed.

GGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGC  >  NC_002947/4862981‑4863278
                                                                                                                                                    |                                                                                                                                                      
ggCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGt                                                                                                                                                       <  1:7688935/150‑1 (MQ=255)
ggCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCGCCCACAAGGGCGt                                                                                                                                                       >  1:2719906/1‑150 (MQ=255)
 gCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTc                                                                                                                                                      <  2:2521317/150‑1 (MQ=255)
           cGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTa                                                                                                                                            <  2:7180115/150‑1 (MQ=255)
             aGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAgg                                                                                                                                          <  1:655523/150‑1 (MQ=255)
             aGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAgg                                                                                                                                          <  1:5638568/150‑1 (MQ=255)
                 tACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCgg                                                                                                                                     >  2:4889568/1‑150 (MQ=255)
                  aCGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGt                                                                                                                                    <  1:1620523/150‑1 (MQ=255)
                  aCGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGt                                                                                                                                    <  2:3602572/150‑1 (MQ=255)
                   cGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTg                                                                                                                                   <  2:7896674/150‑1 (MQ=255)
                    gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGt                                                                                                                                    >  2:420719/1‑148 (MQ=255)
                    gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGt                                                                                                                                    <  1:420719/148‑1 (MQ=255)
                    gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                                  <  2:6452845/150‑1 (MQ=255)
                     ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                                 <  2:4985320/150‑1 (MQ=255)
                     ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                                 <  2:7042405/150‑1 (MQ=255)
                     ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                                 <  2:538603/150‑1 (MQ=255)
                       tttCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGt                                                                                                                               <  2:6742548/150‑1 (MQ=255)
                       tttCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGGGATGt                                                                                                                               >  1:4434709/1‑150 (MQ=255)
                         tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                             <  2:8434926/150‑1 (MQ=255)
                          cggcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGGGCCATTCCCTAGGGCGGTGGTGTAcc                                                                                                                            >  1:526622/1‑150 (MQ=255)
                           ggcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCt                                                                                                                           <  1:8795200/150‑1 (MQ=255)
                             cggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGt                                                                                                                         <  1:2558835/150‑1 (MQ=255)
                               gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                                                       <  2:3377059/150‑1 (MQ=255)
                                  aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGgaggag                                                                                                                    >  2:8855564/1‑150 (MQ=255)
                                   tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCACCCACAAGGGCGTCACATGCCCAAGGGCGGTGATGTACCGGTAGGAGc                                                                                                                   >  2:2248823/1‑150 (MQ=255)
                                   tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAAACCCACCCACAAAGGCGTCACAATCCCAAGGGCGGGGATGTAACAGTATGAGc                                                                                                                   >  1:4649316/1‑150 (MQ=255)
                                     ggaggtGTCTTCTCTTGCAAAGCAAAAGAGCAGCGTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGGTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                 <  2:478728/145‑1 (MQ=255)
                                     ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTGAACCCCCACCCACAACGGCGGCAAATGCCCACGGGCGGGGGTGTAACTGGTGGAGccg                                                                                                                 >  1:8420847/1‑148 (MQ=255)
                                     ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                 <  1:9138867/150‑1 (MQ=255)
                                     ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                 >  2:43006/1‑150 (MQ=255)
                                      gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                <  2:7165788/150‑1 (MQ=255)
                                      gATTTGTATTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCggg                                                                                                                <  1:7592946/150‑1 (MQ=14)
                                       aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGt                                                                                                               <  2:6577975/150‑1 (MQ=255)
                                         ttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGttt                                                                                                             <  2:312793/150‑1 (MQ=14)
                                          tGTCTTCTCTTGCAATGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGGCCTATGGCGGTGTTGTACCTGTAGGAGCGGGTTTa                                                                                                            >  2:1013542/1‑150 (MQ=255)
                                          tGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTa                                                                                                            <  1:3529670/150‑1 (MQ=255)
                                            tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAcc                                                                                                          <  1:388348/150‑1 (MQ=14)
                                            tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAcc                                                                                                          >  2:5039004/1‑150 (MQ=255)
                                            tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCAAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAcc                                                                                                          >  1:9013753/1‑150 (MQ=255)
                                             cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAccc                                                                                                         <  2:2720291/150‑1 (MQ=255)
                                             cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTCCCACAACGGCGTCACAATCCCAAGGGCGGTGATGTAACAGTAGGAGCGGGTTTAacc                                                                                                         >  1:7037423/1‑150 (MQ=255)
                                              ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcg                                                                                                        >  2:6822529/1‑150 (MQ=255)
                                              ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCcg                                                                                                        >  1:3744103/1‑150 (MQ=255)
                                                 tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCCCAAGGGCGTCACAATCCCTAGGGCGGTGATGTACCAGTAGGGGCGGGTTTAACCGCGa                                                                                                     >  2:2730547/1‑150 (MQ=255)
                                                 tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                                     >  2:7349459/1‑150 (MQ=255)
                                                 tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCAAATGCCCTAGGGCGGTGATTTACCTGTAGGAGCGGGGTTAACCGCGa                                                                                                     >  1:2075651/1‑150 (MQ=255)
                                                  ctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                    <  2:7453604/150‑1 (MQ=255)
                                                  ctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                    >  2:2828032/1‑150 (MQ=255)
                                                   ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAag                                                                                                   >  2:7644688/1‑150 (MQ=255)
                                                   ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCAAAATGGCGTGACATGCCCTAGGGCGGTGTTGTACCTGTAGGAGCGGGTTTTACCGCGaat                                                                                                   >  1:6837574/1‑149 (MQ=255)
                                                     gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGGGGTGATGTACCTGTAGGAGCGGTTTTACCCGCGAagag                                                                                                 <  2:1102258/150‑1 (MQ=255)
                                                       aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTcac                                                                                                                                                    >  2:6300446/1‑98 (MQ=255)
                                                       aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                               >  2:6076299/1‑150 (MQ=255)
                                                       aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                               >  2:6502262/1‑150 (MQ=255)
                                                       aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCGGGAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                               >  1:6695053/1‑150 (MQ=255)
                                                           caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGTAAACCGCACCCACAAGGGCGGCACAAGCCCAAGGGCGGGGGTGGACCAGGATGGGCGGGGTTTCCCCCGAAGAGGCCGGc                                                                                           >  1:6215678/1‑150 (MQ=255)
                                                           caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAAACCCCTCCCACAAGGGCGGCCCAATCCCGAGGGCGGGGATGTACCGGTGGGAGCGGGTTTAACCGCGAAGAGGCCGGc                                                                                           >  1:3575423/1‑150 (MQ=255)
                                                             aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAACCCGCTCCCACAAGGGCGGCACATGCCCTAGGGCGGTGATGTACCGGTTGGAGCGGGTTTTACCCCGCAGAGGCCGGcac                                                                                         >  1:5157332/1‑150 (MQ=255)
                                                             aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcac                                                                                         >  2:4632369/1‑150 (MQ=255)
                                                             aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGGAAACCCCGCCCCCAAGGGCGGCCCCGGCCCCAGGGCGGGGGTGTTCCGGGAGGGGCGGGGTTACACCCGAAGAGggcggccc                                                                                         >  1:121939/1‑148 (MQ=255)
                                                                  gctgcTTGGGGTGGTGCGTGGGGTTCCCGGTGTCTGTTGCCTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACACGGGCGTCACTGGCCCTTGGGGGGTGATGTACCTGTCGGTGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCa                                                                                    <  2:3575423/147‑1 (MQ=37)
                                                                  gcagcTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCGAGGGCGGTGATGTACCAGTAGGAGCGGGTTTACCCGCGAAGTGGCCGGCACAGGCa                                                                                    >  2:822957/1‑150 (MQ=255)
                                                                   cagcTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCCAAAGGGCGTCACAATGCCTAGGGCGGTGTTGTTACGGTATGAGCGGGTTTAACCGCGAAGAGGGCGGCACAGGcca                                                                                   >  1:6350033/1‑148 (MQ=255)
                                                                    agcTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAac                                                                                  >  2:4740935/1‑150 (MQ=255)
                                                                                        gTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTATACCCGGTGCCACAAGGGCGTCACATGGCCTTGGGGGGTGTTGTTGGTGTTGGTGCGGGTTTTCCCGCGAAGTGGGCGGGACAGGGAACGTCGTTAGGCGAAGCgggct                                                              >  2:2941831/1‑147 (MQ=255)
                                                                                                   cttttcaCTGTCCCACCCTCTACCCGTCTAAACCCGCACCCACAAGGCCTTCACATCCCATAGGGCGGTTATTTACATTTAGTAGCGTTTTTACCCGCCAAGAGGCCGGAACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCa                                                   <  2:1999714/144‑1 (MQ=255)
                                                                                                           ggggCCTGCGTCTTCGGGGGTAACCCGGCTCCCACCGGGGGGTCACTGGCCCTTGGGGGGTTCTTTACGTTTTGGAGGGGGTTCACCCGCTAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGctc                                           <  2:4009050/147‑1 (MQ=255)
                                                                                                                    ctctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAACCGCGAAGAGGCCGGCACAGGCAACGACGATAACCGAACCGGCCACCACGCCCTCAGGCTGCTCTTTCGTCTg                                  >  1:6072460/1‑150 (MQ=255)
                                                                                                                       ttCGGGGGCAAACCCTCTCCCACAAGGGCGTCACATGCCCTAGGGCGTTTATTTACCTGTAGTAGCGGGTTTACCCGCGAAGTGGCCGGCACATGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAt                               <  2:9135782/150‑1 (MQ=255)
                                                                                                                       ttCGCGGGTAAACCCGCGCCCCCCAGGGCGTCACATGCCCTAGGGCGGTGAGGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAt                               <  2:5156529/150‑1 (MQ=255)
                                                                                                                         cgggggAAACCCCGCCCCCACACGGGCGACACAGGCCCTTGGGCGGTTATGTACCTGTCGGAGGGGGTTAACCCGCGAAGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                             <  2:1585001/147‑1 (MQ=255)
                                                                                                                          gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                            >  1:3135993/1‑150 (MQ=255)
                                                                                                                          gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                            >  1:7222634/1‑150 (MQ=255)
                                                                                                                           gggtaaaaCCCGCCCCCCCACGGGCGCCACAGGCCCTGGGGGGGTGCTGGACCTGTCGGTGGGTTTTCCCCCGCGAAGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGt                           <  2:1720215/145‑1 (MQ=255)
                                                                                                                            gggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTc                          >  1:7209063/1‑150 (MQ=255)
                                                                                                                             gttAAACCCGCCCCCACCCGGGCGTCACTGGCCCTTGGGGGGTGTTGTACCTGTCGGAGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                         <  2:2677690/148‑1 (MQ=255)
                                                                                                                             ggaacaCCCGCCCCCACCCGGGCGTCACAGGCCCTCGGGGGTTGTTGACCGTGTAGGCGGGGTTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                         <  2:1574888/145‑1 (MQ=255)
                                                                                                                             ggaaaccccGCTCCCACACGGGCGTCACTGGCCCTCGGGCGGTGATGTCCCTGTTGGCGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                         <  2:6176391/144‑1 (MQ=38)
                                                                                                                             ggTAACCCCGCTCCCACACGGGCGTCACATGCCCTAGGGGGGTGATGTACCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                         <  2:8709801/150‑1 (MQ=255)
                                                                                                                               tAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGGGGTGTTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgc                       <  2:4524914/150‑1 (MQ=255)
                                                                                                                                aaaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcg                      >  1:7314507/1‑150 (MQ=255)
                                                                                                                                 aaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcgc                     >  1:5267095/1‑150 (MQ=255)
                                                                                                                                  aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                    >  2:5609756/1‑150 (MQ=14)
                                                                                                                                  aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                    >  2:578276/1‑150 (MQ=255)
                                                                                                                                   cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                   >  2:1154496/1‑150 (MQ=255)
                                                                                                                                   accGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAACTGTAGGAGCGGGTTTACCCGCGAAGATGCAGGCACAGGCAACGACGATACACGAACCGGCCAGCAATCCCTCAGGATGCTCTTTCGTCTGGATCTGTCAGCGCtt                   >  2:2763937/2‑150 (MQ=38)
                                                                                                                                    ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                  >  2:3263958/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                  >  2:8832685/1‑150 (MQ=255)
                                                                                                                                       cTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTgg               >  2:9028512/1‑150 (MQ=255)
                                                                                                                                        tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg              >  1:6218933/1‑150 (MQ=255)
                                                                                                                                          ccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGtt            >  1:2167764/1‑150 (MQ=255)
                                                                                                                                           cacaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAgg                                                 <  1:2877530/112‑1 (MQ=255)
                                                                                                                                            acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg              <  1:7816053/146‑1 (MQ=255)
                                                                                                                                            acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg              >  2:7816053/1‑146 (MQ=255)
                                                                                                                                            acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg          >  1:3381607/1‑150 (MQ=255)
                                                                                                                                            acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg          >  1:5665547/1‑150 (MQ=255)
                                                                                                                                            acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg          >  1:8291531/1‑150 (MQ=255)
                                                                                                                                               aGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgcgc       >  2:885040/1‑150 (MQ=255)
                                                                                                                                                gggCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCa      >  1:8194928/1‑150 (MQ=255)
                                                                                                                                                    gTCACATGCCCTAGGGCGGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGgc  >  1:1878332/1‑150 (MQ=38)
                                                                                                                                                    |                                                                                                                                                      
GGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGC  >  NC_002947/4862981‑4863278

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: