Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,830,893	A→G	N146S (AAT→AGT)	yicH →	AsmA family protein YicH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,830,893	0	A	G	98.0%	183.5 / NA	51	N146S (AAT→AGT)	yicH	AsmA family protein YicH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (26/24);  minor base T (0/1);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGC  >  NC_000913/3830660‑3831134                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  gTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCAAAAAGTGTCGCCATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGCCCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGCAGGGGGCGGCGTAGTTCCGt                                                                                                                                                                                                                                  >  2:243317/1‑251 (MQ=255)          tGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACCGAGCGCGCAGCGGGTAAGTGGCGGCGGAGTTCCGTGGTCACCa                                                                                                                                                                                                                          >  2:228951/1‑250 (MQ=255)                ccAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGCCCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCg                                                                                                                                                                                                                   <  2:366697/251‑1 (MQ=255)                    cccGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGACGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTa                                                                                                                                                                                                                >  2:31945/1‑250 (MQ=255)                              tggtggCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGCGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGt                                                                                                                                                                                                     >  1:196043/1‑251 (MQ=255)                              tggtggCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGt                                                                                                                                                                                                     >  1:56923/1‑251 (MQ=255)                               ggtggCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGt                                                                                                                                                                                                     <  1:352447/250‑1 (MQ=255)                                         aGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCa                                                                                                                                                                                          <  1:31945/251‑1 (MQ=255)                                          gtgtCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAg                                                                                                                                                                                         <  2:153179/251‑1 (MQ=255)                                                aCATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTc                                                                                                                                                                                     >  2:192763/1‑249 (MQ=255)                                                              aGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTTCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATc                                                                                                                                                                     <  1:16368/251‑1 (MQ=255)                                                                     aaCNGACCGAACCACGCCATGTCGCTAGCATGCTGCTGGAAAACGGGACGCGGAATCTCACCGACCAGACCGCCCCGCTCCCGTTCAAAGCCGGTCGTCTGCAACGGCGTGATGTGGCGTTGAATAGCTGGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGTCGTAGTTCCGTGGTCACCAAAAGGCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCg                                                                                                                                                              <  2:370527/251‑1 (MQ=255)                                                                         gACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCa                                                                                                                                                          <  1:324508/251‑1 (MQ=255)                                                                                   cGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTcc                                                                                                                                                >  1:42891/1‑251 (MQ=255)                                                                                     ccATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTg                                                                                                                                              <  1:239507/251‑1 (MQ=255)                                                                                               aCCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGt                                                                                                                                    >  1:222091/1‑251 (MQ=255)                                                                                                            aaaaCGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGc                                                                                                                       >  2:19235/1‑251 (MQ=255)                                                                                                              aaCGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAg                                                                                                                     <  2:354688/251‑1 (MQ=255)                                                                                                                 gggACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTAt                                                                                                                  <  1:312508/251‑1 (MQ=255)                                                                                                                    aCGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGa                                                                                                               <  1:286040/251‑1 (MQ=255)                                                                                                                         gAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACg                                                                                                          >  2:276649/1‑251 (MQ=255)                                                                                                                               cACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATcgcg                                                                                                    >  2:398143/1‑251 (MQ=255)                                                                                                                                   gACCAGACAGCGGAGCTACCGCACCAAGCCGTTCGTCTGCAACTGCGTCAAATGGAGTTTAATAGCCCGAATAGCGAATGGAAACTGACCCCGCAGAGGGTAAGTGGCGGCGTAGTTCCGTGATCACCACAAGCCGGTAAAGTGCCGGGTACGAAGCCGCCGATTCAGTTTAGGGCCGGATCGCTTGCGCTCAATGATTTTCCTGCTACCAATGTACTGATTGAAGGCAGCCTTGATAACGATCGCGTTAc                                                                                                >  2:392095/1‑251 (MQ=255)                                                                                                                                          ccGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACt                                                                                         <  2:165560/251‑1 (MQ=255)                                                                                                                                             cgccgcTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAAc                                                                                      <  1:276649/251‑1 (MQ=255)                                                                                                                                             cgccgcTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAAc                                                                                      >  2:18481/1‑251 (MQ=255)                                                                                                                                                            aaGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATc                                                                       <  1:286393/251‑1 (MQ=255)                                                                                                                                                            aaGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATc                                                                       >  1:200561/1‑251 (MQ=255)                                                                                                                                                               ccGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGcc                                                                    <  2:413975/251‑1 (MQ=255)                                                                                                                                                                cgctCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGAGTGAAGGCAGTATTGATGTCGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGccc                                                                   <  2:165482/248‑1 (MQ=255)                                                                                                                                                                             cTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGAGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCACTGTCCTGATTGAAGGCAGTATTGATAACGATCGCCTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAAc                                                      >  2:50804/1‑251 (MQ=255)                                                                                                                                                                                  tGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGa                                                  >  1:304256/1‑250 (MQ=255)                                                                                                                                                                                   gATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGaa                                                 >  1:323098/1‑250 (MQ=255)                                                                                                                                                                                    atatGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAAc                                               <  1:81260/251‑1 (MQ=255)                                                                                                                                                                                       tGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGAAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGGGCCGACATCGCCCGCGGGACATTAACCGGAAAcgcg                                            >  2:22952/1‑251 (MQ=255)                                                                                                                                                                                             ttAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACAtt                                                         <  1:311259/232‑1 (MQ=255)                                                                                                                                                                                             ttAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGGGCCGACATCGCCCGCGGGACAtt                                                         >  2:311259/1‑232 (MQ=255)                                                                                                                                                                                              tAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCATTATTGATAACGATCGCGTTACGCTGACTATCCTGGGTGCCGACATCGCCCGCGTGACATTAACCGGAAACGCGCAGCGt                                      <  2:196043/250‑1 (MQ=255)                                                                                                                                                                                               aaTAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTa                                     >  1:130768/1‑250 (MQ=255)                                                                                                                                                                                                 tAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGCCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGCTTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAAACGGAAACGCGCAGCGCAAcc                                  >  2:166104/1‑250 (MQ=255)                                                                                                                                                                                                       aaaTAGCGAATGGAACCTGAGCGCGCAGCGGGGAATGGGCGGCGTAGGTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGAGTCAGTTTAGTGCCGGATCGCTGTCGCTCAATGAGGTGCCGGCCACCAATGTACTGGTTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAGCCTGGGTGCCGACATCGCGCGCGGGACATTAAGCGGAAACGCTGAGCGTAACGCCGACg                            <  2:44088/250‑1 (MQ=255)                                                                                                                                                                                                        aaTAGCGAATGGAACCTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTACGTGGACACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGAGCAGATTCCGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGCAACGCGCAGCGTAACGCCGACgg                           >  2:185971/1‑251 (MQ=255)                                                                                                                                                                                                         aTAGCGAATGGAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCATCAAAAGCCGGTACAGGGCTGGGTACGAAGGAGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGGCATTAACCGGAAACGCGCAGCGTAACGCCGACGGc                          >  2:141303/1‑251 (MQ=255)                                                                                                                                                                                                                  ggAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCa                  <  2:303558/250‑1 (MQ=255)                                                                                                                                                                                                                   gAAACTGAGCGCGCAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGAGGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAg                <  2:222091/251‑1 (MQ=255)                                                                                                                                                                                                                               gcAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACATgc    >  2:292544/1‑251 (MQ=255)                                                                                                                                                                                                                                cAGCGGGTAAGTTGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTGCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGTGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTgcg   <  2:89584/251‑1 (MQ=255)                                                                                                                                                                                                                                cAGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTgcg   <  1:398143/251‑1 (MQ=255)                                                                                                                                                                                                                                 aGCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGGGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTgcgc  <  1:58825/251‑1 (MQ=255)                                                                                                                                                                                                                                  gCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTgcgc  >  1:346720/1‑250 (MQ=255)                                                                                                                                                                                                                                  gCGGGTAAGTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTgcgc  >  1:115108/1‑250 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  GTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGC  >  NC_000913/3830660‑3831134

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