Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,897,126	T→C	T249A (ACG→GCG)	ygcS ←	putative transporter YgcS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,897,126	0	T	C	97.3%	191.0 / NA	73	T249A (ACG→GCG)	ygcS	putative transporter YgcS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (38/33);  minor base A (2/0);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

AAACAGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCGGGTTTTCGGAAAT  >  NC_000913/2896882‑2897372                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       aaaaagcccatgtctacgcgaattGCTGCCAGCGGCGAAACACTTCCCAGCAAAGATTTGAGATGCGCCACCCGCTGCGTACGAACTAATCCCAGCAGCGCGCCCCCAATTACCAACGCATTAAGCATCTGACTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGAGAGTTGGCGGCCAGCCATAAATCACAAACCATGGGAGTACTAGGCAGAAAAAGACGACGCTGTTCAACGACGCGcgccg                                                                                                                                                                                                                                                  >  2:285145/23‑251 (MQ=255) aaaCAGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGcgccg                                                                                                                                                                                                                                                  <  1:397250/251‑1 (MQ=255)    cAGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCgccgcc                                                                                                                                                                                                                                                >  1:342207/1‑250 (MQ=255)      ggCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCAAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCCCGCGCCGCCAgg                                                                                                                                                                                                                                             >  2:184574/1‑250 (MQ=255)        ccATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAAGTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTaa                                                                                                                                                                                                                                           <  1:120666/251‑1 (MQ=255)         cATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAAGCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAAc                                                                                                                                                                                                                                          <  1:42162/251‑1 (MQ=255)            gACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAACTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAAcgc                                                                                                                                                                                                                                        <  1:518595/250‑1 (MQ=255)                     cGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCTCAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACa                                                                                                                                                                                                                              <  1:420061/251‑1 (MQ=255)                       ttGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACa                                                                                                                                                                                                                              >  1:328619/1‑249 (MQ=255)                              cagcagCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGATGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTg                                                                                                                                                                                                                      >  1:562636/1‑250 (MQ=255)                                 cagcaAAAAACTTCCCAGCAACCATTTGCGATGTGCAAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATAAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGAtgt                                                                                                                                                                                                                  >  1:439851/1‑251 (MQ=255)                                     aaaaaaCTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGttt                                                                                                                                                                                                              <  1:531279/251‑1 (MQ=255)                                                   cAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTa                                                                                                                                                                                                <  2:181016/251‑1 (MQ=255)                                                       aaTTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCAc                                                                                                                                                                                            >  1:475176/1‑251 (MQ=255)                                                             cGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATc                                                                                                                                                                                      <  1:189049/251‑1 (MQ=255)                                                                       gcagGTGCGTCGGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCCTCCAGACCAAGAGTCTGGGCGAGAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACTAGGCAGACAAAGAAGACGCTGGTAAACGCAGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTaaa                                                                                                                                                                            <  2:281304/251‑1 (MQ=255)                                                                                 cAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGAc                                                                                                                                                                  >  1:210204/1‑251 (MQ=255)                                                                                      cTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCaaa                                                                                                                                                              >  1:232427/1‑250 (MQ=255)                                                                                       tAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTAATCGCCCAGTAAAACATGGGGACCAAAAt                                                                                                                                                            >  1:473031/1‑251 (MQ=255)                                                                                         aTCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGAAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGGTGGCAGCCGGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCTGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTGTATGGGTCGCCGTTACCACTTCACAGCCCAGTAAAACATGGGGACCAAAAt                                                                                                                                                            >  2:162071/1‑249 (MQ=255)                                                                                           cccAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCg                                                                                                                                                        >  1:615520/1‑251 (MQ=255)                                                                                            ccAGCAGCGCGCCCACAAGAAACAACGCATGAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCGGTGCGATAGTTGGCAGCCAGGTATAAATCTCAAACCATGGGATGACGAGGCAGACAAAGAAGACGCTTTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTGGATTTGTTTATGGGACGCCGTGACCACTTCAGCGCCCAGGAAAGCATGGTGACCAGAATAGCGa                                                                                                                                                       <  2:61834/251‑1 (MQ=255)                                                                                                gaccgcgcCCCAAATTAACCACGCATGAAGGCTCAAGCTGGCAGTCAGCGCACTCTTCC‑GAATAAGAGTCTGTGCGATAGTTGGCAGCCAGGTATAAAACACAATCCAGGGGAATACGAGGCATACAAAGACGACACTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTGTTGATGTGTTGATGGGTCGCCGTTCCCACATCATCGCCCAGTAAAACATGGTGACCAAAAGAGCGATGCa                                                                                                                                                   <  2:249693/248‑1 (MQ=255)                                                                                                 atcacgcCCACAATTAACAACGCATTAAGCATCAGGATGGAAGTCAGAGCA‑TCTTCAAGAACAATAGTCTGGGCGATAGTTGGAAGCCAGGGATAAATCACAAACAATGGGATTACGAGGAAGACAAAGAAGACGCTGTTAAACGCCGCGAGCAGCAAGTAACGCGAAGAGAACAAGGTCTTGATGAGTTTATGGGTCGCCGTTACCACTACATCGACCAGAAAAACATGGGGACCAAAAAAGCGATGCa                                                                                                                                                   >  1:498387/5‑250 (MQ=255)                                                                                                        ccACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCg                                                                                                                                           >  2:190037/1‑251 (MQ=255)                                                                                                                 acaacGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCt                                                                                                                                  >  1:444413/1‑251 (MQ=255)                                                                                                                          tAAGCATCAGGCTGGCAGCCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACAAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGc                                                                                                                         >  1:112978/1‑251 (MQ=255)                                                                                                                            aGCATCAGGCGGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCCCTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGGAAAACGccc                                                                                                                       <  2:519616/251‑1 (MQ=255)                                                                                                                             gCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGcccc                                                                                                                      >  2:4955/1‑251 (MQ=255)                                                                                                                               agcAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTg                                                                                                                    <  1:319399/249‑1 (MQ=255)                                                                                                                                  aGGCTGGCAGTCAGCGCC‑TCTTACAGACCAATCGTCTGGGCGCTAGTTGGCAGCCAGGTATAAATCACAAACAATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGAGACCAAACTAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGgcg                                                                                                                 >  1:430820/1‑251 (MQ=255)                                                                                                                                  aGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGgcg                                                                                                                 >  2:234579/1‑251 (MQ=255)                                                                                                                                  aGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGgcg                                                                                                                 >  1:457774/1‑251 (MQ=255)                                                                                                                                  aGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGCCAGCCAGGTATAAATCACACCCCCTGGGATTGCGAGGCAGACACAGAAAACGCGGTTAAACGCCGCGCGCCGCCAGTACCGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGgcg                                                                                                                 >  1:250185/1‑251 (MQ=255)                                                                                                                                      tGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGATTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAAAGCCCCTGGCGCAGt                                                                                                             >  1:429562/1‑251 (MQ=255)                                                                                                                                             cAGCTCA‑TCTCCCAGACCAATAGTCTGGTGGATAGTGGGCAGCCAGGTAAAAATCACAAACCAGGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGAGGCACGATAGCGTGAGCTTCGGCAAAACGCCCCTGGCGCAGTAGCCAGc                                                                                                      <  2:125022/251‑1 (MQ=255)                                                                                                                                                cgcA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCgtg                                                                                                   <  2:232427/251‑1 (MQ=255)                                                                                                                                                 gcA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGAAGCCAGCgtgg                                                                                                  >  1:373874/1‑251 (MQ=255)                                                                                                                                                      cTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGAt                                                                                              >  1:200937/1‑251 (MQ=255)                                                                                                                                                         ccAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTc                                                                                            <  2:44523/250‑1 (MQ=255)                                                                                                                                                          cAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTg                                                                                          >  1:212753/1‑251 (MQ=255)                                                                                                                                                                        gggCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTGTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATATCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCg                                                                             <  1:407516/250‑1 (MQ=255)                                                                                                                                                                            gATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCaaaa                                                                                                                            <  1:168332/199‑1 (MQ=255)                                                                                                                                                                            gATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCaaaa                                                                                                                            >  2:168332/1‑199 (MQ=255)                                                                                                                                                                             aTAGTCGGCAGCCAGAAATAACCCACAAACCATGGGATTACGAGGCAGATACCGCAGACGCTGTTCAACGACGAGCGACGCCAGTAAAGCGCAGCGAACAAGGTTTTGATGCGTTTACGGGTCGCCGTTACCACTTCATCGCACAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTATGCAAAACGCCCCTGGCGCAGTAGCCAGCGCGGGGCTTATGGAGTTCCCCAGCgcaaga                                                                        >  2:407912/1‑245 (MQ=255)                                                                                                                                                                                gTTGGCAGCCAGGAATACATCACAAACCATGGGATTACGAGGCAGACCAAGCAGACGCTGTTAAACGCCGCGCGCCGCCAGCAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCACCTGGCGCAGTAGCCAGCGTGGTGCTTCTGGCGTTCCCCAGCGTAATAACGt                                                                    >  1:423995/1‑251 (MQ=255)                                                                                                                                                                                 ttGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGt                                                                    >  2:171081/1‑250 (MQ=255)                                                                                                                                                                                                 aaTCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGc                                                    >  1:144675/1‑250 (MQ=255)                                                                                                                                                                                                 aaTCACAAACCATGGGATTACGAAGCAAACAAAGCAGACGCTGTTAAACGCCGAGCGCCGCCAGTAACGCGAAGAGAACCAAGTTTTGATGCGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGCCCCAAATAGCGATGAACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGGTTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGc                                                    >  2:577302/1‑250 (MQ=255)                                                                                                                                                                                                      caAACCATGGGATTACGAGGCAGACAAAGAAGACGCTTTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCt                                              <  1:376751/251‑1 (MQ=255)                                                                                                                                                                                                      caAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCt                                              >  1:406923/1‑251 (MQ=255)                                                                                                                                                                                                        aaCCATGGGATTACGAGGCAGACAAAGAAGCCGTTGTTAAACTCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGc                                            <  2:25297/251‑1 (MQ=255)                                                                                                                                                                                                        aaCCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGACGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGc                                            >  1:285832/1‑251 (MQ=255)                                                                                                                                                                                                         aCCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGcc                                           <  1:336398/251‑1 (MQ=255)                                                                                                                                                                                                            aTGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGCAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGAt                                        >  1:252120/1‑251 (MQ=255)                                                                                                                                                                                                             tGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAAAATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATg                                       >  2:178965/1‑251 (MQ=255)                                                                                                                                                                                                                 aTTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTG‑TGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGt                                  >  2:358382/1‑251 (MQ=255)                                                                                                                                                                                                                      gAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTa                                                                                                                                                                                                             >  1:596419/1‑76 (MQ=255)                                                                                                                                                                                                                      gAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTa                                                                                                                                                                                                             <  2:596419/76‑1 (MQ=255)                                                                                                                                                                                                                           aGACAAAGAAGGCGCTTTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAATGTTTTGATGTGTTTATGGGTCGCCGTTACCTCTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAgcgc                         <  1:217791/251‑1 (MQ=255)                                                                                                                                                                                                                            gACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGcc                        >  2:149526/1‑251 (MQ=255)                                                                                                                                                                                                                               aaaaaagaCGCTGGTAACCGCCGCGCGCCGCCAGTAACGAGAAGAGAACAAGGTTTTGATGTGTTTATGGGTAGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAAGCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAgt                     <  2:23767/247‑2 (MQ=255)                                                                                                                                                                                                                                 agaagaCGCTGGTAAACGCCGCCCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGcc                   <  2:328619/251‑1 (MQ=255)                                                                                                                                                                                                                                  gaagaCGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGTTGCCAGTAGCCAGCGCCAGGCCt                  <  1:184574/251‑1 (MQ=255)                                                                                                                                                                                                                                   aagaCGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTc                 <  2:200937/251‑1 (MQ=255)                                                                                                                                                                                                                                    agaCGCGGTTAATCGCGGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTGGATGGGTGGTCGTTCCCACTTCATCGCCGAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGGCCCTGGCGCAGTAGCCAGCTTGGTGATTCCGGCGCTCCCCAGCGTAATATCGTAATCAACAACTCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCTGGCCTcc                <  2:250185/251‑1 (MQ=255)                                                                                                                                                                                                                                      aCGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCgg              <  1:207422/251‑1 (MQ=255)                                                                                                                                                                                                                                           gTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCGGGtttt         >  2:121225/1‑251 (MQ=255)                                                                                                                                                                                                                                              aaaCGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCGGGTTTTCgg      >  1:335734/1‑251 (MQ=255)                                                                                                                                                                                                                                              aaaCGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCGGGTTTTCg       <  1:337988/250‑1 (MQ=255)                                                                                                                                                                                                                                                aCGCCGCGCGCCGCCAGTATCGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCTCTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAAGAACGTAACCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCGGGTTTTCGGat    <  2:511898/251‑2 (MQ=255)                                                                                                                                                                                                                                                aCGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCGGGTTTTCGGaa    <  1:231159/251‑1 (MQ=255)                                                                                                                                                                                                                                                  gCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCGGGTTTTCGGAAAt  <  1:4955/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       AAACAGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCA‑TCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCGGGTTTTCGGAAAT  >  NC_000913/2896882‑2897372

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