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→GAChflDputative 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

GTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTG  >  NC_000913/1192148‑1192576
                                                                                                                                                                                                                                           |                                                                                                                                                                                                    
gtgGTCTTTGCAGTACAGCCGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCg                                                                                                                                                                                                   <  2:9972/236‑1 (MQ=255)
 tgGGCTTTGCAGTTCAACAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCgg                                                                                                                                                                                                  <  1:62191/236‑1 (MQ=255)
     cTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAACGGGTTAATGCGGCGtt                                                                                                                                                                                              >  2:632150/1‑236 (MQ=255)
        tGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTa                                                                                                                                                                                            >  2:457053/1‑235 (MQ=255)
           aGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGcc                                                                                                                                                                                        <  2:506297/236‑1 (MQ=255)
               cagcagGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGAACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGACCTg                                                                                                                                                                                     >  1:557990/1‑235 (MQ=255)
                  cagGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTg                                                                                                                                                                                     >  1:649254/1‑231 (MQ=255)
                         aCCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGActgc                                                                                                                                                                             <  2:401159/233‑1 (MQ=255)
                            ggTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCAt                                                                                                                                                                       >  1:477085/1‑236 (MQ=255)
                                 ccTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAg                                                                                                                                                                   <  1:457053/235‑1 (MQ=255)
                                     aaTGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGc                                                                                                                                                              >  1:4244/1‑236 (MQ=255)
                                          gcgGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCg                                                                                                                                                                                  <  2:158059/211‑1 (MQ=255)
                                          gcgGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCg                                                                                                                                                                                  >  1:158059/1‑211 (MQ=255)
                                             ggCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAg                                                                                                                                                       <  2:4244/235‑1 (MQ=255)
                                             ggCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAg                                                                                                                                                       <  1:11568/235‑1 (MQ=255)
                                                    cGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCa                                                                                                                                               <  1:759726/236‑1 (MQ=255)
                                                             cacaTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGttt                                                                                                                                                                                             >  1:283625/1‑181 (MQ=255)
                                                             cacaTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGttt                                                                                                                                                                                             <  2:283625/181‑1 (MQ=255)
                                                              acaTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGAcgcg                                                                                                                                     <  1:49970/236‑1 (MQ=255)
                                                                 tCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAg                                                                                                                                  >  1:213576/1‑236 (MQ=255)
                                                                     cATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCg                                                                                                                                                                                                   <  1:215001/167‑1 (MQ=255)
                                                                     cATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCg                                                                                                                                                                                                   >  2:215001/1‑167 (MQ=255)
                                                                          atagCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAgg                                                                                                                                 >  1:436947/1‑228 (MQ=255)
                                                                          atagCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAgg                                                                                                                                 <  2:436947/228‑1 (MQ=255)
                                                                              cagcCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCt                                                                                                                      >  1:121170/1‑235 (MQ=255)
                                                                               agcCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGa                                                                                                                                                           <  2:509538/197‑1 (MQ=255)
                                                                               agcCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGa                                                                                                                                                           >  1:509538/1‑197 (MQ=255)
                                                                                                 ctttgcGGACTGTAAGTCGATGTGTCGAGGCGGCGGTTGCGGGCGGTGATCCGGGTTCCCCGGGGTGTCGAGCCCGCCTTCCGCTGAGGGGAGTTTGCGCTAAAGCACCATCAAGCTGGGTGTGTGGCGGTTTAAGTCGGCGTTTAAGCCCTGGCGACTGCGGGCATTGAGCCCGCCGAGCGGGGTTCCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCg                                                                                                  <  2:217722/231‑1 (MQ=255)
                                                                                                      cGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGCGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCTGCCCGACACGCAGGTTGGTTTCGCGACAGCCACAACACCACAGCGTCgtg                                                                                             >  2:667956/1‑234 (MQ=255)
                                                                                                      cGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGt                                                                                                 >  1:217722/1‑232 (MQ=255)
                                                                                                              aacgcGAAGTGCTCGAGTGGGCGTCGCAGCCCGTGGATCCGGTGCCCAAGGTGGCCGGGCGCTCCTTCGGAGGGGGAGGTTTTGCGCGAAAGCACCATCAAGCTGGTGGTGTGGGGGGTTAAGGCGGGGTTTAAGCCCGGGGGACGGCTGGCATTGGGCACGCCGGGCGGGGTTTCCAGCCCGACGCGAAGGTTGGTTTCGCTCCCGCCAAAAACCGCCAGCGTCGAGCTGGGGtt                                                                                     <  2:693300/232‑1 (MQ=255)
                                                                                                                  cGACGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCGGGTTGCCCAGGGTGCCGAGCGCCCTTGTCGCGGAGGAGAGTTTGCGCTAAAGCACAATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCGGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGTTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCat                                                                                  <  2:121170/235‑1 (MQ=255)
                                                                                                                     aGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCGGACGCGCAGGTTGGCTTCGCTACCGCCAACAACCGCCAGCGCCGAGCTGCGGTGCCTATCa                                                                              >  2:598337/1‑236 (MQ=255)
                                                                                                                       tgtTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGActgc                                                                                                                                                                             >  2:426035/1‑139 (MQ=255)
                                                                                                                       tgtTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGActgc                                                                                                                                                                             <  1:426035/139‑1 (MQ=255)
                                                                                                                                   cGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTtgag                                                                >  1:654481/1‑236 (MQ=255)
                                                                                                                                                            gCCCAGGGTGTGGAGCGCGCTGTGCGCTGAGGAGAGTTTGCGCTAAAGCCCCATCAAGCGGGGTGTGTGGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATGGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGTTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCaata                                                                           <  2:474395/202‑1 (MQ=255)
                                                                                                                                                                   gtgtCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAAGTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGt                                                                    >  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)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                    
GTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCACCATCAAGCTGAGTGTGTAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTG  >  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: