Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,993,470	A→T	9.0%	*459Y (TAA→TAT)	ygeH →	putative transcriptional regulator YgeH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,993,470	0	A	T	9.0%	124.1 / 3.0	46	*459Y (TAA→TAT)	ygeH	putative transcriptional regulator YgeH
Reads supporting (aligned to +/- strand):  ref base A (22/19);  new base T (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

GCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTAAAATGCTCTAGAGTTAAAAGTATGATGTAAAGACCATTGGCCTGGACAATTAAAAAAGGGAGCTTTATGACTAATCCAATCGGTATCAATAATTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCAT  >  NC_000913/2993207‑2993716                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       gCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTGGTGATAAAAGAATGTGTAGATATGTATGTGCCTACCCCGCTCACAAACAACATTAAACTCTACTATCAAGAAACTGAGAGTGAAAGCCAACGAGTTATAATTGACAACATTTTGAATTTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGa                                                                                                                                                                                                                                             >  1:203093/1‑275 (MQ=255)           caaactaTAGGGACCTAGTCCACTTTACTAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGGAAGCCATCGAGTTATAATTGTCAACATTTTGAAATTAAAGCAGCTGACGATAATTTGTATGCGATAAATTTTTATTGACGCAATTtta                                                                                                                                                                                                                                   <  2:203092/272‑1 (MQ=255)           caaaaaataGGGAACTAGTCAACCTTCATAATGAAAAAATTGAAGAGGCAAGGAGTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACCATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATCGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTtta                                                                                                                                                                                                                                   <  2:123933/269‑1 (MQ=255)           actactaTAGGGCACTAGTCCACTTTAATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGGATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTtta                                                                                                                                                                                                                                   <  1:219920/275‑1 (MQ=255)                       cACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTACACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAAt                                                                                                                                                                                                                                       >  1:123933/1‑259 (MQ=255)                                                                      ataGACAAATCACTACAACTCGAACCCAGAAGAAGAAAAGCAGTTGTGATAAAAGAAGGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCAAGTTAGAATTGACAACAATTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTAGTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATAtt                                                                                                                                                                        <  1:349397/275‑1 (MQ=255)                                                                               tCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAAAAACATTAACCTCTACTATAAAGAACCTGAGAGTGAAAGCCATCGAGTTATAATTGACAAAATTTTGAAATTAAAGCAGCTGACGAGAATTGGTATGCGATAAATTTTTATGTACGCAATTTTATTCTTTAGATGTATGCAACTTg                                                                                                                                                                                                             >  1:139161/1‑229 (MQ=255)                                                                                        cTCGAACCCAGAAGACGAACAGCAGTTGTGATATAAGAATGTGTAGATATGCATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAACGCCATCGAGTCATAATTGACAACCTTTTGAAATTATAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGCCGCAATTTTATTATTTAGATGTATGCATCTTGTTATTTata                                                                                                                                                                                                    >  1:77349/1‑227 (MQ=255)                                                                                            aaCCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGGTATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTAGAAAGAGACGGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAGATTTTTAGTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAaacacc                                                                                                                                                                                                <  2:77349/229‑3 (MQ=255)                                                                                                   aagaCGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATc                                                                                                                                                                                 >  1:63149/1‑237 (MQ=255)                                                                                                      aCGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAAt                                                                                                                                                                                                                                       >  1:131038/1‑180 (MQ=255)                                                                                                       cGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAg                                                                                                                                                                             <  2:63149/237‑1 (MQ=255)                                                                                                         aaaaGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTAAAATGCTCTAGAGt                                                                                                                                     >  1:316057/1‑275 (MQ=255)                                                                                                          aaaGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTtta                                                                                                                                                                                                                                   <  2:131037/180‑1 (MQ=255)                                                                                                                                             tgtgCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGCTGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTAAAATTCTCTAGAGTTAAAAGTAATATTTAAAAACCAtttttttttttttt                                                                                                 >  2:401237/1‑263 (MQ=255)                                                                                                                                                gCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAata                                                                                                                                                                                   >  1:11453/1‑190 (MQ=255)                                                                                                                                                 ccTAACCCGCTCAAAAACAACATTAAACCCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTAAAATGCTCTAGAGTTAAAAGTATGATGt                                                                                                                       <  1:401238/249‑1 (MQ=255)                                                                                                                                                     aCCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTAAAATGCTCTAGAGTTAAAAGTATGATGTAAAGACCATTGGCCTGGACAATTaaaaa                                                                                           >  1:104658/1‑273 (MQ=255)                                                                                                                                                     aCCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTAAATATGCTCTAGAGTTAAAAGTATGATGTAAAg                                                                                                                   >  1:387296/1‑249 (MQ=255)                                                                                                                                                       ccacaaaaaaaaaaaaTTAAACTAGACAATAAATTAGCTTAGAGTGAAAGACTTCGAGTTATAATTGAAAACATTTTGAAATGAAAGGTGCTGACGAGAATTTTGATGCGATAAATTTGTAATGACGAAATTTTATTATTTAGGTGTAAGCAACCTGTTTTGTAAACAACAGAGATCAAAATATCTAAGATATTCACCTTAAGGAAAGATGTTTAAAATGCTCTAGAGTTAAAAGTATGTTGTAAACACCATTGGCCTGGACAATTAAAAAAggg                                                                                       <  2:104658/269‑1 (MQ=255)                                                                                                                                                          cccAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGAGATAAATTTTTATTGACGCAATTTtattat                                                                                                                                                                                                                               >  2:192807/3‑136 (MQ=255)                                                                                                                                                           tCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATATATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAata                                                                                                                                                                                   >  2:330042/1‑179 (MQ=255)                                                                                                                                                           tCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAata                                                                                                                                                                                   >  1:306874/1‑179 (MQ=255)                                                                                                                                                           tCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATAtt                                                                                                                                                     >  2:25233/1‑209 (MQ=255)                                                                                                                                                           tCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCAGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAAt                                                                                                                                                                                                                                       >  2:249841/1‑127 (MQ=255)                                                                                                                                                              aaaaCAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAg                                                                                                                                                                                                                           <  1:192808/136‑1 (MQ=255)                                                                                                                                                               aaaCAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATATATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTa                                                                                                                                                                               <  1:330043/179‑1 (MQ=255)                                                                                                                                                               aaaCAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTaaa                                                                                                                                                 <  1:25233/209‑1 (MQ=255)                                                                                                                                                               aaaCAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCGGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTtta                                                                                                                                                                                                                                   <  1:249842/127‑1 (MQ=255)                                                                                                                                                                               tactaTAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTAAAATGCTCTAGAGTTAAAAGTATGATGTAAAGACCATTGGCCTGGACAATTaaaaa                                                                                           >  1:174529/1‑247 (MQ=255)                                                                                                                                                                               tactaTAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGAAACTTGTTATTTAAACCACAGAGATAAAAATCTCTAAGATATTCACCTTATTGCAAGATATTTAAAATGCACTAGAGTTAAAAGTAAGATGTAAAg                                                                                                                   >  1:175416/1‑223 (MQ=255)                                                                                                                                                                                   ataAAGAAACCGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTAAAATGCTCTATAGTTAAAAGTATGATGTAAAGACCATTGGCCTGGACAATTAAAAAAggg                                                                                       <  2:174528/247‑1 (MQ=255)                                                                                                                                                                                   aaaaaGAAACTGAGCGGTAAAGCCATCGAGTTATAATTGACAACAGTTTTAAAGTAAAGCAGCTGATGCGAATTTGTATAGGATGAATTTTTATTGACGCAATTCTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATGAAAACATCTAAGATATTCACCTGTTTGCAAGTTATTTAAACTGCTCTAGAGTTAAAAGTATGATGTAAAGACCa                                                                                                               <  2:175415/221‑1 (MQ=255)                                                                                                                                                                                         aaaCTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAata                                                                                                                                                                                   >  1:379502/1‑149 (MQ=255)                                                                                                                                                                                             tGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAAGTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTa                                                                                                                                                                               <  2:379501/149‑1 (MQ=255)                                                                                                                                                                                                                          acaacaTTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATATATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTAAAATGCTCTAGAGTTAAAAGTATGATGTAAAGACCATTGGCCTGGACAATTAAAAAAGGGAGCttc                                                                                 >  1:363931/1‑213 (MQ=255)                                                                                                                                                                                                                          acaacaTTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAACAATATCTAAGa                                                                                                                                                                            >  1:244004/1‑123 (MQ=255)                                                                                                                                                                                                                              caTTTTGAAATTAAGGCAGCTGACGATAAGTTGTTTGCGATAAATTTTTTTTGACGTAATTTTATTAATTAGATGTATTCAACTTGTTATTAGAACGACAGAGATAAAAATATCTAAGATATTCATCTTATTGCAAGAGATTGAAAATGCTCTAGAGTTAAAAGTATGATGTAAAGACCATTGGCCTGGACAATTAAAAAAGGGAGCtt                                                                                  <  1:386012/209‑1 (MQ=255)                                                                                                                                                                                                                              caTTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATATATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTAAAATGCTCTAGAGTTAAAAGTATGATGTAAAGACCATTGGCCTGGACAATTAAAAAAGGGAGCTTCATGa                                                                             <  2:363930/214‑1 (MQ=255)                                                                                                                                                                                                                              caTTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATAtt                                                                                                                                                                        <  2:244003/123‑1 (MQ=255)                                                                                                                                                                                                                                            agaagttaaCGAAAATTTGTATGCGATAAATTTTTATTTACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTAAAATGCTCTAGAGTTAAAAGTATGATGTAAAGACCATTGGCCTGGACAATTAAAAAAGGGAGCTTTATGACTAATCCAATCGGTATCAATAATTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCAt  <  2:316056/267‑1 (MQ=255)                                                                                                                                                                                                                                                 cTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTa                                                                                                                                                                                                          >  1:10265/1‑70 (MQ=255)                                                                                                                                                                                                                                                     cGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTa                                                                                                                                                                                                      <  2:10265/70‑1 (MQ=255)                                                                                                                                                                                                                                                           tttGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACACCAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTAAAATGCTCTAGAGTTAAAAGTATGATGTAAAGACCATTGGCCTGGACCATTAAAAAAGGGAGCTTTATGACTAATCCa                                                                     >  1:82121/1‑193 (MQ=255)                                                                                                                                                                                                                                                               tATTCGATAAAGTTTTATTGACGCAATTTTATTATTTAGGTGTATGCAACCTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATGGCAAGATATTTAAAATGCTCTAGAGTTAAAAGTATGATGTAAAGACCATTGGCCTGGACAATTAAAAAAGGGAGCTTTATGACTAATCCAATCg                                                                 <  2:82121/193‑1 (MQ=255)                                                                                                                                                                                                                                                                   cGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGAt                                                                                                                                                                                          >  2:6421/1‑68 (MQ=255)                                                                                                                                                                                                                                                                       aaaTTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATaaaa                                                                                                                                                                                      <  1:6421/68‑1 (MQ=255)                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       GCCTCTCTCTACTACTATAGGGCACTAGTCCACTTTCATAATGAAAAAATTGAAGAGGCAAGGATTTGTATAGACAAATCACTACAACTCGAACCCAGAAGACGAAAAGCAGTTGTGATAAAAGAATGTGTAGATATGTATGTGCCTAACCCGCTCAAAAACAACATTAAACTCTACTATAAAGAAACTGAGAGTGAAAGCCATCGAGTTATAATTGACAACATTTTGAAATTAAAGCAGCTGACGAGAATTTGTATGCGATAAATTTTTATTGACGCAATTTTATTATTTAGATGTATGCAACTTGTTATTTAAACAACAGAGATAAAAATATCTAAGATATTCACCTTATTGCAAGATATTTAAAATGCTCTAGAGTTAAAAGTATGATGTAAAGACCATTGGCCTGGACAATTAAAAAAGGGAGCTTTATGACTAATCCAATCGGTATCAATAATTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCAT  >  NC_000913/2993207‑2993716

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