Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,716,030 T→G 10.2% intergenic (+20/+36) eamB → / ← grcA cysteine and O‑acetylserine exporter/autonomous glycyl radical cofactor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,716,0300TG10.2% 221.1 / 21.0 79intergenic (+20/+36)eamB/grcAcysteine and O‑acetylserine exporter/autonomous glycyl radical cofactor
Reads supporting (aligned to +/- strand):  ref base T (36/35);  new base G (3/5);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-03

TTTGCAGTTTGTGAACGTCAAAATCATTTTGTACGGTGTTACGGCACTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGGTGTTGACCACCTTCAACGCGAACTTCTGGTTTCACTTCTACT  >  NC_000913/2715758‑2716276
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                           
tttgcATTTTGTGAACGTCAAAATCATTTTGTACGGTGTTACGGCAGTTTGGACGTTTGTTCGGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGGGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCg                                                                                                                                                                                                                                                           <  2:6183/275‑1 (MQ=255)
    cAGTTTGTGAACGTCAAAATCATTTTGTACGGTGTTACGGCACTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGACAGTATGGTCGCAAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGcccg                                                                                                                                                                                                                                                       >  2:453639/1‑274 (MQ=255)
                      aTCATTTTGTACGTGGTTACGGCACTGTGGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGGGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTg                                                                                                                                                                                                                                     <  2:119250/275‑1 (MQ=255)
                                         cGGCACTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGTTCGCCCTCTTCCGCTTAGTCACTTGC‑TACTTAATCCTTACAgg                                                                                                                                                                                                                  >  2:230100/1‑274 (MQ=255)
                                          ggCACTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGGGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑t                                                                                                                                                                                                               <  2:339594/275‑1 (MQ=255)
                                            cccTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAAGTGGGGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑ttt                                                                                                                                                                                                             <  2:50383/273‑1 (MQ=255)
                                              cTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCGGGGTAGTTGGCGCCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTTCTTTCTTGCCGGTCTACTTCGCGGTACGCATCTTCTATTAACGAAACAACGCGGAAGAGGTCGCTGTCTGCCGTGTAGTCA‑TTGCGTACCTAAGCCTTACAGAC‑tttta                                                                                                                                                                                                           >  4:437534/1‑273 (MQ=255)
                                                  cGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCctc                                                                                                                                                                                                                                                      <  2:284993/230‑1 (MQ=255)
                                                  cGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGGGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTa                                                                                                                                                                                                                      <  2:437539/261‑1 (MQ=255)
                                                      gTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGAGTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCGTTTTCTATTAAGGAAAAAAAGCGGAAGAGGTCGGCCTCTTCTGCTTAGTAACTGGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAgt                                                                                                                                                                                                    >  1:16010/1‑273 (MQ=255)
                                                          agacagCCGCAAACACAGGCGTTAAGCGGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACg                                                                                                                                                                                                <  2:471995/269‑1 (MQ=255)
                                                           ttCTGCCGCAAACACACGCGTTAAGCTGGTTAGTTGGCGTTAGCGCTTCGCTGGCGATGATTGGGACCTTTGGCAATGTGTGCTGGGCGCTCGCGTGGCCTCTGTTTGCGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTTGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAATGCGCAATAGGTCGCCCTCTTCCGCATAGTAACTTGC‑TACTTGAGCACTA‑AGGCATTTCAGTaaggcaaaac                                                                                                                                                                                               >  4:262816/1‑266 (MQ=255)
                                                                                   ggggggTAGTTGGCGGTAGCGTTTGGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTGCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATGGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCg                                                                                                                                                                      <  4:72265/272‑1 (MQ=255)
                                                                                      gggTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGNGTTTNGGCAATGTGTGCTGGGCGCTGGCNGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGcc                                                                                                                                                                                                                                                         <  4:505751/191‑1 (MQ=255)
                                                                                      gggTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGACAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTCAAGTTACGACCGATAACGTCGGGCTGCTGTTCCTGa                                                                                                                                                                    >  1:40926/1‑274 (MQ=255)
                                                                                      gggTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTTCTGTTCCGGa                                                                                                                                                                    >  2:122363/1‑274 (MQ=255)
                                                                                       ggTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTAGCCCTCTTCCGCGTAGTAAATTGC‑TACTTAAGCCTTaaa                                                                                                                                                                                                                    >  1:85197/1‑224 (MQ=255)
                                                                                          aGTTGGCTTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTTGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCa                                                                                                                                                                <  1:99698/274‑1 (MQ=255)
                                                                                                       gTTTGGGTGGCGATTATTGGTACGTTTGGCAAGTTGTGATGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTATATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCg                                                                                                                                                   <  1:136061/274‑1 (MQ=255)
                                                                                                          ttGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCGGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGCTAACGTCGCGCTGCTGTTCCGGAGTTGGATGATTAAAGCAAAca                                                                                                                                               >  1:114001/1‑274 (MQ=255)
                                                                                                                               tttGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGCCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGAGGAAGAGGTCGCACACTTCCGCTTGGTAACTGGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAGACGAAACGCAGCACCGGAAACACGGGTg                                                                                                                          >  4:131736/1‑275 (MQ=255)
                                                                                                                                    cAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCctct                                                                                                                                                                                                                                                     >  8:23589/1‑149 (MQ=255)
                                                                                                                                      aTGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGGCGACCTCTTc                                                                                                                                                                                                                                                   <  6:296722/149‑1 (MQ=255)
                                                                                                                                       tgtgtgCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTcc                                                                                                                                                                                                                                                  <  7:23589/149‑1 (MQ=255)
                                                                                                                                       tgtgtgCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGGCGACCTCTTcc                                                                                                                                                                                                                                                  <  8:378263/149‑1 (MQ=255)
                                                                                                                                       tgtgtgCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTACGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGGCGACCTCTTcc                                                                                                                                                                                                                                                  <  8:195632/149‑1 (MQ=255)
                                                                                                                                            gCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTa                                                                                                                                                                                                                                             <  5:56637/149‑1 (MQ=255)
                                                                                                                                                ggCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGt                                                                                                                                                                                                                                                             >  4:438166/1‑129 (MQ=255)
                                                                                                                                                  cgcTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAACGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAAAGTTAAAGCGAACTGCATAACCGGATCCAGGGACGGTCA‑‑‑GGTGAGGATGtttt                                                                                                     >  2:491998/1‑274 (MQ=255)
                                                                                                                                                    cgggCGGGGCAGCTGTTTCAGCGATTGTTTCGCCAGTAGGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCGGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTcc                                                                                                  <  1:53468/273‑1 (MQ=255)
                                                                                                                                                    cTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGcc                                                                                                                                                                                                                                                         <  2:438171/129‑1 (MQ=255)
                                                                                                                                                     tGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCTGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATCACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCg                                                                                                 >  2:275366/1‑275 (MQ=255)
                                                                                                                                                       ggggtgaattggTTTCAGGGTTTGTGTNCNAAGTATGGTCGCCAGTTAAATAGCGGTATTGCCCTATTGCTGATTNATTGCTCGGTATGCAATTTCTATTAACGAAAAAAAGCGGAAGAGGTAGGACTATGCCGCTGAGTTACTGNC‑TACCTAAGCCGTGCAGGC‑TTTCAGTAAAGGTACGAGCGATAACGGCGCGCTGCTGTTCCGGAGTCAGAGAGTGAAAGCGAACTGCAGAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCgg                                                                                                <  3:40926/263‑1 (MQ=255)
                                                                                                                                                               tCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGcc                                                                                                                                                                                                                         >  8:468132/1‑149 (MQ=255)
                                                                                                                                                                 tGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCtt                                                                                                                                                                                                                       >  6:55106/1‑149 (MQ=255)
                                                                                                                                                                    ttCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGACGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACa                                                                                                                                                                                                                    <  7:124624/149‑1 (MQ=255)
                                                                                                                                                                          gATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑ttt                                                                                                                                                                                                             <  8:351604/149‑1 (MQ=255)
                                                                                                                                                                             tGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGCGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATAGTTTTCCGTATGCTTAACTGCATCTTACagag                                                                         >  1:72031/1‑275 (MQ=255)
                                                                                                                                                                                 tCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTaaa                                                                                                                                                                                                      <  7:468132/149‑1 (MQ=255)
                                                                                                                                                                                          tGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGCa                                                             <  2:265600/274‑1 (MQ=255)
                                                                                                                                                                                           ggTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCg                                                                                                                                                                                            <  5:55106/149‑1 (MQ=255)
                                                                                                                                                                                                ccAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGa                                                                                                                                                                                           <  7:172817/145‑1 (MQ=255)
                                                                                                                                                                                                ccAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGa                                                                                                                                                                                           >  8:172817/1‑145 (MQ=255)
                                                                                                                                                                                                    ttAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTcg                                                                                                                                                                                   >  8:438822/1‑149 (MQ=255)
                                                                                                                                                                                                             gTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGt                                                                                                                                                                          >  8:193941/1‑149 (MQ=255)
                                                                                                                                                                                                             gTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGgtgt                                         <  2:131737/275‑1 (MQ=255)
                                                                                                                                                                                                                cTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAgg                                                                                                                                                                                                                  >  2:168363/1‑107 (MQ=255)
                                                                                                                                                                                                                   gCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGGTGTTGac                                     >  2:370720/1‑273 (MQ=255)
                                                                                                                                                                                                                    cccTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑ttt                                                                                                                                                                                                             <  4:168361/107‑1 (MQ=255)
                                                                                                                                                                                                                         ttGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTTACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGGTGTTGACCACCCTCa                             >  2:298340/1‑275 (MQ=255)
                                                                                                                                                                                                                             tGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGt                                                                                                                                                          >  7:204677/1‑149 (MQ=255)
                                                                                                                                                                                                                               gTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTc                                                                                                                                                                                                            >  4:402890/1‑97 (MQ=255)
                                                                                                                                                                                                                                tCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCaga                                                                                                                                                              <  7:300417/142‑1 (MQ=255)
                                                                                                                                                                                                                                tCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCaga                                                                                                                                                              >  8:300417/1‑142 (MQ=255)
                                                                                                                                                                                                                                 cTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTaaa                                                                                                                                                      <  7:193941/149‑1 (MQ=255)
                                                                                                                                                                                                                                   aTTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTa                                                                                                                                                                                                        <  2:402895/97‑1 (MQ=255)
                                                                                                                                                                                                                                   aTTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACCCAGAACGTTAACGTTCAGGTGTTGACCACCTTCAACGCGaa                      >  2:392520/1‑275 (MQ=255)
                                                                                                                                                                                                                                   aTTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTaa                                                                                        >  4:428924/1‑206 (MQ=255)
                                                                                                                                                                                                                                      gcgcGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGCGTCAGAGAGTTAAAGCGAAc                                                                                                                                                >  4:36090/1‑150 (MQ=255)
                                                                                                                                                                                                                                       cgcgGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGc                                                                                    <  2:428929/206‑1 (MQ=255)
                                                                                                                                                                                                                                          ggTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCa                                                                                                                                            <  2:36090/150‑1 (MQ=255)
                                                                                                                                                                                                                                            tACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGGTGTTGACCACCTTCAACGCGAACTTCTGGTTTCa          >  2:483853/1‑275 (MQ=255)
                                                                                                                                                                                                                                            tACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAAATTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GACGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGACTAACGTTAACGTTCAGGTGTTGACCCCCTTCAACGCCAACTTCTGGTTTCa          >  4:385396/1‑275 (MQ=255)
                                                                                                                                                                                                                                               gCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGa                                                                                                                                                                    <  8:8738/121‑1 (MQ=255)
                                                                                                                                                                                                                                               gCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGa                                                                                                                                                                    >  7:8738/1‑121 (MQ=255)
                                                                                                                                                                                                                                                      cTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGATTAACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGGTGTTGACCACCTTCAACCCGCACTTCCGGTTTCACTTCTACt  >  2:308695/1‑273 (MQ=255)
                                                                                                                                                                                                                                                       tATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAgg                                                                                                                                                                                                                  >  8:192189/1‑68 (MQ=255)
                                                                                                                                                                                                                                                       tATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAgg                                                                                                                                                                                                                  <  7:192189/68‑1 (MQ=255)
                                                                                                                                                                                                                                                          tAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTGCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGGACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGTATCTTCCAGAGTGTCGCGACGCAGAACGTTAACGTTCAGGTGTTGacc                                    <  1:7078/235‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGAt                                                                                                                           >  6:440680/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            aCGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGAt                                                                                                                           <  7:438822/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             cGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATg                                                                                                                          >  8:273270/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        ggAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGa                                                                                                                                                                                           <  7:380260/73‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ggAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGa                                                                                                                                                                                           >  8:380260/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                        ggAAGAGGGCGACCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCgg                                                                                                            >  7:356039/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        ggAAGAGGGCGACCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCgg                                                                                                            >  7:26784/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         gAAGAGGGCGACCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTaaa                                                                                                                                                                                                      <  8:299810/61‑1 (MQ=37)
                                                                                                                                                                                                                                                                         gAAGAGGGCGACCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTaaa                                                                                                                                                                                                      >  7:299810/1‑61 (MQ=255)
                                                                                                                                                                                                                                                                          aaGAGGGCGACCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGat                                                                                                          <  8:356039/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gagGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGGGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGatat                                                                                                        <  5:440680/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                           
TTTGCAGTTTGTGAACGTCAAAATCATTTTGTACGGTGTTACGGCACTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGC‑TACTTAAGCCTTACAGGC‑TTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGGTGTTGACCACCTTCAACGCGAACTTCTGGTTTCACTTCTACT  >  NC_000913/2715758‑2716276

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑