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

TTTGCAGTTTGTGAACGTCAAAATCATTTTGTACGGTGTTACGGCACTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑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)
                                         cGGCACTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGTTCGCCCTCTTCCGCTTAGTCACTTGCTACTTAATCCTTACAgg                                                                                                                                                                                                                  >  2:230100/1‑274 (MQ=255)
                                          ggCACTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGGGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCt                                                                                                                                                                                                               <  2:339594/275‑1 (MQ=255)
                                            cccTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAAGTGGGGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCttt                                                                                                                                                                                                             <  2:50383/273‑1 (MQ=255)
                                              cTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCGGGGTAGTTGGCGCCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTTCTTTCTTGCCGGTCTACTTCGCGGTACGCATCTTCTATTAACGAAACAACGCGGAAGAGGTCGCTGTCTGCCGTGTAGTCATTGCGTACCTAAGCCTTACAGACtttta                                                                                                                                                                                                           >  4:437534/1‑273 (MQ=255)
                                                  cGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCctc                                                                                                                                                                                                                                                      <  2:284993/230‑1 (MQ=255)
                                                  cGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGGGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTa                                                                                                                                                                                                                      <  2:437539/261‑1 (MQ=255)
                                                      gTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGAGTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCGTTTTCTATTAAGGAAAAAAAGCGGAAGAGGTCGGCCTCTTCTGCTTAGTAACTGGCTACTTAAGCCTTACAGGCTTTCAGTAAAgt                                                                                                                                                                                                    >  1:16010/1‑273 (MQ=255)
                                                          agacagCCGCAAACACAGGCGTTAAGCGGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACg                                                                                                                                                                                                <  2:471995/269‑1 (MQ=255)
                                                           ttCTGCCGCAAACACACGCGTTAAGCTGGTTAGTTGGCGTTAGCGCTTCGCTGGCGATGATTGGGACCTTTGGCAATGTGTGCTGGGCGCTCGCGTGGCCTCTGTTTGCGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTTGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAATGCGCAATAGGTCGCCCTCTTCCGCATAGTAACTTGCTACTTGAGCACTAAGGCATTTCAGTaaggcaaaac                                                                                                                                                                                               >  4:262816/1‑266 (MQ=255)
                                                                                   ggggggTAGTTGGCGGTAGCGTTTGGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTGCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATGGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCg                                                                                                                                                                      <  4:72265/272‑1 (MQ=255)
                                                                                      gggTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGNGTTTNGGCAATGTGTGCTGGGCGCTGGCNGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGcc                                                                                                                                                                                                                                                         <  4:505751/191‑1 (MQ=255)
                                                                                      gggTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGACAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTCAAGTTACGACCGATAACGTCGGGCTGCTGTTCCTGa                                                                                                                                                                    >  1:40926/1‑274 (MQ=255)
                                                                                      gggTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTTCTGTTCCGGa                                                                                                                                                                    >  2:122363/1‑274 (MQ=255)
                                                                                       ggTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTAGCCCTCTTCCGCGTAGTAAATTGCTACTTAAGCCTTaaa                                                                                                                                                                                                                    >  1:85197/1‑224 (MQ=255)
                                                                                          aGTTGGCTTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTTGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCa                                                                                                                                                                <  1:99698/274‑1 (MQ=255)
                                                                                                       gTTTGGGTGGCGATTATTGGTACGTTTGGCAAGTTGTGATGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTATATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCg                                                                                                                                                   <  1:136061/274‑1 (MQ=255)
                                                                                                          ttGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCGGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGCTAACGTCGCGCTGCTGTTCCGGAGTTGGATGATTAAAGCAAAca                                                                                                                                               >  1:114001/1‑274 (MQ=255)
                                                                                                                               tttGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGCCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGAGGAAGAGGTCGCACACTTCCGCTTGGTAACTGGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAGACGAAACGCAGCACCGGAAACACGGGTg                                                                                                                          >  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)
                                                                                                                                                  cgcTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAACGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAAAGTTAAAGCGAACTGCATAACCGGATCCAGGGACGGTCA‑‑‑GGTGAGGATGtttt                                                                                                     >  2:491998/1‑274 (MQ=255)
                                                                                                                                                    cgggCGGGGCAGCTGTTTCAGCGATTGTTTCGCCAGTAGGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCGGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTcc                                                                                                  <  1:53468/273‑1 (MQ=255)
                                                                                                                                                    cTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGcc                                                                                                                                                                                                                                                         <  2:438171/129‑1 (MQ=255)
                                                                                                                                                     tGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCTGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATCACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCg                                                                                                 >  2:275366/1‑275 (MQ=255)
                                                                                                                                                       ggggtgaattggTTTCAGGGTTTGTGTNCNAAGTATGGTCGCCAGTTAAATAGCGGTATTGCCCTATTGCTGATTNATTGCTCGGTATGCAATTTCTATTAACGAAAAAAAGCGGAAGAGGTAGGACTATGCCGCTGAGTTACTGNCTACCTAAGCCGTGCAGGCTTTCAGTAAAGGTACGAGCGATAACGGCGCGCTGCTGTTCCGGAGTCAGAGAGTGAAAGCGAACTGCAGAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCgg                                                                                                <  3:40926/263‑1 (MQ=255)
                                                                                                                                                               tCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGcc                                                                                                                                                                                                                         >  8:468132/1‑149 (MQ=255)
                                                                                                                                                                 tGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCtt                                                                                                                                                                                                                       >  6:55106/1‑149 (MQ=255)
                                                                                                                                                                    ttCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGACGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACa                                                                                                                                                                                                                    <  7:124624/149‑1 (MQ=255)
                                                                                                                                                                          gATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCttt                                                                                                                                                                                                             <  8:351604/149‑1 (MQ=255)
                                                                                                                                                                             tGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGCGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATAGTTTTCCGTATGCTTAACTGCATCTTACagag                                                                         >  1:72031/1‑275 (MQ=255)
                                                                                                                                                                                 tCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTaaa                                                                                                                                                                                                      <  7:468132/149‑1 (MQ=255)
                                                                                                                                                                                          tGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGCa                                                             <  2:265600/274‑1 (MQ=255)
                                                                                                                                                                                           ggTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCg                                                                                                                                                                                            <  5:55106/149‑1 (MQ=255)
                                                                                                                                                                                                ccAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGa                                                                                                                                                                                           <  7:172817/145‑1 (MQ=255)
                                                                                                                                                                                                ccAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGa                                                                                                                                                                                           >  8:172817/1‑145 (MQ=255)
                                                                                                                                                                                                    ttAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTcg                                                                                                                                                                                   >  8:438822/1‑149 (MQ=255)
                                                                                                                                                                                                             gTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGt                                                                                                                                                                          >  8:193941/1‑149 (MQ=255)
                                                                                                                                                                                                             gTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGgtgt                                         <  2:131737/275‑1 (MQ=255)
                                                                                                                                                                                                                cTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAgg                                                                                                                                                                                                                  >  2:168363/1‑107 (MQ=255)
                                                                                                                                                                                                                   gCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGGTGTTGac                                     >  2:370720/1‑273 (MQ=255)
                                                                                                                                                                                                                    cccTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCttt                                                                                                                                                                                                             <  4:168361/107‑1 (MQ=255)
                                                                                                                                                                                                                         ttGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTTACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGGTGTTGACCACCCTCa                             >  2:298340/1‑275 (MQ=255)
                                                                                                                                                                                                                             tGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGt                                                                                                                                                          >  7:204677/1‑149 (MQ=255)
                                                                                                                                                                                                                               gTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTc                                                                                                                                                                                                            >  4:402890/1‑97 (MQ=255)
                                                                                                                                                                                                                                tCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCaga                                                                                                                                                              <  7:300417/142‑1 (MQ=255)
                                                                                                                                                                                                                                tCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCaga                                                                                                                                                              >  8:300417/1‑142 (MQ=255)
                                                                                                                                                                                                                                 cTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTaaa                                                                                                                                                      <  7:193941/149‑1 (MQ=255)
                                                                                                                                                                                                                                   aTTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTa                                                                                                                                                                                                        <  2:402895/97‑1 (MQ=255)
                                                                                                                                                                                                                                   aTTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACCCAGAACGTTAACGTTCAGGTGTTGACCACCTTCAACGCGaa                      >  2:392520/1‑275 (MQ=255)
                                                                                                                                                                                                                                   aTTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTaa                                                                                        >  4:428924/1‑206 (MQ=255)
                                                                                                                                                                                                                                      gcgcGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGCGTCAGAGAGTTAAAGCGAAc                                                                                                                                                >  4:36090/1‑150 (MQ=255)
                                                                                                                                                                                                                                       cgcgGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGc                                                                                    <  2:428929/206‑1 (MQ=255)
                                                                                                                                                                                                                                          ggTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCa                                                                                                                                            <  2:36090/150‑1 (MQ=255)
                                                                                                                                                                                                                                            tACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGGTGTTGACCACCTTCAACGCGAACTTCTGGTTTCa          >  2:483853/1‑275 (MQ=255)
                                                                                                                                                                                                                                            tACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAAATTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GACGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGACTAACGTTAACGTTCAGGTGTTGACCCCCTTCAACGCCAACTTCTGGTTTCa          >  4:385396/1‑275 (MQ=255)
                                                                                                                                                                                                                                               gCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGa                                                                                                                                                                    <  8:8738/121‑1 (MQ=255)
                                                                                                                                                                                                                                               gCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGa                                                                                                                                                                    >  7:8738/1‑121 (MQ=255)
                                                                                                                                                                                                                                                      cTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGATTAACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGGTGTTGACCACCTTCAACCCGCACTTCCGGTTTCACTTCTACt  >  2:308695/1‑273 (MQ=255)
                                                                                                                                                                                                                                                       tATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAgg                                                                                                                                                                                                                  >  8:192189/1‑68 (MQ=255)
                                                                                                                                                                                                                                                       tATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAgg                                                                                                                                                                                                                  <  7:192189/68‑1 (MQ=255)
                                                                                                                                                                                                                                                          tAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTGCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGGACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGTATCTTCCAGAGTGTCGCGACGCAGAACGTTAACGTTCAGGTGTTGacc                                    <  1:7078/235‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGAt                                                                                                                           >  6:440680/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            aCGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGAt                                                                                                                           <  7:438822/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             cGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATg                                                                                                                          >  8:273270/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        ggAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGa                                                                                                                                                                                           <  7:380260/73‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ggAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGa                                                                                                                                                                                           >  8:380260/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                        ggAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCgg                                                                                                            >  7:356039/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        ggAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCgg                                                                                                            >  7:26784/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         gAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTaaa                                                                                                                                                                                                      <  8:299810/61‑1 (MQ=37)
                                                                                                                                                                                                                                                                         gAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTaaa                                                                                                                                                                                                      >  7:299810/1‑61 (MQ=255)
                                                                                                                                                                                                                                                                          aaGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGat                                                                                                          <  8:356039/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gagGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGGGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGatat                                                                                                        <  5:440680/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                           
TTTGCAGTTTGTGAACGTCAAAATCATTTTGTACGGTGTTACGGCACTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑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: