Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913803,8690GT32.5% 60.3 / 12.9 40V123L (GTG→TTG) ybhJputative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base G (4/23);  new base T (12/1);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 3.68e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCTGGC  >  NC_000913/803623‑804119
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                              
aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGcc                                                                                                                                                                                                                                                            <  1:438179/251‑1 (MQ=255)
   cactctACGCATAATACGAGCAGAAATATGAATAAACATAAAGTTAAGTTTGATTAATTAGCCTCTCAAGATATTTGCTTTTTTGGTATTGTACAGACAGGTAAGGCTGCCGGTATTGACCGTGTCCCGGTGCCGGATGTGCTGACCAACTGCGATAACTCTCTCTGCGCCGTCGGCGGTACTGTTAGCGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCTGTATTTTTGTGCCTcc                                                                                                                                                                                                                                                         <  2:418694/249‑1 (MQ=255)
     ctctTCTCATAATTCGTCCGGAAATATGGATAAACTTAAAATTAAGTTTTATTCATTAGCCTCTCACGATATTACCTTTGTTGTTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTc                                                                                                                                                                                                                                                          <  2:103924/248‑1 (MQ=255)
     cactTCGCATAATACGTCCGGAAATATGTATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACTATATTTCCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCTCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTTATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGTCGGTATTTTTGTGCCTcc                                                                                                                                                                                                                                                         <  1:29086/247‑1 (MQ=255)
            cATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGc                                                                                                                                                                                                                                                <  1:139792/251‑1 (MQ=255)
             ataataCGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCTCCGCATAtttcg                                                                                                                                                                                                                                               >  1:454300/1‑248 (MQ=255)
                              atGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAatatat                                                                                                                                                                                                                              <  2:501069/251‑1 (MQ=255)
                                                 aGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCCCCTTTTTCGGTCCTCCAACAATATATGTGTGTGATGGTGGGAggg                                                                                                                                                                                                         >  1:77584/1‑250 (MQ=255)
                                                 aGTTTGATTCATTAGCCTCTCACGATAGTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAggc                                                                                                                                                                                                         <  1:277308/251‑1 (MQ=255)
                                                  gTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCTCCGCCTATTGGGGTCATCCAACAatatatttgtggggtggtgggagggg                                                                                                                                                                                                          >  1:457929/1‑231 (MQ=255)
                                                   tttGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCCCCTATTGCGGGCATCCCTCAATATATTCGTGGGATGATGGCAggcgg                                                                                                                                                                                                       >  1:239345/1‑251 (MQ=255)
                                                           aTTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCCTATTGCGGGCATCCCTCAATATTAGCGGGAGATGATTGCAggcggcgggcca                                                                                                                                                                                                >  1:392111/1‑246 (MQ=255)
                                                             tAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAAt                                                                                                                                                                                              <  1:258502/250‑1 (MQ=255)
                                                              aGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGa                                                                                                                                                                                            <  1:418812/251‑1 (MQ=255)
                                                                 ctctCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCGCTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTTTCTTCCCCTATTTCGGGCCTCCATCAAAATATGCGGGGGATGGAGGCAGGCGGGGGCaaaaagaatc                                                                                                                                                                                         >  1:135771/1‑245 (MQ=255)
                                                                         atatTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCCTTTTTCGGGCATCCATCCATCTATGCGGTGGAGTGTGGCCGGCGGCGGCCAAATGATCCTCgggggc                                                                                                                                                                                 >  1:17654/1‑248 (MQ=255)
                                                                                  tggttgGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCca                                                                                                                                                                         <  2:239345/248‑1 (MQ=255)
                                                                                       ggTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTTCCCATATTTCGGTCATTCCTTAATATTTTCGGGGGATGATTGGAGGCGGGGGCAAAAAGGTCCTCGGGGTAAAAAGCCACAcccc                                                                                                                                                                  >  1:139458/1‑250 (MQ=255)
                                                                                        gTATTGTACACACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGt                                                                                                                                                                 <  1:138511/251‑1 (MQ=255)
                                                                                          aTTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGTTGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGtt                                                                                                                                                                <  1:291428/250‑1 (MQ=255)
                                                                                               aCAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTTCCTCCGCATATTTCGGGCATCCATCCATATTTTCGGGAGATGATGGCAGGCGGGGGGAAAATGATCCTCGGGTCCGACAGCCCAACCCCTTACgggg                                                                                                                                                          >  1:135727/1‑249 (MQ=255)
                                                                                                         aagtcttCCGTTATGGAACGTTTCCCGCTGCCGTATGGGCTGACCATCTTCCATAACTCACTCTTCGCCGTCGGCTGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTGCGGTCATCCAGCAATATATGCGTGAGATGATGGCAGTCGGCGGCAAAATGATCCTCGGATCAGACAGCCCCACCCTTTACGGTGCATTAggg                                                                                                                                                  <  2:279938/243‑1 (MQ=255)
                                                                                                                                     tGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTTCGGGCATCCATCCATATATGCGTGAGATGATTGCAGGCGGCGGCAAAATGATCCTCCGGGCAGACAGCCACACCCGTTACGGGGCATTTGGGACAATGGCAGGCGGGGAGGGcggcgggg                                                                                                                      >  1:333817/1‑247 (MQ=255)
                                                                                                                                             gtgcgtaCCAACTGCCATAACTCACGCTGCGCCGGCGGCGGCACTAGTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCAGCAAAATGATCCTCGGGTCAGACAGCCACTCCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTaa                                                                                                             <  2:333971/244‑1 (MQ=255)
                                                                                                                                                     cAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCCCCTATTTCGGGCATCCCTCAAAATATGCGTGGGATGAATGCAGGGGGCGGGAAAATGATACCCGGGGCCGACCGCCACACCCCGTACGGGGCAATAGGGGCAATGGCAGTCGGGGAGGGGGGGGGGGGGGTGGGAAAACAgctgg                                                                                                    >  1:349703/1‑250 (MQ=255)
                                                                                                                                                           ccATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACATCTGCTTAATg                                                                                              <  1:54242/251‑1 (MQ=255)
                                                                                                                                                            cATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAAtta                                                                                             <  1:86834/251‑3 (MQ=255)
                                                                                                                                                                   cactcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCCTATTTCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAAAGATCCTCGGGTCAGAAAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGGGAGGGGGGCGGTGAGGTGGGAAAACAGCTGCTTAATGACACCTgg                                                                                      >  1:206680/1‑251 (MQ=255)
                                                                                                                                                                       ctGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTGTGCCTCCGCATATTTCGGGCATTCATCAATATATGCGGGAGATGATGGGAGGGGGCGGGAAAATGAACCTCGGGGGAGAAAGCCACACCCGGTACGGGGGATTTGGGGCAATGGGAGGCGGTGAGGGCGGGGGGGGGGTGGGAAAAAAGATGCTTAATGAAAACTggggg                                                                                   >  1:236941/1‑248 (MQ=255)
                                                                                                                                                                            ccGTCGTCGGCACTATTAACTGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCTTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTTGGATATCGAc                                                                             <  2:139458/251‑1 (MQ=255)
                                                                                                                                                                                                   gatgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGc                                                      <  1:211179/250‑1 (MQ=255)
                                                                                                                                                                                                       tccATGTTTTTGGTGTAGCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCt                                                  <  2:352472/249‑1 (MQ=255)
                                                                                                                                                                                                            gTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCCCATATTTCGGTCATCCCTCAATATATGCGTGAGGTGATGGCAGGGGGCGGCAAAAAGAACCCCCGGTCAAACAGCCACACCCGCTACGGGTCATTAGGGAAAAAGGGAGGCGGGGAGGGGGGGGGGGGGGTTGTAAAAAAGCCGCTTAATGACACCCGGGGTATTGACTATTCCggggggggtgcggggcattttgcccg                                             >  1:462733/1‑228 (MQ=255)
                                                                                                                                                                                                                ttGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTTCGGGCATCCATCAATATATGCGTGAGATGATGGCAGGGGGCGGCAAAATGATCCTCGGGGCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGGGAGGGCGGCGGTGAGGTGGTAAAACAGCTGCTTAATGAAAACCGGGGTATCGACTATCCGGGCGTGGTTGCGGGGCATCTGACCGGaaa                                         >  1:384220/1‑250 (MQ=255)
                                                                                                                                                                                                                                ccAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCAGCTGACCGGAAAACCTGCGCCGTATGTgg                        <  2:315551/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGCGTTATGGCGGTATTTTTGTTCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTggg                       <  2:463841/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTTCGGCCATCCATCCATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATACCCGGGTCAGACAGCCACCACCCTTACGGAGCATTAGGGAAAATGGCAGGCGGCGAGGGCGGGGGCGAGTTGGGAAAAAAGCTGCCTAATGAAACCCGGGGTATCGACCATTCGGGCGCGGGTGCGGTGCATCTTGCCGGAAAACCTTCCCCGTATGCggg                       >  1:394731/1‑251 (MQ=255)
                                                                                                                                                                                                                                       tATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCcac                  <  1:99717/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cgttATTTTTGTGCCTCCGCATAGTGCGTTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGAtg            <  2:427092/248‑1 (MQ=255)
                                                                                                                                                                                                                                             ggTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGAtgt           <  2:454300/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      gtgCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCTGGc  >  1:313746/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                              
ATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCTGGC  >  NC_000913/803623‑804119

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