Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,825,061	0	G	A	37.0%	77.8 / 17.0	46	intergenic (‑124/+79)	ves/spy	HutD family protein Ves/ATP‑independent periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base G (27/2);  new base A (0/17);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 9.77e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCG  >  NC_000913/1824816‑1825311                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                           gCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAAc                                                                                                                                                                                                                                                       <  2:146391/251‑1 (MQ=255)      ggtggaaCGCCAATAAAAATCCCGTTTTGCCGGTGGGAACGTGCAATTTTCGCGCGTTTCCCCGGCAGCATTCCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCTTTTTTTT‑CTTTGTTTTCTTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTGGGGTGATTTACGTAAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAAAA‑GAGAAAAAAAAAccgccg                                                                                                                                                                                                                                                  <  2:13902/245‑1 (MQ=255)          gCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATcc                                                                                                                                                                                                                                              >  1:419005/1‑251 (MQ=255)              gCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTc                                                                                                                                                                                                                                          >  1:586581/1‑251 (MQ=255)                 aGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                                                                                                                                       >  2:159638/1‑251 (MQ=255)                                ttGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACCTTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                                                                                                                                        >  2:268150/1‑251 (MQ=255)                                   ccGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                                                                                                                                                                                                       >  1:741085/1‑249 (MQ=255)                                                   aaTTTGGCGCGTTACCCCGGCAGCATTCCGCCAAAGATTCCCCCACATTTTCCGCATATAAAAGTATTCCATACTCGCTCCCGGCCTGTTTTTTT‑CATTTTTATCGTTCAGGTACCGGTGGATGGTT‑AGCGGCA‑ACTCGGGTGATTTCCGTTAGGGGTGATCGGGTTTTTTTT‑AC‑AACTGTTGGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg                                                                                                                                                                                                     <  1:530442/251‑1 (MQ=255)                                                          cgcgttccccgggccgcttttCCCCACAGTTTCCCCAACTTTTTCCGCATACCAAATTTTTCCATCCTCGCTCCCGCCCTTTTTTTTT‑TATTTTTATCGTTAAGGTAACGGTGGATGTTT‑ATCGGCA‑CCTGGGGTGATTTACGTTAGGGTTGATC‑GGGTTTTTTT‑ACCACCTTTTTGAAAA‑GAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAg                                                                                                                                                                                              <  1:167747/233‑1 (MQ=255)                                                              tttCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttt                                                                                                                                                                                          >  1:109093/1‑251 (MQ=255)                                                                                   acaGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTaagc                                                                                                                                                                     >  1:286047/1‑251 (MQ=255)                                                                                        tcccccGACATTTTACTAATATAAAATTTTTCCATCCTCTCTCCCGGCCCTTTTTTTT‑CTTTTTTATCTTTCAGGTACCGGTGGAGTGTT‑AT‑GGCATCCTGGGTTGTTTTGCGTTAGGGGTGATCAGGTTTTTTTT‑AC‑ATCTGTGTGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGAGTATTCAGCAGtt                                                                                                                                                                 <  1:562198/246‑1 (MQ=255)                                                                                         cccccaacattttcccaatatAAAAGTTTCCCATCCTCGC‑CCCCGCCCGTTTTTTTTCTTTTTTATCGTTAAGGTACCGGTGGTTGTTT‑ATCGGAA‑ACTGGGGTTATTTACTTTGGTGGTGATCAGTTTTTTTTT‑AC‑ATCTGTTGGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc                                                                                                                                                               <  1:192810/234‑1 (MQ=255)                                                                                                 aTTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAACATTATTTAGCCATTGCAGGCAttt                                                                                                                                                       >  2:169829/1‑251 (MQ=255)                                                                                                    ttacgtaaaaaaaattATTCCATCCCCGCCCCCGGCCCTTTCTTTT‑TTTTTTTTTCGTTCAGGTACCCGTGGTTGGTT‑ATCGGAA‑ACTGGGGTGATTAACGTTGGTGGTAATCAGGTTTTTTTT‑AA‑ACCTGTTTGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                                                                                                                                     <  2:109093/235‑1 (MQ=255)                                                                                                    ttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAc                                                                                                                                                    >  1:724767/1‑251 (MQ=255)                                                                                                      aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCt                                                                                                                                                  >  1:489352/1‑251 (MQ=255)                                                                                                          atatCAAAGTATTCCATACTAGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTTGCTGAAAGATTATTCAGCAGTTGTAGGCATTTTTCCttttt                                                                                                                                              >  2:661171/1‑250 (MQ=255)                                                                                                            ataaaagttttccttcctcccccccGGCCCGTTTTTTT‑CATTTTTATCGTTCAGTTACCCGTGG‑TGGTTTGTCGGCA‑ACTGGGGTGATTTACGTTAGTGGTGATC‑GGCTTTTTTTTAC‑AACTTTTGGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGcc                                                                                                                                            <  1:268150/229‑1 (MQ=255)                                                                                                               aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTGCCGCt                                                                                                                                         >  1:582395/1‑251 (MQ=255)                                                                                                                 aGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCCTTTAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGGATTTTACCTTTTGCCGCTgg                                                                                                                                       >  1:689461/1‑251 (MQ=255)                                                                                                                                     ggCCCTTTCTTTT‑CATTGTTATCGTTCAGGTACCCGTGGATGGTT‑ATCGGCA‑ACTGGGGTGATTTCCGTTGGTGGTGATCAGGTTTTTTTT‑AC‑AATTGTTTGAGAA‑AAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtct                                                                                                                   <  1:100162/251‑1 (MQ=255)                                                                                                                                     ggCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTTTCAGAAGCTtct                                                                                                                   >  2:100162/1‑251 (MQ=255)                                                                                                                                           ttttttt‑CATTGTTATCGTCGAGTTACCGGGGGTTGTTT‑ATCGACA‑ACTGGGGTTTTTTACGTTAGTGGTGATCGGGCTTTTTTT‑AC‑AAC‑GTTTGAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaa                                                                                                              <  1:169829/247‑1 (MQ=255)                                                                                                                                              tttt‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAACGAAAACCGACGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTTTTCAACAGTTGCAGGCCTTTTTCCTTTTGCCGCTGGACTTTTTGTCAGACGCTTCTCAAAAtt                                                                                                            >  1:259286/1‑249 (MQ=255)                                                                                                                                                        ttATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGTCGGTCTGTCAGACGCGTCTCAAAATTAGCATTTAAtt                                                                                                 >  1:256438/1‑251 (MQ=255)                                                                                                                                                         tATCGTTCAGGTACCGGTGGATGGTT‑ATCGGCA‑ACTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTT‑AC‑AACTGTTGGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAAtt                                                                                                 <  1:159638/250‑1 (MQ=255)                                                                                                                                                          atgtttCAGGTAGCGGTGGTTGTTT‑ATCGGCA‑ACTGTGGTGATTTCCTTTAGTGGTGAACGGTTTTTTTTT‑AC‑ACCTGTTTGAGAG‑AAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATGAGCATTAAATTg                                                                                                <  1:561619/246‑1 (MQ=255)                                                                                                                                                                 aGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCCTAAGTTTTTTTGGTGAAAGATTATTCAACAGTTGCAGGCATTTTACCTTTTGCCCCTGGTGGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAAtttgttttttt                                                                                        >  2:52628/1‑242 (MQ=255)                                                                                                                                                                  ggTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCtt                                                                                       >  2:186931/1‑250 (MQ=255)                                                                                                                                                                    ccgtatctgttgttt‑ATCGCCA‑CCTGGGGTGATTTACGTTAGGGGTGCTCCGGTTTTTTTT‑AC‑AATTGTTTGAGAA‑AAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCt                                                                                        <  2:717899/235‑1 (MQ=255)                                                                                                                                                                            gATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCt                                                                                        >  1:717899/1‑240 (MQ=255)                                                                                                                                                                             gttgtttaccGGCA‑ACTGGGGTGATTTACTTTATTGGTGATCAGGTTTTTTTT‑AC‑AACTTTTTGAAAA‑AAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTc                                                                             <  1:779643/242‑1 (MQ=255)                                                                                                                                                                             gattgttttagcaa‑aCTGGGGTGATTTACGTTAGTGGTGAGCAGGTTTTTTTT‑CC‑ACCTGTTTGAAAA‑GAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTc                                                                             <  2:489352/238‑1 (MQ=255)                                                                                                                                                                                gTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGAGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GCGAAAAGAAAACAGCCGATACTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGGCGTTCTGTCAGAAGCTTCTCAAAATTAACATTAAATTGCTTTTTCTGTTCCGGcttctggg                                                                         >  1:50181/1‑244 (MQ=255)                                                                                                                                                                                              gggggttatttagtttAGTGGTGATCAGTTTTTTTTT‑AC‑AACGTTTTGAGAA‑AAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAAttt                                                             <  2:582395/237‑1 (MQ=255)                                                                                                                                                                                               gggttgttttACGTTTGTGGTGACCAGTTTTTTTTT‑AC‑ATCTGTTGGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAAtttt                                                            <  2:724767/244‑1 (MQ=255)                                                                                                                                                                                                           gTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTc                                                                                                                                 >  1:185107/1‑170 (MQ=255)                                                                                                                                                                                                               gtggtgATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCAtg                                            >  1:94922/1‑251 (MQ=255)                                                                                                                                                                                                                       cagttttttttt‑AC‑AACTGTTGGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCa                                     <  2:94922/245‑1 (MQ=255)                                                                                                                                                                                                                            ttttttt‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTTTTTTGGGTTTCCATGTGCGCCAGCATGt                               >  2:128513/1‑251 (MQ=255)                                                                                                                                                                                                                               tttt‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCCTGTGCGCCAGCATGtt                              >  1:224051/1‑249 (MQ=255)                                                                                                                                                                                                                                       aaCTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTgc                      >  1:683421/1‑251 (MQ=255)                                                                                                                                                                                                                                           gTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTg                  >  1:259652/1‑251 (MQ=255)                                                                                                                                                                                                                                               gagaa‑gagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCCTTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGttct              >  1:623680/1‑251 (MQ=255)                                                                                                                                                                                                                                                aaaa‑gaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGttctt             <  2:23799/249‑1 (MQ=255)                                                                                                                                                                                                                                                     aaaaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCAt         <  1:430162/248‑1 (MQ=255)                                                                                                                                                                                                                                                       gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCAttt       >  1:588318/1‑251 (MQ=255)                                                                                                                                                                                                                                                        naaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGGTGTTGTAAATTTTGTTCTGGGTTTCCATGTTCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCAtttt      >  2:347111/2‑251 (MQ=255)                                                                                                                                                                                                                                                         aaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTACGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCCTTTTTCCTTTTTCCGCCGGACGATCTGTCAGACGCTTCTCCAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCCGCCGGATGTTGTAAATTTTGTTCTGGGTTTTCATGTTCGCCAGCCTGTTAGCTTTTCGCTTTTCTTCCCttttt     >  2:286405/1‑251 (MQ=255)                                                                                                                                                                                                                                                            gaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCGTTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGGTTTCATGTGCGCCAGCATGTTAGCTTTTCGCTGTTCTTCCATTTTTGCg  >  2:751031/1‑251 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                           GCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTT‑CATTGTTATCGTTCAGGTAACCGTGGATGGTT‑ATCGGCA‑ACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑AC‑AACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCG  >  NC_000913/1824816‑1825311

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