Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,0490AC28.8% 70.3 / 11.6 52G13G (GGT→GGGykfBprotein YkfB
Reads supporting (aligned to +/- strand):  ref base A (24/13);  new base C (1/14);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.43e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATGTACTGCTGGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAAGGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCC  >  NC_000913/265800‑266290
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
aatGTACTGCTGGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAAGGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTAGTCCGTATTGCCCTGACCGTAACCCGGCTGCCAGAACCCCGGGATAGCAGGGGCGTTAAAGTTCGCGAGCAGTaca                                                                                                                                                                                                                                                   <  2:386778/250‑1 (MQ=255)
 atGTACTGCTGGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAAGGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACAcc                                                                                                                                                                                                                                                 <  1:549529/251‑1 (MQ=255)
  tGTACTGCTGGTGATCCATCACTACGGTGATAGTAGTCTCGTCATCGTACGAACTGACCACCTTATCAACGGCACGGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGATGGTAAACGACTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACAccc                                                                                                                                                                                                                                                >  1:534880/1‑251 (MQ=255)
   gTACTGCTGGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAAGGGTAAGAAAGACTGAATGCTGATAGAATCCATGCTGGTCCGGGTTCCCTGTGCAGTGGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTAGTCCGTATTGCCCTGACCGTAACCCTGCGGTCAGAACCCGGGGATAGCAGATGCATTAAAGCTCGCGAGCAGTACAccc                                                                                                                                                                                                                                                <  1:210638/250‑1 (MQ=255)
                  cATCACTACGGTGATAGTAGTGACGTCATCGTGCGAACTAACCATATTGTCATCGGCAAGGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTACGGGTTCCCTGTGCAGTTGATGGTAAACGTCTTTACGCTGGCTTCGGTCACGATGTATTCCGTATTGCACTGACCGTAACCCTGCTGCCAGAACACCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGAc                                                                                                                                                                                                                                >  1:84494/1‑251 (MQ=255)
                         aCGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAAGGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTg                                                                                                                                                                                                                          >  2:491294/1‑250 (MQ=255)
                              aaTAGTAGTGTCTTCATCGTGCGAACTGCCCATCTTGTCATCGGCAAGGGTAAGAAAGACTGAATTCTGATAAAACCCATTCTGTTCCGGTTTCCCTGTGAAGTTGATGGTAAACGTCTTTCCGCGGCCTTCGGTCACGCTGTATTCCGTATTGCCCTGCCCGTACCCCTGCTGCCAGACCCCGGGTATAGCAGAGGATTTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGCCTTAGTGAAAGt                                                                                                                                                                                                                     <  1:407427/249‑1 (MQ=255)
                                      gggcgtaaTCGTGCGACCTGACCATCTTTTCATGGGCAAGGGTAAGAAAGACTGAATGCTAAAAGAAACCATTCTGGCCCGGGTTCCGTTTGCAGTGGATGGAAACCTTCTTTCCGCGGGCTTCGTTCACGCTGTATTCCGGTTTGCCCTGACCGAAACCCTGCTGCCAGACCCCCGGTATAGCAGGGGGATTAAAGCTCGCGGGAGGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCAtt                                                                                                                                                                                                            <  1:284608/244‑1 (MQ=255)
                                            cATCGTGCGAACTGACCATCTTGTCATCGGCAAGGGTAAGCAAGACTGAATGCTGATCGCAACCATTCTGGTCCGGGATCCCTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGtc                                                                                                                                                                                                      >  2:49857/1‑251 (MQ=255)
                                                  gCGACCTGACCATCTTGTCATCGGCAAGGGTAAGAAAGACTGAATGCTGATAGAACCCATTCTGGTCCGGGTTCCCTGTGCAGTTGATGGAAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCtt                                                                                                                                                                                                 <  1:221843/250‑1 (MQ=255)
                                                     accggACCATCTGGCCACCGGAAAGGGAAAGAAAGACTGAATGCGGATAGAACCCATTCTGGTCGGGGTTCCCTGTGCAGTTGATGGAAAACGTCTTCCCGCGGGCTGCGGTCACGCTGTATCCCGTATTGCCCTGACCGTAACCCTGCTGCCAGACCCCGGGGATAGCAGAGGCATAAAATCTCGCGAGCAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTc                                                                                                                                                                                             <  1:226649/247‑1 (MQ=255)
                                                     aaCTGACCATCTTGTCATCGGCAAGGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTTTTCCTTATTGCCTTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTGGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTc                                                                                                                                                                                             <  1:509346/251‑1 (MQ=255)
                                                                 tGTCATCGGAAAGGGTAAGAAAGACTGAAGGCGGATAGAAACCATTCTGGTCCGGGTTCCCTGTGAAGTTGATGGAAACCGTCTTCCCGCTGGCTTCGGTCCCGCTGTATTCCGTAGTGCCCTGACCGAAACCCGGCTGCAAAACCCCCGGGATAGCAGAGCCATAAAAGCTCGCGAGAAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTAt                                                                                                                                                                                   <  1:364430/249‑1 (MQ=255)
                                                                        ggCAGGGGTAAAAAAGCCTGAATGCTGATAGAACCCTTCCGGTCCCGGGTTCCCTTTGCATTTGATGGTAGACGTTTTGCCGTGGGCTTCGGTCCCGCTGTTTCCCTTTTTGCCCTGCCGGTACCCTTGCCGACAGACCCCCGGAATAGCAGTGGCATTAAAGCTCGCAGCAAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGGATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTg                                                                                                                                                                          <  1:64444/251‑1 (MQ=255)
                                                                                   gaaagaCTGAATGCTGATAGAAACCATTCTGGTCCGGGTACCATGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCg                                                                                                                                                                >  1:499586/1‑250 (MQ=255)
                                                                                      agaCTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGAAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGAGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAggg                                                                                                                                                            >  1:292860/1‑251 (MQ=255)
                                                                                                  gATAGAAACCATTCTGGTCCGGGTTCCCGGTGCAGTTGATGGAAAACGTCTTCCCGCTGGCTTCGGTCCCGCTGTATTCCTTATTGCCCTGCCCTTAACCCTGCTGCCAGACCCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACCCCCGCCAGCATAAACCGa                                                                                                                                                                                                                                 <  1:520977/170‑1 (MQ=255)
                                                                                                  gATAGAAACCATTCTGGTCAGGGTTCCCTGTGCAGTTGATGGTAAACGTCTTTCTGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGa                                                                                                                                                                                                                                 >  2:520977/1‑170 (MQ=255)
                                                                                                      gAAACCATTCTGTCCCGGTTTCCCTGTGCAGTTGTTGGTAACCGTCTTCCCGCTGGCTTCGGTCGCGTTGTTTCCCTTATTGCCTTGCCCGTACCCCTGCTGCCAGAACCCCGGGATAGCAGAGCAATTAAAGCTCGCGAGAGTTCCACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGcccc                                                                                                                                            <  1:320154/251‑1 (MQ=255)
                                                                                                          ccATTCTGGTCCGGTTTCCCTGTGCAGTTGATGGAAAACGTCTTTCCGCTGGCTGGGGTCACGCTGTATTCCGTATTGCCCTGACCGTACCCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCCCGCGAGAAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCa                                                                                                                                        <  1:491294/251‑1 (MQ=255)
                                                                                                          ccATCCTGGTCGGGGTTCCCGGGGCTGTTGGGGTAACACGTCTTCCGGCGGTCTTCGGTCACGTTGTATTCCGTATGCCCCGGCCGGAACCCCTGCTACCAGACCCCGGGGATAGCAGAGTCATTAAAGCCCGCGAGCGGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCa                                                                                                                                        <  2:125274/251‑1 (MQ=255)
                                                                                                               cgggcccGGGTCCCCTTGGCATTTGATGGAAACCGTCTTTCCGCTGGCTTCGGTCCCGCTGTATCCCGTATTCCCCTGACCGTACCCCGGCTGCCAGAACCCCGGGATAGCAGAGGAATTAAAGCCCGCGAGCAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGttt                                                                                                                                    <  1:366793/245‑1 (MQ=255)
                                                                                                                  gTCCGGGTTCCCTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCg                                                                                                                                >  1:289993/1‑251 (MQ=255)
                                                                                                                       ggggtccttgtgCAGTTGATGGTAAACTTCTTTCCGCTGGCTTCGGTCACGCTGTATCCCGTATTGCCCTGACCGGACCCCGGCTCCCAGAACCCGGGGAAGGCAGAGGCATTAAAGCTCGCTGGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCGGATTGTCAGTGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCAt                                                                                                                            <  1:49857/243‑1 (MQ=255)
                                                                                                                        ggTTCCTTGTGCTGTTGATGGTAACCGTCTTTCCGCGGGCTTCGGTCCCGCTGTATTCCGTTTTGCCCTGACCGTAACCCTGCTGCCAGACCCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCGGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATg                                                                                                                           <  2:21664/251‑1 (MQ=255)
                                                                                                                         gTTCCCTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGa                                                                                                                                                                                                                         >  1:364798/1‑156 (MQ=255)
                                                                                                                          ttCCCTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATggg                                                                                                                         >  2:52151/1‑251 (MQ=255)
                                                                                                                          ttCCCTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATggg                                                                                                                         >  1:215868/1‑251 (MQ=255)
                                                                                                                          ttCCCTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGAATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATggg                                                                                                                         >  2:249691/1‑251 (MQ=255)
                                                                                                                            cccTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCGCGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCg                                                                                                                       >  2:406321/1‑251 (MQ=255)
                                                                                                                            cccTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCg                                                                                                                       >  2:440523/1‑251 (MQ=255)
                                                                                                                                 tgCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTc                                                                                                                  >  1:314351/1‑251 (MQ=255)
                                                                                                                                                   gTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTgg                                                                                                 >  1:415096/1‑251 (MQ=255)
                                                                                                                                                      tttCCGCTGGCTTCGGTCCCGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGtat                                                                                              >  1:482028/1‑251 (MQ=255)
                                                                                                                                                         ccgctgctttCGGTCACGCTGTCTTCCGTATTGCCCTGACCGTACCCCTGCTGCCAGACCCCCGGGATGGCAGAGGATTTAAGGCTCGCGAGCAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTTCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACg                                                                                           <  1:112957/243‑1 (MQ=255)
                                                                                                                                                                  gcggccACGCTGTATCCCTTATTGCCCGACCCGTAACCCTGCTGCCAGACCCCCGGGAAAGAAGAGGAATTAAAGCTCGCGAGCAGTCCACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTc                                                                                   <  1:118580/245‑1 (MQ=255)
                                                                                                                                                                    ggTCACGCGGTATCCCGTATCGCCCAGACCGTAACCCTGCTGCCAAACCCCGGGGATAACAGGGGCATTAAAGCGCGCGAGCAGTCCCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACTGGTCCGTCAcc                                                                                <  2:59327/251‑1 (MQ=255)
                                                                                                                                                                     gTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCAcct                                                                               >  1:482487/1‑251 (MQ=255)
                                                                                                                                                                        ccgcggTATTCCGTATTGCCCTGCCCGTAACCCTGCTGCCAGACCCCGGGGATAGCAGAGGCATAAAAGCTCGCGAGCAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCAcctcct                                                                            <  2:430574/246‑1 (MQ=255)
                                                                                                                                                                         cgcgttatcccGTATTGCCCTGACCGTAACCCGGCTGCCAGACCCCCGGGATAGCAGAGACATAAAAGCCCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCACTCTTGTCGCCGTTGTCGTTTTTTTATTCCTGATTGGCAGGGTTCGAGGGTATCAGTAGCTTCCCCATCGGTTTGCCGTCAGGTGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTg                                                                           <  1:502801/242‑1 (MQ=255)
                                                                                                                                                                                             ccgaacagaacccgtgcggccacacccccGGGATAGCAGAGGAATAAAGGCACGCGAGCAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGtt                                                         <  2:160402/226‑1 (MQ=255)
                                                                                                                                                                                                         ctgctgACAGACCCCCGGGATAGCAGAGGCATTACAGCTCGCGAGCAGTACACCCGCCAGCATCAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTACTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTAATCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTACAGTTTGATGTTATCGGCAATacac                                            >  2:239116/1‑251 (MQ=255)
                                                                                                                                                                                                          cgtctgcCATACCCCCGGGAGAGGAGAGGCATTAAGGCTCGCTGGCTGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATacac                                            <  1:106367/248‑1 (MQ=255)
                                                                                                                                                                                                            gctgcCAGACCCCCGGGATAGCAGAGGCTTAAAGGCTCGAGAGGAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTGTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATacac                                            <  1:52151/249‑1 (MQ=255)
                                                                                                                                                                                                                  aGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGTAAGCATAAACCGACTTAGTGAAACTATCGTCTTTCTTGCCTCTTTTGTCGTTGTTTTATTACTGATTGACCGGGTTCGAGGGTATCAGTAGCTGTCCCATCAGTTTGCCGTCCTGGGCGTACTCAAAGTATTTTTCTTTGGTGTACGGGTCCGTCACCACCTGGTATTCCAGTTTGATGTTATAGGCAATACACAGTGCAtt                                    >  1:550723/1‑251 (MQ=255)
                                                                                                                                                                                                                              aTAGCAGAGGAATTAAGGCTCGGGAGCAGTTCACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCGGGGATTCCAGTGTGATGTTAGCGGCAATACACAGCGCATTCATCAGTGGCt                         <  1:73606/251‑1 (MQ=255)
                                                                                                                                                                                                                                 agcagaggcattaaccccggcggtcAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCGTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTgg                       <  2:215868/227‑1 (MQ=255)
                                                                                                                                                                                                                                   agagGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACg                    >  2:85602/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cATTAAAGCCCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGtttt               <  1:421756/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    gcgAGCAGTACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCatat    >  2:77126/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      gAGCAGTACACCCGCCAGCATAAAACGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATcc  >  1:126018/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       aGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCACTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTgactatcc  >  2:279/1‑243 (MQ=255)
                                                                                                                                                                                                                                                          aGTACACCCGCCAGCATACACCGACTTAGTGAAAGTATTGTCATTCCTGTATCCTTTGTCGTTGTTTTATTCCTGCTTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTAtt                                                                                                          >  2:503558/1‑143 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
AATGTACTGCTGGTGATCCATCACTACGGTGATAGTAGTGTCGTCATCGTGCGAACTGACCATCTTGTCATCGGCAAGGGTAAGAAAGACTGAATGCTGATAGAAACCATTCTGGTCCGGGTTCCCTGTGCAGTTGATGGTAAACGTCTTTCCGCTGGCTTCGGTCACGCTGTATTCCGTATTGCCCTGACCGTAACCCTGCTGCCAGAACCCCGGGATAGCAGAGGCATTAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCC  >  NC_000913/265800‑266290

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