Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,716,033 C→A 10.8% intergenic (+23/+33) 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,0330CA10.8% 227.7 / 17.3 76intergenic (+23/+33)eamB/grcAcysteine and O‑acetylserine exporter/autonomous glycyl radical cofactor
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base A (3/5);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-03

CAGTTTGTGAACGTCAAAATCATTTTGTACGGTGTTACGGCACTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGGTGTTGACCACCTTCAACGCGAACTTCTGGTTTCACTTCTACT  >  NC_000913/2715762‑2716276
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        
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)
                                                                                                                                              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                                                                                                                                                                                           >  8:172817/1‑145 (MQ=255)
                                                                                                                                                                                            ccAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGa                                                                                                                                                                                           <  7:172817/145‑1 (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                                                                                                                                                              >  8:300417/1‑142 (MQ=255)
                                                                                                                                                                                                                            tCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCaga                                                                                                                                                              <  7:300417/142‑1 (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                                                                                                                                                                                           >  8:380260/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                    ggAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGa                                                                                                                                                                                           <  7:380260/73‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ggAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCgg                                                                                                            >  7:26784/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    ggAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCgg                                                                                                            >  7:356039/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)
                                                                                                                                                                                                                                                                              gCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTcc                                                                                                  <  8:204677/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        
CAGTTTGTGAACGTCAAAATCATTTTGTACGGTGTTACGGCACTGTCGACGTTTGTTCTGCCGCAAACACAGGCGTTAAGCTGGGTAGTTGGCGTCAGCGTTTTGCTGGCGATGATTGGGACGTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCA‑‑‑GCTGCGGATATTTTTCCGGATGCTTAACTGCATCTTCCAGAGTTTCGCGACGCAGAACGTTAACGTTCAGGTGTTGACCACCTTCAACGCGAACTTCTGGTTTCACTTCTACT  >  NC_000913/2715762‑2716276

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