Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,001,911	0	G	A	33.3%	52.6 / 14.2	36	intergenic (‑122/+1)	fliA/fliX	RNA polymerase sigma factor FliA/small RNA FliX
Reads supporting (aligned to +/- strand):  ref base G (20/4);  new base A (0/12);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.45e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCT  >  NC_000913/2001667‑2002155                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       aaGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGaaaaa                                                                                                                                                                                                                                                   >  2:487039/1‑250 (MQ=255)         cGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCcc‑gcc                                                                                                                                                                                                                                        <  2:117621/251‑1 (MQ=255)          gCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCcc‑gccg                                                                                                                                                                                                                                       >  1:468442/1‑251 (MQ=255)                gggCAGTGGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCg                                                                                                                                                                                                                                 >  2:200404/1‑251 (MQ=255)                       cgacccTGCAGGCGCAATGCTTCGTGACGCCCCAGGGGGAAAAAACGCTGCCACAGCGATTGTTTATCCTTTCCCCCTTCAGCGGTATAGAGTGAATTCACGAAAAACGGCCCTGCTTTATATGAGTTACCGGCATGTTTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGG‑TAAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATCCC‑AGGCAAAAAAA‑ACCCCC‑GCCGGTGGCGGGGAAGc                                                                                                                                                                                                                          <  1:11063/247‑1 (MQ=255)                          aCCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACg                                                                                                                                                                                                                       >  1:273500/1‑251 (MQ=255)                             tGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTg                                                                                                                                                                                                                    >  1:165353/1‑251 (MQ=255)                                     cAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGAACGTTGCTGACaa                                                                                                                                                                                                             >  2:272566/1‑250 (MQ=255)                                                    cACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                                                                                                             >  1:119742/1‑251 (MQ=255)                                                                      cTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                                                                                                           >  1:329356/1‑251 (MQ=255)                                                                                            aTTACCCCTTCAGGGGTATAGAGTAAATTCACAAAAAACACCCCTGCTTTATATGAGTTATCGGAATGTTTATCCTTTTCTGCAGGTTTTTTAATCGGACGATTAGGGGG‑GAAAATG‑GGGGTTTTTTGGGGGGTT‑A‑CAGGTAAATTCC‑AGGCAAAAAAA‑ACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                                                                                                                     <  1:487039/251‑1 (MQ=255)                                                                                            aTTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                                                                                                                     >  1:354022/1‑251 (MQ=255)                                                                                                          ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATACC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                                                                                                                        >  2:12265/1‑250 (MQ=255)                                                                                                              taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                                                                                                                   >  1:523103/1‑251 (MQ=255)                                                                                                                    gaatccacaaaaaaCACCCCTGGTTTATATGAGTTATCGGCATGATTTTCCGTTTCTGCAGGTTTTTAAAGCGGACGATTAGTGGG‑TGAAAGGGGGGTTTATTGGGGGT‑TT‑A‑CGGGTAAATTCC‑AGGCAAAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                                                                                                              <  1:285341/239‑1 (MQ=255)                                                                                                                     aaTTCACAATAAACAGCCCGGGTTTATATGATTTAGCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGGGGG‑GAAAAGGGGGGGTTTTTTGGGGG‑TG‑A‑AAGGTAAATCCA‑AGGAAAAAAAA‑ACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                                                                                                            <  2:68372/251‑1 (MQ=255)                                                                                                                           cGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                                                                                                                      >  2:405304/1‑251 (MQ=255)                                                                                                                           agataaaacgCCCTGCGTTATATGAGTTATCGGCTTGATTATCCGTTTCTGCGGGGTTTTAAATGGGACGATTAGTGGT‑TGAAAGGAGGGGTTTTTTGGGGG‑TT‑A‑AGGGTAAATGCC‑AGGAAAAAAAA‑ACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                                                                                                                      <  2:448394/242‑1 (MQ=255)                                                                                                                                            ttATATGAGTTATGGGCATGATTATCCTTTTCGGCAGGTTTTTTAATCGGACGATTAGTGGG‑TAAAAGGAGGG‑TTTTTTGGGGG‑TTTA‑CAGAAAAATCCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCGAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca                                                                                                      <  1:546553/250‑1 (MQ=37)                                                                                                                                             aataggatttATCGACTTGATTATCCGTTTCTGCAGGTTTTTTAAGCGGAGGATAGGGGGG‑TGAAAGGAGGGTTTTTTTGGGGG‑TT‑A‑CAG‑AAAAATCC‑AAGAAAAAAAA‑AACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagc                                                                                                    <  2:196755/243‑1 (MQ=25)                                                                                                                                                  gAGTTATGGGCAGTATTATCCTTTTCGGCAGGGTTTTAAATCGGACGATTAGTGGG‑AGAAAGGAGGGGTTATTGGGGGG‑TT‑A‑CAGGTAAATCCC‑AGGAAAAAAAA‑ACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc                                                                                               <  1:139393/251‑1 (MQ=37)                                                                                                                                                   aGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc                                                                                               >  1:191535/1‑250 (MQ=255)                                                                                                                                                      tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGGAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                           >  1:81786/1‑251 (MQ=255)                                                                                                                                                       agcgacttGATTATCCGTTTTTGCGGGTTTTTTAAT‑GGAGTATTAGTGGGGTGAAAGGAGGGGTTTTTTGGGGT‑TG‑A‑AAGGAAAATTCA‑GGGCAAAAAAA‑ACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                          <  2:87129/244‑1 (MQ=25)                                                                                                                                                                  tATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTg                                                                               >  2:237051/1‑251 (MQ=255)                                                                                                                                                                                     tttAATGGGACGATTAGTGGG‑TAAAAGGGGGGTTTTTTTGGGGG‑TT‑A‑AAGGTAAATTCC‑AGGAAAAAAAA‑ACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAAt                                                             <  2:510504/250‑1 (MQ=37)                                                                                                                                                                                      taaagggaacaaTTAGTGGG‑TGAAAGGGGGGGTT‑GTTGGGGG‑TT‑A‑ACGGTAAAATCCGAGGCAGAAAAA‑ACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATc                                                            <  1:185247/239‑1 (MQ=37)                                                                                                                                                                                                        gg‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGAAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTAATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAAccc                                         >  2:27243/1‑251 (MQ=255)                                                                                                                                                                                                                  ggggtttATTTGGGGG‑TT‑A‑AGGGAAAATTCC‑AGGAAGAAAAA‑CCCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCaga                                <  2:154050/246‑1 (MQ=38)                                                                                                                                                                                                                     ggTTTTTTGGGGG‑TT‑A‑AAGGTAA‑TTCC‑CGGGAAAAAAACAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGAc                             <  2:398953/251‑1 (MQ=38)                                                                                                                                                                                                                     ggTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGAc                             >  1:263070/1‑251 (MQ=255)                                                                                                                                                                                                                      gttttttGGGGG‑TT‑AACAG‑TAAATTCC‑AGGCAAAAAAA‑AACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGAc                             <  2:80250/247‑1 (MQ=255)                                                                                                                                                                                                                             tgggg‑TT‑AACAGTTAAATCCC‑AGGCAAAAAAA‑ACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAAcc                       <  2:273500/249‑1 (MQ=255)                                                                                                                                                                                                                                gg‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACctgctg                  >  1:428827/1‑251 (MQ=255)                                                                                                                                                                                                                                                   c‑AGGCAGAAAAA‑AACCCC‑GCCGCTGGCGGGGAAGCACGTTGCTGACAAACCACGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCCAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACCACGACTAACCATGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTCAGCt  >  2:192735/1‑251 (MQ=38)                                                                                                                                                                                                                                                       gCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGAAAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                                                                                                                    >  1:508947/1‑102 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       AAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGG‑TGAAATGAGGGGTTATTTGGGGG‑TT‑A‑CAGGTAAATTCC‑AGGCAGAAAAA‑AACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCT  >  NC_000913/2001667‑2002155

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