Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,458,946 Δ1 bp 6.7% intergenic (+23/‑39) leuS → / → PP_4795 leucine‑‑tRNA ligase/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,458,9460T. 6.7% 268.7 / 13.8 89intergenic (+23/‑39)leuS/PP_4795leucine‑‑tRNA ligase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (42/41);  new base . (2/4);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

GCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGA  >  NC_002947/5458798‑5459085
                                                                                                                                                     |                                                                                                                                            
gCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGAtt                                                                                                                                             <  1:2047505/150‑1 (MQ=255)
  cAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGc                                                                                                                                           >  1:1057855/1‑150 (MQ=255)
    gCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCgg                                                                                                                                        >  2:1139795/1‑150 (MQ=255)
      cgcgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATCGCCTGACCATCCGCAAGGTCACCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATTAAAACGCCCACCCGCATCGATTGCGGGt                                                                                                                                      <  2:94881/150‑1 (MQ=255)
       gcgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTg                                                                                                                                     <  1:2087331/150‑1 (MQ=255)
       gcgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTg                                                                                                                                     <  2:936723/150‑1 (MQ=255)
        cgacgaAGTCGAAGGCGCGGACAGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGATCTTCCGCAAGTTTATCGTGGTACCGGCCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCTTCGATTGCGGGTgg                                                                                                                                    <  2:786094/150‑1 (MQ=255)
        cgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                    <  1:86491/150‑1 (MQ=255)
           agaagTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCg                                                                                                                                 <  1:390291/150‑1 (MQ=255)
               gTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAAc                                                                                                                             >  1:1462236/1‑150 (MQ=255)
                  gAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAgg                                                                                                                          >  1:352265/1‑150 (MQ=255)
                  gAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAgg                                                                                                                          <  1:824356/150‑1 (MQ=255)
                      ccGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGccc                                                                                                                       >  2:1429613/1‑149 (MQ=255)
                      ccGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGccc                                                                                                                       <  1:1429613/149‑1 (MQ=255)
                      ccGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCa                                                                                                                      >  1:2084980/1‑150 (MQ=255)
                          tGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAg                                                                                                                  <  1:2189017/150‑1 (MQ=255)
                           gCCCGCAGCAACGAGAATGTCATGAGCTACATCTATGGCCTGACCATCCGCAATGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAgg                                                                                                                 <  2:1775943/150‑1 (MQ=255)
                            cccgttgcaACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAggg                                                                                                                >  2:1203776/7‑150 (MQ=255)
                            cccGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAggg                                                                                                                >  1:2192815/1‑150 (MQ=255)
                             ccGCAGCAACGAGAATGTCCTGCGCTTCATTGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGTGCAAGCTGGTCAACATCGTAGCCAACTGATGACAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                               <  1:2293486/150‑1 (MQ=255)
                                  gcaACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaac                                                                                                          <  2:606931/150‑1 (MQ=255)
                                    aaCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCg                                                                                                                                 <  2:1516387/125‑1 (MQ=255)
                                    aaCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCg                                                                                                                                 >  1:1516387/1‑125 (MQ=255)
                                     aCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCaacaac                                                                                                       >  1:616057/1‑150 (MQ=255)
                                       gagaATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCTGCGCCAACTGATGAACACGCCCCCCCGCAGCTATCGCGGTTCGGCGGAACAGGCCCAAAAGGGAGCAACAACAt                                                                                                     >  1:2271466/1‑150 (MQ=255)
                                       gagaATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGGCGCCAACTGATGACCACGATACTCTGCGGGATGACGTGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACAt                                                                                                     >  1:2271834/1‑150 (MQ=255)
                                         gaATGTCCTGCGCTTCATCGATGGCCTGACGATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTACGGGTGGGAAGAACAGGCCCCAAAGGGAGAAACAAGAtcg                                                                                                   >  2:403291/1‑148 (MQ=255)
                                             gTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCaa                                                                                               <  1:596103/150‑1 (MQ=255)
                                              tCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCaaa                                                                                              >  1:1662432/1‑150 (MQ=255)
                                              tCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATAAAAACGCCCAACCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCaaa                                                                                              <  1:939584/150‑1 (MQ=255)
                                                cTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACg                                                                                            <  2:616057/150‑1 (MQ=255)
                                                 tGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCaaa                                                                                              >  2:673388/1‑147 (MQ=255)
                                                 tGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCaaa                                                                                              <  1:673388/147‑1 (MQ=255)
                                                   cgcTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCaa                                                                                         <  1:2050259/150‑1 (MQ=255)
                                                      ttCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAAttt                                                                                      >  1:1003073/1‑150 (MQ=255)
                                                      ttCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAAttt                                                                                      >  1:2258468/1‑150 (MQ=255)
                                                       tCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTg                                                                                     <  1:1100662/150‑1 (MQ=255)
                                                          tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATAGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                                  >  2:207096/1‑150 (MQ=255)
                                                            gATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGt                                                                                >  2:1247941/1‑150 (MQ=255)
                                                                     aCCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAg                                                                                                              >  2:90751/1‑111 (MQ=255)
                                                                     aCCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAg                                                                                                              <  1:90751/111‑1 (MQ=255)
                                                                      ccATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTg                                                                      >  2:1953356/1‑150 (MQ=255)
                                                                        aTCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGc                                                                    >  1:1987082/1‑150 (MQ=255)
                                                                             cAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCa                                                               <  2:1987082/150‑1 (MQ=255)
                                                                               aGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATg                                                             <  1:207096/150‑1 (MQ=255)
                                                                               aGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATg                                                             <  2:442377/150‑1 (MQ=255)
                                                                                   cATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGCACTGGCGGTCATGCTCa                                                         <  1:1247941/150‑1 (MQ=255)
                                                                                     tCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAgc                                                       <  2:1003073/150‑1 (MQ=255)
                                                                                       gTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACAGGATCAAACGCAATTTGCTGGTAATCGGGCTGGCGGTCATGCTCAgcgc                                                     >  1:596069/1‑150 (MQ=255)
                                                                                         ggTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTgg                                                                                                                                    <  2:171225/69‑1 (MQ=255)
                                                                                         ggTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTgg                                                                                                                                    >  1:171225/1‑69 (MQ=255)
                                                                                          gTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGAtt                                                                                                                                             <  2:1533964/60‑1 (MQ=255)
                                                                                          gTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGAtt                                                                                                                                             >  1:1533964/1‑60 (MQ=255)
                                                                                              cGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGt                                              <  2:2258468/150‑1 (MQ=255)
                                                                                               gggCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGtt                                             >  1:945623/1‑150 (MQ=255)
                                                                                                 gCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTc                                           <  2:2192815/150‑1 (MQ=255)
                                                                                                   aaGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCa                                         <  2:352265/150‑1 (MQ=255)
                                                                                                     gCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAACAGGGTTCAACAACATGATCAAACGCCCTTTGCTACGCATCGGCCTTTCGATCATGCTCAGCGCCTGCGGTTTCCAGc                                       >  1:2270332/1‑150 (MQ=255)
                                                                                                     gCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGc                                       >  2:1984856/1‑150 (MQ=255)
                                                                                                       tGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTg                                     >  1:2245693/1‑150 (MQ=255)
                                                                                                         gTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCt                                                                                   <  1:76132/102‑1 (MQ=255)
                                                                                                         gTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCt                                                                                   >  2:76132/1‑102 (MQ=255)
                                                                                                         gTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTgcg                                   >  1:789557/1‑150 (MQ=255)
                                                                                                            aaCATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGTCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCgg                                >  2:774131/1‑150 (MQ=255)
                                                                                                            aaCATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAgc                                                       <  2:1379685/127‑1 (MQ=255)
                                                                                                            aaCATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAgc                                                       >  1:1379685/1‑127 (MQ=255)
                                                                                                              cATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCa                              >  2:811320/1‑150 (MQ=255)
                                                                                                               aTCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGctg                             <  2:1057855/150‑3 (MQ=255)
                                                                                                                tcgtcgCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCTCGGCAcc                            >  2:795334/1‑150 (MQ=255)
                                                                                                                tcgtcgCCAACTGATGAAAACGCCCACCCGCATCGATTGANGGTGGGCGTAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCAcc                            >  1:1343294/1‑150 (MQ=255)
                                                                                                                 cgcctccAACTGCCGAAAGCCCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCg                           <  2:1662432/145‑1 (MQ=255)
                                                                                                                  gtcgCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCgg                          >  2:1923257/1‑150 (MQ=255)
                                                                                                                   tcgCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGTAACAGGCCCAAAAGTGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCTGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGgc                         <  1:1139795/150‑1 (MQ=255)
                                                                                                                     gCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGgcgc                       <  2:789557/150‑1 (MQ=255)
                                                                                                                     gCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGgcgc                       >  1:273296/1‑150 (MQ=255)
                                                                                                                      ccAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGc                                                                         >  1:2337882/1‑99 (MQ=255)
                                                                                                                      ccAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGc                                                                         <  2:2337882/99‑1 (MQ=255)
                                                                                                                      ccAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGAGCCTGCGGTTTCCAGCTGAGAGGCACAGGCGcc                      >  1:1731662/1‑150 (MQ=255)
                                                                                                                        aaCTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCAc                    >  1:421847/1‑150 (MQ=255)
                                                                                                                          cTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGTCACCg                  >  1:1725550/1‑150 (MQ=255)
                                                                                                                           tgatgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGa                 >  1:1654561/1‑150 (MQ=255)
                                                                                                                             atgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGTCACCGAGc               <  2:1725550/150‑1 (MQ=255)
                                                                                                                             atgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGc               <  2:945623/150‑1 (MQ=255)
                                                                                                                              tgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCt              <  1:1953356/150‑1 (MQ=255)
                                                                                                                                  aaCGCCCACCCTCACCGATTGCGGGTGGGCGGAACAGGCCCAAAAGTGAGCAACAACATGATCAAACGCAATTTGCTGTTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGTGGCACCGGCGCCACCGAGCTGAGc          <  2:15158/150‑1 (MQ=255)
                                                                                                                                     gCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTg       <  1:811320/150‑1 (MQ=255)
                                                                                                                                      cccACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGa      <  2:2084980/150‑1 (MQ=255)
                                                                                                                                       ccACCCGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACAt                                                                                                     >  2:2316451/1‑54 (MQ=21)
                                                                                                                                       ccACCCGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACAt                                                                                                     <  1:2316451/54‑1 (MQ=21)
                                                                                                                                       ccACCCGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCg                                                                   >  1:1601010/1‑88 (MQ=38)
                                                                                                                                       ccACCCGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCg                                                                   <  2:1601010/88‑1 (MQ=38)
                                                                                                                                          cccGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGa  <  2:2245693/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGA  >  NC_002947/5458798‑5459085

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