Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,863,127 G→A 7.2% intergenic (‑111/+126) 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,1270GA 7.2% 742.4 / 30.5 238intergenic (‑111/+126)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base G (126/94);  new base A (2/15);  total (128/110)
Fisher's exact test for biased strand distribution p-value = 2.63e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 240 total aligned reads displayed.

TCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCA  >  NC_002947/4862979‑4863274
                                                                                                                                                    |                                                                                                                                                    
tCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGc                                                                                                                                                     <  1:2282303/150‑1 (MQ=255)
tCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGc                                                                                                                                                     <  1:3383763/150‑1 (MQ=255)
  ggCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGt                                                                                                                                                   <  1:7688935/150‑1 (MQ=255)
   gCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTc                                                                                                                                                  <  2:2521317/150‑1 (MQ=255)
             cGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTa                                                                                                                                        <  2:7180115/150‑1 (MQ=255)
              gAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTCCCACCAGGGCGTCACATGCCCTAg                                                                                                                                       >  1:2674931/1‑150 (MQ=255)
               aGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAgg                                                                                                                                      <  1:655523/150‑1 (MQ=255)
               aGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAgg                                                                                                                                      <  1:5638568/150‑1 (MQ=255)
                  aTACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCCGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCg                                                                                                                                  <  1:3170280/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                                                                                                                                <  1:420719/148‑1 (MQ=255)
                      gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGt                                                                                                                                >  2:420719/1‑148 (MQ=255)
                      gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                              <  2:6452845/150‑1 (MQ=255)
                       ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                             <  2:7042405/150‑1 (MQ=255)
                       ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                             <  2:538603/150‑1 (MQ=255)
                       ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                             <  2:4985320/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)
                                 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)
                                                    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:6502262/1‑150 (MQ=255)
                                                         aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                           >  2:6076299/1‑150 (MQ=255)
                                                         aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCGGGAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                           >  1:6695053/1‑150 (MQ=255)
                                                         aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGGCACAAGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTAACCGCGAAGAGGc                                                                                           >  1:2052401/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)
                                                                             ggtggtGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTGAACCCCCGCCCACAAGGGCGGCACATGCCCAAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCCCAGGCAACCACGATa                                                                       >  1:2798588/1‑150 (MQ=38)
                                                                                          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)
                                                                                                                           ggcgGTTAACCCCGCCCCCACTGGGGCTCCACAGGCCCTTGGGGGGTTTTTTCCCTGTTGGTGGGGTTTAACCGGCTAAGCGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                         <  1:2937173/149‑1 (MQ=255)
                                                                                                                           cgggggAAACCCCGCCCCCACACGGGCGACACAGGCCCTTGGGCGGTTATGTACCTGTCGGAGGGGGTTAACCCGCGAAGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                         <  2:1585001/147‑1 (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)
                                                                                                                                              acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg          >  2:7816053/1‑146 (MQ=255)
                                                                                                                                              acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg          <  1:7816053/146‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
TCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCA  >  NC_002947/4862979‑4863274

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: