Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,875,521	0	T	C	25.8%	89.7 / 16.0	63	A28A (GCT→GCC)	yidX	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (20/26);  new base C (15/1);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 3.82e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAG  >  NC_000913/3875236‑3875816                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                        aGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGCCTGTATCTCGTATGCTCTGGTTTCCCATAc                                                                                                                                                                                                                                                                                           >  1:318162/1‑300 (MQ=255)   ctagatgaAGTTCGCGTTTGCAGTGCGTGTAATGTAGTACAAACGTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCAGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATATCGTATGCTCTGGTTTCCCATAc                                                                                                                                                                                                                                                                                           <  1:299121/295‑1 (MQ=255)        tgaCGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAgg                                                                                                                                                                                                                                                                                   >  1:1927/1‑301 (MQ=255)                        cgagcgtgTAATGTAGGACAAACTTATATTGTTGTACTACAAATTAGATCACAACAAGAACAATGCAGTAAAAAATGACATACGTCGGGCAGAAATCTGAAAAAGGCTATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCGTTCAACGGCTTAGCTACTCGTGTAAAGGATTAATCATGAAGTTGGATGTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCGCTGGTTTCCCATGCCACAAAGGGTAGTTCAGGAAAGTa                                                                                                                                                                                                                                                                   <  1:613456/298‑1 (MQ=255)                         gtgcgtgTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTAt                                                                                                                                                                                                                                                                  >  1:278851/1‑301 (MQ=255)                          ggcgagtgatgtaggattgaattATATTGTTGTTCTCCTACTTATTTCACGAAAAGAAAAATGAATAAAAAAATATATTGCGGAAAGAAAAATTATGAAAAAAAAAATAAGGCGCAAAAAAGAATGAAAAGAAGCGTATACAATCAACAAAGTAGCTACTCGTTTAAAGAATTAATCATGAAGTTGAATTTTAAGCAATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTAt                                                                                                                                                                                                                                                                  <  2:117402/279‑1 (MQ=255)                               gtAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTTTGCTTTACCGAAAAGGGTAGTTCAGGAAAGTATCAATCg                                                                                                                                                                                                                                                            >  2:631809/1‑301 (MQ=255)                                tAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTCTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCg                                                                                                                                                                                                                                                            >  2:578726/1‑300 (MQ=255)                                 aaTGTAGTACAAACTTATATTATTGTACTACAATTTAGGTCACAAAAAGAAGAATGCATAAAAAAATGACATGCGTCATGCAGAAATCTCAAAAGAGATATCAAACGCTAAAAATGAGATAAAGAAGAGTATGCTTTCAAAGGATTAGCTACTCGTTTAAAGGATTAATAATGAGGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGCTTTATTTCCACCGGCGTTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAAGCGCa                                                                                                                                                                                                                                                          <  2:127458/301‑1 (MQ=255)                                                          tagtacagattaggacactAAAAGAACAGTGAATAAAAGAATGATTTGCGGCGGGCTGAAATCTGAAAAGAGATAATAGGCGCTAAACAGAAGGGAAAGAAGAAAATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTAGGATTTTAAGGGATTTTTTAAGGCGGCCGGTTTATTCTCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCg                                                                                                                                                                                                                                 <  2:318165/286‑1 (MQ=255)                                                               ttcagacgtttgatcatccaatgttGCATAAAAAAATGACATTCGTCGAGGAGAAATCTGAAAAGGGATATCAGGAGCTAGACGGGAGAGAAAGAAGAGTATGCGTTCAAAGGCATAAATACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGAATTTTTTAAGGCTGCCGGTTTAATCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAg                                                                                                                                                                                                                              <  2:315542/275‑1 (MQ=255)                                                                tggagttcagacgtgtgctcttccgatctaaaaaaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCa                                                                                                                                                                                                                           <  2:216594/272‑1 (MQ=255)                                                                aaTTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAggaa                                                                                                                                                                                                                           >  1:386389/1‑299 (MQ=255)                                                                      gATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTCCCATACCACAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAgtt                                                                                                                                                                                                                     >  2:28099/1‑299 (MQ=255)                                                                              aaaGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCGTTTGCCTTACCACAAAGGGTAGTTCAGGAAAGTATCAATCGCATTCAGACACCATCACTTGCCTATCCCAGGCAAACGTTAGTACTgg                                                                                                                                                                                                             >  2:648962/1‑301 (MQ=255)                                                                                 tccaatctcATAAAAAAAGGACATGCGTCGGGCAGAAAGCTGAAAAGGTATATCAGGAGATAAACAGGAGGGAAAGAAGAGTATTCTTTCAACGGCATAGGTACTCGTTAAAAGGATTAATCATGAAGTTGAATTGTAAGGGATGTTTGAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTAGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAAc                                                                                                                                                                                                          <  1:526938/293‑1 (MQ=255)                                                                                  aagaaTGAATACAAAAATGACATGCTTCGGGAAGAAATCTGAAAAGGGATGTCAGGCGATAAACGGGAGAGAAAGAAGAGTATGCTGTCAACGGGTTAGATACTCGGTTAAAGGATTAATCATGAAGTTAAATTTTAAGGGATTTTTTAGGGCTGCCGGGTTATTCCCACTGGCGATAGTGCTTTCAGGCTGTATCGCGTGTGCTCTGGTTTCCCATACCGCAAGGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCGGGCAAAGGATAGTAATGGGACa                                                                                                                                                                                                         <  1:502933/298‑1 (MQ=255)                                                                                       tgaaaaaaaaaaTGTCAGGGATCGGACATAAGTCTGTAAATGGATATAAGTCTCTAAACAAGATAGAAAAAAAAGTTTGTTTTCAACGGCTTATCTACTCGTTAAAAGAAATAATCATGAAATTGAATTTTAAGTGATTTTTTAAGGCTACCGGTTTATTCCCGCTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAgg                                                                                                                                                                                                    <  2:278854/296‑1 (MQ=255)                                                                                       tGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAgagaccggaagagcgccgtgtagggaa                                                                                                                                                                                                    >  1:315539/1‑275 (MQ=255)                                                                                        tgtgctcttccgatcttGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGc                                                                                                                                                                                                   <  2:79441/285‑1 (MQ=255)                                                                                         cATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGTTTTCCCATAACGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCACCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACagatcgga                                                                                                                                                                                                  >  2:526942/1‑294 (MQ=255)                                                                                          aTAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTCGTTTGCCATAACGCAAAGGTTAGTTCAGGAAAGTATCAATCGCAGTCAGACACTATCACTTGGCTATCGCAGGCAAAAGATCGATATGGAACACAAGGCTa                                                                                                                                                                                                 >  2:327639/1‑301 (MQ=255)                                                                                             aaaaaaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGTTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGt                                                                                                                                                                                              >  2:276097/1‑301 (MQ=255)                                                                                             aaaaaaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCaagatcggaagagcgtcgtgtagggaaaga                                                                                                                                                                                              >  1:216591/1‑273 (MQ=255)                                                                                                        tGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGcagatcggaagagcgt                                                                                                                                                                                    >  1:79441/1‑285 (MQ=255)                                                                                                                       tCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGTTTTCCCTTACCGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACACAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTAc                                                                                                                                                                     >  2:153229/1‑300 (MQ=255)                                                                                                                          gAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTCTTCTGCTTTAAAAAAAAAATGAGTTCAGGAAAGTATCAATTGCAGTCAGACACCCTCACTGGGCTATAGCAGGCCAAAGATCGTAATGGAAAAAAAGGCTATGCTTTTGTAGGGGAATCGATAGAtcggaaga                                                                                                                                                                 >  2:340300/1‑294 (MQ=255)                                                                                                                          gAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCTTGACCGAAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCCAAAGATAGTAATGGAACAAAAGGCTATGTCTTTGTAAGGGAATCGCTGGATTACCttt                                                                                                                                                                 >  2:63504/1‑300 (MQ=255)                                                                                                                                tgttagcAGGCGCTAAACAGGAGAGAAAGAAGAGTATGATATAAACGGCTTAGCTACTCGTTTAAAGGAGTAATCATGAAGTTTAATTCTAAGGGATTTTTTAAGGCTTCCGGTTTATTCCCACTGGCGGTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCAATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAAGCGCTGGATTACCTTATCAc                                                                                                                                                             <  1:578720/293‑1 (MQ=255)                                                                                                                                 gATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAAGCTGCCGGTTTATTCCCACTGGCGCTTATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCTTAACGCAAAAGGTAGGTCAGGAAAAGTTCAATCCCTGTCCGGCCCCCTCACTTGGCCATTCCAAGCCAAAGATTAGAATTGAACCAAAAGCTATGTTTTTGTAAGGGCATTGGTTGATTTCCTTAttcattt                                                                                                                                                          >  2:465571/1‑295 (MQ=255)                                                                                                                                        ggCGCTAAACAGAAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATAACGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGTACAAAAGGATATGTTTTTGTAGGCGAATCGCTGGATTACCTTATCACTGATGGttc                                                                                                                                                    >  2:54383/1‑298 (MQ=255)                                                                                                                                           gcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCCCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGTGCAATCGCAGGCAAAAGATAGTAATGGAACAAACGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCCCTTAGTGTGcctag                                                                                                                                                >  1:495350/1‑298 (MQ=255)                                                                                                                                            cTAAACAGGAGGGAAAGAAGAGTATGCGTTCAACGGCTTAGATACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGAt                                                                                                                                                <  1:276094/300‑1 (MQ=255)                                                                                                                                              aaaCAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGCCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTTCCGATGa                                                                                                                                              >  2:362971/1‑300 (MQ=255)                                                                                                                                              aaaCAGGAGGGAAAGAAGAGTATGCATTCAACGTGTTATCTACTAGTTTAAAGGATTAATCATTAAGTTGAATTTTAAGTGATATTTTAAGGCTGCCGGTTTAGTCCCAAAGGAGCTGATGCTTTCAGGCTGTATAGCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATAGCAGTCAGACACCATCACTGTGCTATCGAAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGAc                                                                                                                                             <  1:631799/301‑1 (MQ=255)                                                                                                                                                 cAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTTACCACAAAGGGTAGTTCATGAAAGTATCAATCGCAGTAAGACTCCATCACTGGGCTATCGCTGGCAAAAGATAGTACTGGAACAAAAGGCTATGTTTTCGTATGGGAAGGGCTGGTTTACCTTCTCACTGATGGTTCCTATGAGAt                                                                                                                                           >  2:681959/1‑300 (MQ=255)                                                                                                                                                      gggAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTCTGCCAGACCGCAAAAGGTAGTTCAGGAAAGTATCATTCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGACTATGTTTTTGTAGGGGACTCGCTGGATTACCTTATCACTGATGGTGCCGAAGACATCGGTaa                                                                                                                                     >  2:423878/1‑301 (MQ=255)                                                                                                                                                              agagTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCa                                                                                                                             >  2:135720/1‑301 (MQ=255)                                                                                                                                                                agTATGATGTCTAGGGCTTAGCTTATCAGTTAAAAGATTAAGCATGAAGTTGAATTTTAAGTGAGTATTGAGGGCTGAGGGTTTAATCCCACTTGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAAGAAAGTATCAATCGCAGTTAAAAACCATCACTGAGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCGCTGATGGTGCCGATGACATCGTTAAGATGCTCa                                                                                                                             <  1:362968/299‑1 (MQ=255)                                                                                                                                                                     gCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATc                                                                                                                       >  1:378941/1‑300 (MQ=255)                                                                                                                                                                      cTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTCTTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGTTCTGCCATAAAAAAAAGGGTAGTTAAGTCAACTATAACTCACAGTCCGCTACGCTCCCTGGGCTTTCGCCGGCAAAATCTTGCAATGGAACACAAGGCTATTTCCTTGTAGGGGAATCGCTGCATTACCTTATCACTGCTGCCTCCAATTACATCCTTCTGAGCCTCAATGCTctt                                                                                                                     >  2:172480/1‑299 (MQ=255)                                                                                                                                                                       tttCAAGGGCTTAGGTACTAGTGTAAATGATTAATCATGAAGGTGAATTTTAAGGGATTGTTAAAGGCTGCCGATTTATTCCCACTGGGGCTGAGGCATACAGGAGGTATCGCGTATGATCTGGTTGCCCATGCCGCGAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGAAACCATCACTGGGCTATAGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAg                                                                                                                    <  1:648952/301‑1 (MQ=255)                                                                                                                                                                                        cTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGCATCGCTGGATTACCTTATCACTGATTGTTCCGATGACATCGTTAAGATGCTCAATGATCCAGCCCTTAACCGGCACAat                                                                                                   >  1:302733/1‑301 (MQ=255)                                                                                                                                                                                        cTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAat                                                                                                   >  2:454923/1‑301 (MQ=255)                                                                                                                                                                                           gTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACACAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGACTATGACATAGTTAAGATGACCAATGATACAGCACTTAAACGGccaactatt                                                                                                >  2:298902/1‑293 (MQ=255)                                                                                                                                                                                           gTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTCTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCCCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTGTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGCTGACATCGTTTAGATGCTCAATGATCCATCACTTAACCGGCACAATAtt                                                                                                >  2:234866/1‑301 (MQ=255)                                                                                                                                                                                            tttAAAGGATTAATCATGAAGTTGAATGTTAAGGTGTTTTTTAAGGCTGCCGGTTTATTCCCACGGGCGCTGATGCTTTCAGTCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAAGCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTc                                                                                               <  2:386392/301‑1 (MQ=255)                                                                                                                                                                                               aagtggttAATCATGAAGTTGAATTTTAAGGTATTTTTTAAGGCTGCCGGTTTATTCCCACTGGGGCTGATGCTTTCAGGCTGTATCTCGTATGTTCTGATTTCCCATACCGCAAAGAGTAGTTAAGGAAAGGATCCATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAATAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAGGATTCTCAATGTTCCAGCACTTAAACAGCACAGTATTCAgg                                                                                            <  1:34172/295‑1 (MQ=255)                                                                                                                                                                                                   gATTAATCATGAAGTCGAATTTTAAGTGATTTTTTAAGTCTGCCGGTTTATTCCCACTAACGCTGATGCTTTCAGGCTGTATCTCGTATGCTTTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTg                                                                                         <  2:1927/300‑1 (MQ=255)                                                                                                                                                                                                     ttAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCCAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCg                                                                                      >  1:192855/1‑301 (MQ=255)                                                                                                                                                                                                                  tGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTTGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGTCTATCGCAGGCAAAAGCTAGTAATGGAACACAAGGCTATGTTTTTGTAGGGGAATCGCTGGCTTACCTTATCACTGCTGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCCCCTAATTCAGGTTGCCGAGGACGCTAGAtt                                                                         >  1:257050/1‑301 (MQ=255)                                                                                                                                                                                                                      tgttAATGGATTCTTGAAGGCATCAGGTGTGTTACGACTGGCGCTGATACTTTGATGCTGTAACTCTAATGATCTGGTTTCACATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAATCAGAAACCAGCAAAGAGCTATCGCAGGCAAAAGATAGTAATGGGACAAAAGGCTATGTTGTTGTAGGGGAATCGCTGGATAAACTTATCACGGATGGTGCCGATGACATCGTTAAGATGCTAAATGATCCAGCACTTAACCGTCACAATATTCAGGTTGCCGATGACGTAAGATTTGtt                                                                     <  1:580637/299‑1 (MQ=255)                                                                                                                                                                                                                             ggATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTaaa                                                                <  2:378944/299‑1 (MQ=255)                                                                                                                                                                                                                                 tttttAAGGCTGCCGGTTTATTCACACTGGCGCTGATGCTTTCAGGCTGGATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCggg                                                          <  2:192858/301‑1 (MQ=255)                                                                                                                                                                                                                                                   tATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCgg                                        <  1:135718/301‑1 (MQ=255)                                                                                                                                                                                                                                                    aTTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTCTTTTTCCAGAAAACAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACTCTCTCACTCTGCTATCCCAGGCAAAAGACAGTAATGGAACAATAGGCTATGCTTTACTAGGGGAATCGCTGGATTACCTTATCACTGATGTTGCTGATGACATCGTTAAGATGCTAAATGATACAACACTTAACCGGCACAATATTCAGGTTGACGATGACGCAAGATTTGTTTTAAATGCGGGGCAAAAGAAATGTAACgca                                       >  2:149365/1‑299 (MQ=255)                                                                                                                                                                                                                                                       cccACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGcaca                                    >  2:521382/1‑301 (MQ=255)                                                                                                                                                                                                                                                       cccACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATTACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGcaca                                    <  2:495353/301‑1 (MQ=255)                                                                                                                                                                                                                                                                   aaTGCGTTTAGGCGGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGAAAACATCACTGGTCTATCGCAGGCAAAAGATAGTAATTGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTtac                        <  1:521378/300‑1 (MQ=255)                                                                                                                                                                                                                                                                    aTGCTTTCAGGTAGTATCTCGTATGCCTTGGTTTACCATACCGCAAAGGGGAGTTCAGGAAAGTATCAAGCGCAGTCAGACACCGTCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTGTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTtact                       <  1:465568/301‑1 (MQ=255)                                                                                                                                                                                                                                                                       cTTTCAGGCTGTCTCTCGTATGCCCTGTTCTCCTATAAAGAAAAGGTTAGTTCAGGAAAGTATCAATTGCAGTCAGACACCATCACTCGGCTATCGCAGGCATAAGATAGTAATGGAACAAAAAGCCATGTTTTTGCAGGGGAATCGCCGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCTAGATTTCTTTTAAATGCGGGGAAAAAGAAATTTAACGGCACAATATCGCTTtacta                      >  2:53370/1‑299 (MQ=255)                                                                                                                                                                                                                                                                       cTTTCAGGCTGTATCTCGTATGATCTGGTTTCCCATACCGCAAAGGGTAGTTAAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTtactac                     <  2:302736/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                          gggCTCTGGATTCCCATACAGCAAAGGATAGTTCAGGAAAGTCGCAATCGAAATCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAAAACAAGGCTAAGTTTTTGTAGCGGAATCACTGGATTAACTTATCACTGATGGAGCCAATGCCGTAGTTAAGATGGTCAATGATCCAGCACTTAACCGGCACAATATTCAAGTTGTCGATGACGCACGTTTTGTTTTAAATGCGGGGAGAAAGAAATTTGCCGCCACCATATCGCTTTACGACTACTGGAATACCgaagaag  <  1:54383/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                        AGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAG  >  NC_000913/3875236‑3875816

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