Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,496,816 T→C 7.8% intergenic (+48/+24) mdeA → / ← PP_1309 methionine gamma‑lyase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,496,8160TC 7.8% 183.6 / 3.6 64intergenic (+48/+24)mdeA/PP_1309methionine gamma‑lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base C (2/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCG  >  NC_002947/1496668‑1496961
                                                                                                                                                     |                                                                                                                                                  
aTCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTCGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCtt                                                                                                                                                    <  6:215394/149‑1 (MQ=255)
        atatCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAATAGGCCTCTTGCGAGGcc                                                                                                                                            <  8:357423/149‑1 (MQ=255)
             ggAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGt                                                                                                                                      <  5:25183/149‑1 (MQ=255)
              gAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCGCAAGAGGCCtttttt                                                                                                                                      <  1:161362/149‑4 (MQ=255)
               aGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGCATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTca                                                                                                                                    >  2:328376/1‑149 (MQ=255)
                 gggTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCtt                                                                                                                                                    <  4:188889/132‑1 (MQ=255)
                 gggTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCtt                                                                                                                                                    >  3:188889/1‑132 (MQ=255)
                  ggTTGGTGCGGCTGGCGGGGGGGCTGGAGGATGTGGAGGCCCTGCTGGCAGATGTCGACCAGGCATTGCAGGCCTGCTGTTAGGTGGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACAt                                                                                                                                 <  2:281477/149‑1 (MQ=255)
                       gTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGa                                                                                                                            <  7:19080/149‑1 (MQ=255)
                       gTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGa                                                                                                                            <  5:382093/149‑1 (MQ=255)
                           gTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCaa                                                                                                                        <  8:32203/149‑1 (MQ=255)
                               gTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTCGCGAGGCCTTTGTCACATTCAGATCAACGCt                                                                                                                    >  4:277799/1‑149 (MQ=255)
                                tCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCtt                                                                                                                   <  8:265106/149‑1 (MQ=255)
                                tCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCAATTGAATCCAGCCCACAAAAAAGGCCTCTTGATAGGCCTTTGTCACATTCAGATAAACGCtt                                                                                                                   >  1:224475/1‑149 (MQ=255)
                                  ggTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTGGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGa                                                                                                                 <  8:295736/149‑1 (MQ=255)
                                   gTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACAtt                                                                                                                                >  8:241218/1‑133 (MQ=255)
                                   gTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCCCTTGCGAGGCCTTTGTCACAtt                                                                                                                                <  7:241218/133‑1 (MQ=255)
                                    tGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGc                                                                                                               <  8:12678/149‑1 (MQ=255)
                                        gCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAg                                                                                                                >  5:223505/1‑144 (MQ=255)
                                        gCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAg                                                                                                                <  6:223505/144‑1 (MQ=255)
                                         cTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCAATGATTTTAGTTTCCCCGCGAAAAGGCAATTGAAACAAGCCAAAAAAAAAGCCATATTGCGAGGCGTTTGGAAATTACAGATAAACCATTGACCCCg                                                                                                            >  1:285023/1‑148 (MQ=255)
                                                     gaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAg                                                                                              >  3:179003/1‑149 (MQ=255)
                                                      aggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGc                                                                                             <  8:112688/149‑1 (MQ=255)
                                                          ccTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGc                                                                                         >  3:291243/1‑149 (MQ=255)
                                                           ctgctgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGcc                                                                                        <  1:328376/149‑1 (MQ=255)
                                                                 gCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcgt                                                                                                  >  2:226134/1‑133 (MQ=255)
                                                                 gCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGa                                                                                  <  4:291243/149‑1 (MQ=255)
                                                                 gCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATACAGATCAACGCTTGAGCCCGTAGCTTtcgt                                                                                                  <  1:226134/133‑1 (MQ=255)
                                                                     aTGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTa                                                                              <  5:211046/149‑1 (MQ=255)
                                                                       gTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAgg                                                                            >  8:108585/1‑149 (MQ=255)
                                                                        tCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGgc                                                                           >  7:100908/1‑149 (MQ=255)
                                                                          gAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGgcgc                                                                         >  7:311520/1‑149 (MQ=255)
                                                                                gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCtt                                                                                                                   >  4:115760/1‑101 (MQ=255)
                                                                                gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCtt                                                                                                                   <  3:115760/101‑1 (MQ=255)
                                                                                gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCgg                                                                   <  1:162109/149‑1 (MQ=255)
                                                                                  aTTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGAt                                                                                                                           >  6:165216/1‑91 (MQ=255)
                                                                                  aTTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTg                                                                 >  5:99470/1‑149 (MQ=255)
                                                                                  aTTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTg                                                                 >  5:187511/1‑149 (MQ=255)
                                                                                    tGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGt                                                                  >  4:61843/1‑146 (MQ=255)
                                                                                    tGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGt                                                                  <  3:61843/146‑1 (MQ=255)
                                                                                      cAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGa                                                                                               <  3:256870/115‑1 (MQ=255)
                                                                                      cAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGa                                                                                               >  4:256870/1‑115 (MQ=255)
                                                                                       aGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAg                                                            >  4:139217/1‑149 (MQ=255)
                                                                                            ctgctgTTAGGTTGCCCGCGAAGAGGCCATTCAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgc                                                        >  3:95226/1‑149 (MQ=255)
                                                                                               ctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcgcg                                                     <  6:99470/149‑1 (MQ=255)
                                                                                                  ttAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAg                                                                                                                                               >  4:24235/1‑57 (MQ=255)
                                                                                                    aGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCAc                                                 <  6:187511/148‑1 (MQ=255)
                                                                                                                aaGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCt                                    <  1:324804/149‑1 (MQ=255)
                                                                                                                 agagGCCATTGAATCCAGCCCACAAAAAAGGCCTCTCGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCtt                                   <  3:277799/149‑1 (MQ=255)
                                                                                                                     gCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAgg                               >  6:248418/1‑149 (MQ=255)
                                                                                                                       cATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGGTCACATTAAGATCAACGCTTGAGCCCGTAGCTTtcgtcg                                                                                                >  1:73450/1‑83 (MQ=255)
                                                                                                                       cATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcgtcg                                                                                                <  2:73450/82‑1 (MQ=255)
                                                                                                                         ttGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcgc                                                      <  4:93650/122‑1 (MQ=255)
                                                                                                                         ttGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATAAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcgc                                                      >  3:93650/1‑122 (MQ=255)
                                                                                                                           gAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGt                         >  4:296424/1‑149 (MQ=255)
                                                                                                                               ccAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAgg                     <  8:311520/149‑1 (MQ=255)
                                                                                                                                 aGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCg                   <  4:179003/149‑1 (MQ=255)
                                                                                                                                    ccACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGt                 >  2:97686/1‑148 (MQ=255)
                                                                                                                                      acaAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAc              >  8:110112/1‑149 (MQ=255)
                                                                                                                                        gacaaaGGCCTCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcgc            >  4:182344/4‑149 (MQ=255)
                                                                                                                                        gacaaaGGCCTCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcgc            <  4:325393/146‑1 (MQ=255)
                                                                                                                                        aaaaaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAc              >  8:173910/1‑147 (MQ=255)
                                                                                                                                          aaaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcggcg                                          >  4:181802/1‑117 (MQ=255)
                                                                                                                                          aaaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcggcg                                          <  3:181802/117‑1 (MQ=255)
                                                                                                                                          aaaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCgg          <  2:224475/149‑1 (MQ=255)
                                                                                                                                                 ctctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGc   >  5:72556/1‑149 (MQ=255)
                                                                                                                                                  tctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCg  >  3:134259/1‑149 (MQ=255)
                                                                                                                                                    ttGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCg  <  4:134259/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCG  >  NC_002947/1496668‑1496961

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: