Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,192,380 T→G E84D (GAA→GAC)  hflD ← putative lysogenization regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,192,3800TG95.2% 126.4 / ‑4.2 42E84D (GAA→GAC) hflDputative lysogenization regulator
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (20/20);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTG  >  NC_000913/1192148‑1192576
                                                                                                                                                                                                                                           |                                                                                                                                                                                                    
gtgGTCTTTGCAGTACAGCCGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCg                                                                                                                                                                                                   <  2:9972/236‑1 (MQ=255)
 tgGGCTTTGCAGTTCAACAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCgg                                                                                                                                                                                                  <  1:62191/236‑1 (MQ=255)
     cTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAACGGGTTAATGCGGCGtt                                                                                                                                                                                              >  2:632150/1‑236 (MQ=255)
        tGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTa                                                                                                                                                                                            >  2:457053/1‑235 (MQ=255)
           aGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGcc                                                                                                                                                                                        <  2:506297/236‑1 (MQ=255)
               cagcagGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGAACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGACCTg                                                                                                                                                                                     >  1:557990/1‑235 (MQ=255)
                  cagGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCG‑TTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTg                                                                                                                                                                                     >  1:649254/1‑231 (MQ=255)
                         aCCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGActgc                                                                                                                                                                             <  2:401159/233‑1 (MQ=255)
                            ggTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCAt                                                                                                                                                                       >  1:477085/1‑236 (MQ=255)
                                 ccTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAg                                                                                                                                                                   <  1:457053/235‑1 (MQ=255)
                                     aaTGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGc                                                                                                                                                              >  1:4244/1‑236 (MQ=255)
                                          gcgGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCg                                                                                                                                                                                  <  2:158059/211‑1 (MQ=255)
                                          gcgGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCg                                                                                                                                                                                  >  1:158059/1‑211 (MQ=255)
                                             ggCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAg                                                                                                                                                       <  2:4244/235‑1 (MQ=255)
                                             ggCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAg                                                                                                                                                       <  1:11568/235‑1 (MQ=255)
                                                    cGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCa                                                                                                                                               <  1:759726/236‑1 (MQ=255)
                                                             cacaTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGttt                                                                                                                                                                                             >  1:283625/1‑181 (MQ=255)
                                                             cacaTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGttt                                                                                                                                                                                             <  2:283625/181‑1 (MQ=255)
                                                              acaTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGAcgcg                                                                                                                                     <  1:49970/236‑1 (MQ=255)
                                                                 tCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAg                                                                                                                                  >  1:213576/1‑236 (MQ=255)
                                                                     cATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCg                                                                                                                                                                                                   <  1:215001/167‑1 (MQ=255)
                                                                     cATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCg                                                                                                                                                                                                   >  2:215001/1‑167 (MQ=255)
                                                                          atagCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAgg                                                                                                                                 >  1:436947/1‑228 (MQ=255)
                                                                          atagCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAgg                                                                                                                                 <  2:436947/228‑1 (MQ=255)
                                                                              cagcCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCt                                                                                                                      >  1:121170/1‑235 (MQ=255)
                                                                               agcCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGa                                                                                                                                                           <  2:509538/197‑1 (MQ=255)
                                                                               agcCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGa                                                                                                                                                           >  1:509538/1‑197 (MQ=255)
                                                                                                 ctttgcGGACTGTAAGTCGATGTG‑TCGAGGCGGCGGTTGC‑GGGCGGTGATCCGGGTTCCCCGGGGT‑GTCGAGCCCGCCTTCCGCTGAGGGGAGTTTGCGCTAAAGCACCATCAAGCTGGGTGTGTGGCGGTTTAAGTCGGCGTTTAAGCCCTGGCGACTGCGGGCATTGAGCCCGCCGAGCGGGGTTCCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCg                                                                                                  <  2:217722/231‑1 (MQ=255)
                                                                                                      cGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGCGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCTGCCCGACACGCAGGTTGGTTTCGCGACAGCCACAACACCACAGCGTCgtg                                                                                             >  2:667956/1‑234 (MQ=255)
                                                                                                      cGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGt                                                                                                 >  1:217722/1‑232 (MQ=255)
                                                                                                              aacgcGAAGTGCTCGAGTGGGC‑GTCGCAGCCCGTGGATCC‑GG‑TGCCCAAGGTGGCCGGGCGCTCCTTCGGAGGGGGAGGTTTTGCGCGAAAGCACCATCAAGCTGGTGGTGTGGGGGGTTAAGGCGGGGTTTAAGCCCGGGGGACGGCTGGCATTGGGCACGCCGGGCGGGGTTTCCAGCCCGACGCGAAGGTTGGTTTCGCTCCCGCCAAAAACCGCCAGCGTCGAGCTGGGGtt                                                                                     <  2:693300/232‑1 (MQ=255)
                                                                                                                  cGACGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCG‑GGTTGCCCAGGGT‑GCCGAGCGCCCTTGTCGCGGAGGAGAGTTTGCGCTAAAGCACAATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCGGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGTTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCat                                                                                  <  2:121170/235‑1 (MQ=255)
                                                                                                                     aGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCGGACGCGCAGGTTGGCTTCGCTACCGCCAACAACCGCCAGCGCCGAGCTGCGGTGCCTATCa                                                                              >  2:598337/1‑236 (MQ=255)
                                                                                                                       tgtTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGActgc                                                                                                                                                                             >  2:426035/1‑139 (MQ=255)
                                                                                                                       tgtTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGActgc                                                                                                                                                                             <  1:426035/139‑1 (MQ=255)
                                                                                                                                   c‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTtgag                                                                >  1:654481/1‑236 (MQ=255)
                                                                                                                                                            gCCCAGGGT‑GTGGAGCGCGCTGTGCGCTGAGGAGAGTTTGCGCTAAAGCCCCATCAAGCGGGGTGTGTGGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATGGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGTTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCaata                                                                           <  2:474395/202‑1 (MQ=255)
                                                                                                                                                                   gt‑gtCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGt                                                                    >  1:474395/1‑202 (MQ=255)
                                                                                                                                                                               gcCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAg                       >  1:430476/1‑236 (MQ=255)
                                                                                                                                                                                                         cAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTg  >  1:453291/1‑231 (MQ=255)
                                                                                                                                                                                                                                    gggTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCg                     <  1:725296/185‑1 (MQ=255)
                                                                                                                                                                                                                                    gggTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCg                     >  2:725296/1‑185 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                    
GTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGC‑GTTGCAGGCCGTTGATCC‑GGTTGCCCAGAGT‑GTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTG  >  NC_000913/1192148‑1192576

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑