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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,863,1200CT 7.4% 707.5 / 31.5 228intergenic (‑104/+133)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base C (117/94);  new base T (2/15);  total (119/109)
Fisher's exact test for biased strand distribution p-value = 5.91e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

Only 100 of 240 total aligned reads displayed.

ATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTA  >  NC_002947/4862972‑4863269
                                                                                                                                                     |                                                                                                                                                      
    aGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAg                                                                                                                                                   <  2:5551380/150‑1 (MQ=255)
     gCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCCGTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCtacaggt                                                                                                                                                  <  1:5070108/150‑4 (MQ=38)
       tCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGc                                                                                                                                                <  1:3383763/150‑1 (MQ=255)
          gCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTc                                                                                                                                             <  2:2521317/150‑1 (MQ=255)
           cACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTca                                                                                                                                            >  2:7357635/1‑150 (MQ=255)
                    cGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTa                                                                                                                                   <  2:7180101/150‑1 (MQ=255)
                     gAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTCCCACCAGGGCGTCACATGCCCTAg                                                                                                                                  >  1:2674931/1‑150 (MQ=255)
                      aGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAgg                                                                                                                                 <  1:5638568/150‑1 (MQ=255)
                         aTACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCCGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCg                                                                                                                             <  1:3170280/150‑1 (MQ=255)
                           ccGGGTTGCGGCGGTCATGGGAGTTGGCGTCTCTTGCCAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGGTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCggg                                                                                                                           <  1:5413859/149‑2 (MQ=39)
                           aCGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGt                                                                                                                           <  1:1620523/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:538603/150‑1 (MQ=255)
                              ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                        <  2:7042405/150‑1 (MQ=255)
                              ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                        <  2:411699/150‑1 (MQ=255)
                              ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                        <  1:2254990/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)
                                      cggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGt                                                                                                                <  1:6210159/150‑1 (MQ=255)
                                        gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                                              <  2:431514/150‑1 (MQ=255)
                                           aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGgaggag                                                                                                           >  2:8855564/1‑150 (MQ=255)
                                            tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                          >  2:7156682/1‑150 (MQ=255)
                                            tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCGCCCACAAAGGCGACACATGCCCAAGGGCGGTGATGTACCGGGAGGAGc                                                                                                          >  2:3041973/1‑150 (MQ=255)
                                            tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAAACCCACCCACAAAGGCGTCACAATCCCAAGGGCGGGGATGTAACAGTATGAGc                                                                                                          >  1:4649316/1‑150 (MQ=255)
                                              ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAACCCCCGCCCACAAAGGCGGCACAAGCCCAAGGGCGGGGTTGTAACAGTAGGAGCgg                                                                                                        >  1:8420756/1‑150 (MQ=255)
                                              ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCACCACGGCGTGACAAGCCCTACGGCGGGGATGTAACAGTATGAGCgg                                                                                                        >  1:8421440/1‑150 (MQ=255)
                                               gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCggg                                                                                                       <  2:4765067/150‑1 (MQ=14)
                                               gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                       <  2:1857871/150‑1 (MQ=255)
                                               gATTTGTATTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCggg                                                                                                       <  1:7592946/150‑1 (MQ=14)
                                                  ttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGttt                                                                                                    <  1:2158915/150‑1 (MQ=255)
                                                   tGTCTTCTCTTGCAATGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGGCCTATGGCGGTGTTGTACCTGTAGGAGCGGGTTTa                                                                                                   >  2:1013542/1‑150 (MQ=255)
                                                   tGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTa                                                                                                   <  1:3529670/150‑1 (MQ=255)
                                                    gTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAc                                                                                                  <  2:1049417/150‑1 (MQ=14)
                                                     tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAcc                                                                                                 <  1:388348/150‑1 (MQ=14)
                                                     tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAcc                                                                                                 >  2:7743231/1‑150 (MQ=255)
                                                     tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAcc                                                                                                 >  1:8272098/1‑150 (MQ=255)
                                                      cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTATCTGTAGGAGCGGGTTTAccc                                                                                                >  2:2728755/1‑150 (MQ=255)
                                                      cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTCCCACAACGGCGTCACAATCCCAAGGGCGGTGATGTAACAGTAGGAGCGGGTTTAacc                                                                                                >  1:7037423/1‑150 (MQ=255)
                                                       ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCcg                                                                                               <  1:9027265/150‑1 (MQ=14)
                                                       ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcg                                                                                               >  2:8241344/1‑150 (MQ=255)
                                                       ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCcg                                                                                               >  1:3744103/1‑150 (MQ=255)
                                                        tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATTCCCTAGGGCGGGGGTGTACCGGTAGGTGCGGGGTTTCCcgc                                                                                              >  1:6572425/1‑150 (MQ=38)
                                                        tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcgc                                                                                              <  2:3378502/150‑1 (MQ=255)
                                                         ctctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAAACAGCGCCCACAAGGGCGTGACATAACCAAGGGCGGGGATGTAAAGGGAGGAGAGGGttttcacccg                                                                                             >  1:7789607/1‑143 (MQ=255)
                                                          tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCCCAAGGGCGTCACAATCCCTAGGGCGGTGATGTACCAGTAGGGGCGGGTTTAACCGCGa                                                                                            >  2:2730547/1‑150 (MQ=255)
                                                          tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                            >  2:6566248/1‑150 (MQ=255)
                                                          tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                            >  2:7349459/1‑150 (MQ=255)
                                                          tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAATGGGGTCACATGCCCTAGGGCGGTGATGTTGCTGTTGGTGGGGGTTTTCCCGCGa                                                                                            >  1:8903979/1‑150 (MQ=255)
                                                           ctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGaa                                                                                           >  2:2828032/1‑150 (MQ=255)
                                                            ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAag                                                                                          >  2:7644688/1‑150 (MQ=255)
                                                            ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCAAAATGGCGTGACATGCCCTAGGGCGGTGTTGTACCTGTAGGAGCGGGTTTTACCGCGaat                                                                                          >  1:6837574/1‑149 (MQ=255)
                                                              gcaaagcaaaATAGCCGCTTGGGGTGGTGCGGGGGGTTCCCGGGGTCTGTTGACTGTGCCAGCCTCTGCGCGGGTAACCCCGCCCCCACACGGGCGTCACATGCCCTCGGGGGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAagag                                                                                        <  2:1102353/150‑1 (MQ=255)
                                                              gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAATCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCTGTAGGAGCGGGTTTACCCGCGAagag                                                                                        >  2:6804933/1‑150 (MQ=255)
                                                                aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTcac                                                                                                                                           <  1:6300446/98‑1 (MQ=255)
                                                                aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                      >  2:6076299/1‑150 (MQ=255)
                                                                aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGGCACAAGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTAACCGCGAAGAGGc                                                                                      >  1:2052401/1‑150 (MQ=255)
                                                                aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATTCCCTAGGGCGGTGTTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                      >  1:5945873/1‑150 (MQ=255)
                                                                    caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAAACCCCTCCCACAAGGGCGGCCCAATCCCGAGGGCGGGGATGTACCGGTGGGAGCGGGTTTAACCGCGAAGAGGCCGGc                                                                                  >  1:3575423/1‑150 (MQ=255)
                                                                    caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAAACCCCGACCACAACGGCGTCACAAGCCCGAGGGGGGGGGGGTAACGGTGGGGGCGGGGTTAACCGCGAAGcggccggc                                                                                  >  1:6216209/1‑150 (MQ=255)
                                                                     aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGca                                                                                 >  2:4956688/1‑150 (MQ=255)
                                                                      aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcac                                                                                >  2:4632369/1‑150 (MQ=255)
                                                                      aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGGAAACCCCGCCCCCAAGGGCGGCCCCGGCCCCAGGGCGGGGGTGTTCCGGGAGGGGCGGGGTTACACCCGAAGAGggcggccc                                                                                >  1:121939/1‑148 (MQ=255)
                                                                      aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAACCCCCTCCCCAAAGGGCGGCACATGCCCAAGGGCGGGGATTTTCCTGTAGGAGGGGGGGTACCCGCCAAGACGCCGGcgc                                                                                >  1:5157373/1‑148 (MQ=255)
                                                                           gctgcTTGGGGTGGTGCGTGGGGTTCCCGGTGTCTGTTGCCTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACACGGGCGTCACTGGCCCTTGGGGGGTGATGTACCTGTCGGTGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCa                                                                           <  2:3575423/147‑1 (MQ=37)
                                                                                    ggtggtGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTGAACCCCCGCCCACAAGGGCGGCACATGCCCAAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCCCAGGCAACCACGATa                                                                  >  1:2798588/1‑150 (MQ=38)
                                                                                                 gTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGCTACCCGAAACGGCCa                                                     >  2:3216792/1‑150 (MQ=255)
                                                                                                     cAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGGCACATGCCCTAGGGCGGTGATGTAACGGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCAc                                                 >  2:502774/1‑150 (MQ=255)
                                                                                                      aGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGAATCCCGAACCGGCCCGCACg                                                >  2:1329628/1‑150 (MQ=255)
                                                                                                            cttttcaCTGTCCCACCCTCTACCCGTCTAAACCCGCACCCACAAGGCCTTCACATCCCATAGGGCGGTTATTTACATTTAGTAGCGTTTTTACCCGCCAAGAGGCCGGAACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCa                                          <  2:1999714/144‑1 (MQ=255)
                                                                                                               ttGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGgc                                       >  2:4665463/1‑150 (MQ=255)
                                                                                                                 gACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGGGATGTACCGGTAGGAGCGGGGTTACCCGCGGAGAGGCCGGCACAGGCCACGACGATACCCCAACCGGCCAGCCCGCCCACAGgcgg                                     >  2:4869954/1‑148 (MQ=38)
                                                                                                                            ccTCTCCGGGGTTAAACCCGCCCCCACACGGGCGACACTTGCCCTCGGGGGGTGATGTCCCTTTAGGAGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCt                          <  2:2776520/150‑1 (MQ=255)
                                                                                                                               ctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                       >  2:5663207/1‑70 (MQ=255)
                                                                                                                                ttCGGGGGCAAACCCTCTCCCACAAGGGCGTCACATGCCCTAGGGCGTTTATTTACCTGTAGTAGCGGGTTTACCCGCGAAGTGGCCGGCACATGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAt                      <  2:9135782/150‑1 (MQ=255)
                                                                                                                                  ggcgGTTAACCCCGCCCCCACTGGGGCTCCACAGGCCCTTGGGGGGTTTTTTCCCTGTTGGTGGGGTTTAACCGGCTAAGCGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                    <  1:2937173/149‑1 (MQ=255)
                                                                                                                                   ggggggAAACCCGGCCCCCACCCGGGCTTCACATGCCTTTGGGGGTTTATGCACGTGTCGGTGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                    <  2:3954385/148‑1 (MQ=255)
                                                                                                                                   gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                   >  1:7222634/1‑150 (MQ=255)
                                                                                                                                    cgGGTAAACCCGCTCCCACATGGGCGTCACTTGCCCTTGGGCGGTGTTGTACCTGTTGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGt                  <  2:5047520/150‑1 (MQ=38)
                                                                                                                                     gggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTc                 >  1:7209063/1‑150 (MQ=255)
                                                                                                                                     cggTAAACCCGGCCCCCCCACGGGCTTCAATGGCCCTGGGGGGGTTTTGTACGTGTTGGCGGGGTTTTACCCGCTAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGt                  <  2:1719709/149‑1 (MQ=255)
                                                                                                                                      ggaacaCCCGCCCCCACCCGGGCGTCACAGGCCCTCGGGGGTTGTTGACCGTGTAGGCGGGGTTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                <  2:1574888/145‑1 (MQ=255)
                                                                                                                                      ggaaaccccGCTCCCACACGGGCGTCACTGGCCCTCGGGCGGTGATGTCCCTGTTGGCGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                <  2:6176391/144‑1 (MQ=38)
                                                                                                                                      ggTAACCCCGCCCCCACCCGGGGGTCACTTGCCCTTGGGGGGTTTTGTACGTGTTGGTGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                <  2:6193788/150‑1 (MQ=255)
                                                                                                                                        tAAACCCGCTCCCCCACGGGCGTCACAGGCCCTTGGGCGGTGATGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgc              <  2:5157880/150‑1 (MQ=39)
                                                                                                                                        tAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGGGGTGTTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgc              <  2:4524914/150‑1 (MQ=255)
                                                                                                                                        tAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgc              >  1:6460338/1‑150 (MQ=255)
                                                                                                                                         aaaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcg             >  1:7314507/1‑150 (MQ=255)
                                                                                                                                           aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt           >  2:23550/1‑150 (MQ=255)
                                                                                                                                           aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt           >  2:578276/1‑150 (MQ=255)
                                                                                                                                           aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt           >  2:8714532/1‑150 (MQ=255)
                                                                                                                                             ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc         >  2:3263958/1‑150 (MQ=255)
                                                                                                                                             ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc         >  2:6489889/1‑150 (MQ=255)
                                                                                                                                             ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc         >  2:8832685/1‑150 (MQ=255)
                                                                                                                                                 cTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTgg      >  2:9028512/1‑150 (MQ=255)
                                                                                                                                                     cacaCGGGCGTCACATGCCCTTGGGGGGTTTTGTACCTTTTGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTa  <  2:2719906/150‑1 (MQ=38)
                                                                                                                                                     |                                                                                                                                                      
ATCGAGCTCGGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTA  >  NC_002947/4862972‑4863269

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: