Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	73,319	T→C	P401P (CCA→CCG)	thiP ←	thiamine ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,319	0	T	C	98.0%	154.7 / ‑4.1	49	P401P (CCA→CCG)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (35/13);  total (36/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

AACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCC  >  NC_000913/73074‑73565                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         aacaacGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGc                                                                                                                                                                                                                                                       >  1:353100/1‑249 (MQ=255)  acaacGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAgg                                                                                                                                                                                                                                                    >  1:19180/1‑251 (MQ=255)       gCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGccaat                                                                                                                                                                                                                                               >  2:221812/1‑248 (MQ=255)                  aaaaTCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGttgttg                                                                                                                                                                                                                                    >  1:345307/1‑250 (MQ=255)                   aaaTCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCTCTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCATGCCAATGG‑TGTTGTTGAg                                                                                                                                                                                                                                  <  2:326340/251‑1 (MQ=255)                          cAATCGACAGAACGCAGGAAAAAGCCAGCGCCTGCGCCAGTGGACGTGTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTTACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGT‑TGGTGTTGAGCAGTaaa                                                                                                                                                                                                                           <  2:166215/251‑1 (MQ=255)                           aaTCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTaaaa                                                                                                                                                                                                                          <  2:117182/251‑1 (MQ=255)                              cGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGGTGTTGAGCAGTAAAaag                                                                                                                                                                                                                        <  2:166012/250‑1 (MQ=255)                               gACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAaaga                                                                                                                                                                                                                       >  2:374877/1‑250 (MQ=255)                                          ggCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCa                                                                                                                                                                                                           <  2:321321/250‑1 (MQ=255)                                                agCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACaa                                                                                                                                                                                                     <  1:174333/250‑1 (MQ=255)                                                                   ggACGTTTCAGGGCGCGCAGCGCGACCACTGTTAAGCGCGACCAGCCTCCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCTCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTGGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCAGCATGGCGAGGAt                                                                                                                                                                                  <  2:270638/251‑1 (MQ=255)                                                                      cGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑Ca                                                                                                                                                                               >  2:269170/1‑250 (MQ=255)                                                                                  cgcAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCAt                                                                                                                                                                  >  2:33571/1‑251 (MQ=255)                                                                                         cgaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCa                                                                                                                                                             <  2:308006/247‑1 (MQ=255)                                                                                            ccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGcgcc                                                                                                                                                        >  2:309716/1‑248 (MQ=255)                                                                                               cTTTTAAGCGCGACCCGCCTTCAATGCCCAGCGACTGACACAAAATACTGTAGAGGGCGGTGATATCGCGCATCGGGTTTTCAAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATACCGTCAACAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACCATCACCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTg                                                                                                                                                      >  1:182538/1‑250 (MQ=255)                                                                                                                     aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCcgcg                                                                                                                               >  1:336400/5‑251 (MQ=255)                                                                                                                     aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGAGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCGCATCACCAGCGCCTGACCCGCCAGCCTTTTCTGCcgcg                                                                                                                               >  1:356498/5‑251 (MQ=255)                                                                                                                           cAGCGACTGAACCAACATACTGGAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGCCAATCACAATGCCGTCCGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGGAGCCAGCACACTCCCCGGCAAGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCCGACCCGCCAGCCTTTTCTGCCGCGCCCgac                                                                                                                         >  2:240693/1‑249 (MQ=255)                                                                                                                               gACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTc                                                                                                                     >  1:262059/1‑251 (MQ=255)                                                                                                                                                      ggcggTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGcatc                                                                                               >  2:124723/1‑250 (MQ=255)                                                                                                                                                        cggTGATATCGCGCGTCGGGTTCTCCAGTACTTTCAGAGCATAAGGGATCGCCATTAAAGCATTGGTGAAAATCAAACTGCCGTCAGCATATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGACAGCACAATCCCCGGCATGGCGAGGATACACCA‑GCCGCTCATCTCCAGAGCCTGACCCGCCAGCATTTTCTTCCGCGCCCGAAGTTCAAGACTGCGCCATCGCCGCATCATg                                                                                            >  2:243883/1‑251 (MQ=255)                                                                                                                                                           tGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATACCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATACCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTc                                                                                         >  2:185244/1‑251 (MQ=255)                                                                                                                                                             atatCGCCCATCGGGTTTGCGAGTACTTGCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGACAGCAGATTGCGGCAGCCCAAGGG‑TGTTGTTGACAAGAAAAAAAAAGCCGTTAGCCAGCACAACCCCCGGCTTGGCGAGGAT‑CAACATGCCGCGCATCTCCAGCCCCTGACCCGCCATCATTTGCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAg                                                                                       <  1:32135/251‑1 (MQ=255)                                                                                                                                                                    gcATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTAcc                                                                                >  2:230326/1‑250 (MQ=255)                                                                                                                                                                       tCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGcc                                                                                                                                  >  1:109276/1‑198 (MQ=255)                                                                                                                                                                        cGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTACCGCATTGGTGAACATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGCT‑AAACATGCCGCTCCTCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCCGCAGAAGTTCGCGACTGCTCCCTAGCAGCATCATGGACAGCACAACGccc                                                                             >  2:176587/1‑248 (MQ=255)                                                                                                                                                                         gggTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGcaca                                                                            >  1:191797/1‑250 (MQ=255)                                                                                                                                                                             tttCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGcc                                                                                                                                               <  1:23811/179‑1 (MQ=255)                                                                                                                                                                             tttCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGcc                                                                                                                                               >  2:23811/1‑179 (MQ=255)                                                                                                                                                                                     tcTTTAAGCGCAAAGGGGACCGCCATTAACGCATTGGTAAAAAACACAATGCCGTCAGCAGATTGCGGCAGGC‑AATGGGTTTTGTTGAGCAGAAAAAAGAAGCCGGAAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTgccg                                                                <  1:195088/249‑1 (MQ=255)                                                                                                                                                                                              gcATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAAt                                                      >  2:347080/1‑251 (MQ=255)                                                                                                                                                                                               cATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCCGCGCAAt                                                      >  2:65863/1‑250 (MQ=255)                                                                                                                                                                                                  aaGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGACAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGCCTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCCGCGCAATACGc                                                  >  2:22228/1‑251 (MQ=255)                                                                                                                                                                                                  aaGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATgg                                                                                                                                                                                         <  2:158285/117‑1 (MQ=255)                                                                                                                                                                                                  aaGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATgg                                                                                                                                                                                         >  1:158285/1‑117 (MQ=255)                                                                                                                                                                                                              ttAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCcac                                      >  2:358699/1‑251 (MQ=255)                                                                                                                                                                                                                aaCGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCACCGAGGTCcac                                      >  2:214889/1‑249 (MQ=255)                                                                                                                                                                                                                   gCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAgcgc                                 >  2:168279/1‑251 (MQ=255)                                                                                                                                                                                                                      ttGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCt                                                                                                                                                       >  1:5255/1‑130 (MQ=255)                                                                                                                                                                                                                      ttGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCt                                                                                                                                                       <  2:5255/130‑1 (MQ=255)                                                                                                                                                                                                                       tGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGc                             >  1:40136/1‑251 (MQ=255)                                                                                                                                                                                                                            aaaaTCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑GGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCCTAGCAGCATCATGGTCCGCACTACACACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCcacac                        >  2:139360/1‑250 (MQ=255)                                                                                                                                                                                                                               aTCACAATGCCGTCAGCAGATTGCGGCAGGCCCATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCCATACGCAACGAGGTCCACAGCGCCTGCCACAGCAc                     >  2:116534/1‑251 (MQ=255)                                                                                                                                                                                                                                      tGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTtgt              >  1:194401/1‑251 (MQ=255)                                                                                                                                                                                                                                             agcagATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGc        >  2:352094/1‑250 (MQ=255)                                                                                                                                                                                                                                                agATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGCCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGCGCCAGCACtt    >  2:266532/1‑251 (MQ=255)                                                                                                                                                                                                                                                  aTTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTcc  >  1:242189/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         AACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCC  >  NC_000913/73074‑73565

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