Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,292,312	C→T	P324P (CCC→CCT)	yraK →	putative fimbrial adhesin YraK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,292,312	0	C	T	98.0%	167.8 / NA	50	P324P (CCC→CCT)	yraK	putative fimbrial adhesin YraK
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (29/20);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCCGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGC  >  NC_000913/3292069‑3292552                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  caccCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATc                                                                                                                                                                                                                                           <  1:138454/251‑1 (MQ=255)    ccGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCtt                                                                                                                                                                                                                                        <  2:22199/251‑1 (MQ=255)        aaTGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGtt                                                                                                                                                                                                                                                <  2:136395/240‑1 (MQ=255)        aaTGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGtt                                                                                                                                                                                                                                                >  1:136395/1‑240 (MQ=255)                gccaGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGGGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACATCCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACTGTAGCCTTCCTGTT‑GATCCTTGCGGGCATGGCa                                                                                                                                                                                                                            <  2:287527/251‑1 (MQ=255)                       ggTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGCGATAATTCAAAGATTAAAAACTCGTTGCCAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCAGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCa                                                                                                                                                                                                                      >  2:400530/1‑250 (MQ=255)                             actctaatAGACTTTGATATCACTAATGACTGTGTTGATACTTCAAAGATTAAAAACTCGTGGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGAACAGTCCAACCTGGTCGCCATGCGAAGAAGTGCAGACAATGTGCCCGTTGTCGGTATTCGTATTGAAAATTTCGGCGGCGGAGTGGCAAATATTCCTTTTTAGAATGGTATACTTCCTGCT‑GATCCTGCCGGGCATGCCAGCGTCAATATgcg                                                                                                                                                                                                               <  2:309060/244‑1 (MQ=255)                             aCTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATgc                                                                                                                                                                                                                <  1:165879/250‑1 (MQ=255)                                    tAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTgg                                                                                                                                                                                                        >  2:166162/1‑251 (MQ=255)                                      gACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGACAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGcc                                                                                                                                                                                                      >  2:119007/1‑251 (MQ=255)                                               atCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAACTGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGCAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCACATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGTGCATGGCACCGTCAATATGCGCGCCTGCCCAGTTAATCt                                                                                                                                                                                             >  1:212090/1‑251 (MQ=255)                                                 cACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTgg                                                                                                                                                                                           <  1:196262/251‑1 (MQ=255)                                                      aTGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCggt                                                                                                                                                                                      >  1:122327/1‑251 (MQ=255)                                                      aTGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGACTGGCCAGTTAATCTGGtaggg                                                                                                                                                                                      >  1:275418/1‑247 (MQ=255)                                                        gACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAAGCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGGCggtg                                                                                                                                                                                     <  2:181295/250‑1 (MQ=255)                                                                    aCTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGCTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTTCAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAAc                                                                                                                                                                        >  1:226555/1‑251 (MQ=255)                                                                    aCTTCAAAGATAAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGaaa                                                                                                                                                                         <  1:5324/250‑1 (MQ=255)                                                                        cAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCCGGCGACGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCATTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAgg                                                                                                                                                                     >  1:170113/1‑250 (MQ=255)                                                                                aaaaaCTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCATTATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTc                                                                                                                                                             >  1:375649/1‑250 (MQ=255)                                                                                             aaaTGCGCATCGACGGGACAACTGGGGTAGTAGACCAGTACAACCGGGTCGCCAGGCGAAGAAGTTCAGACAATGGGCCCGAGGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTGTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGGGGTGAGCTGGAAACAGGAAAATTGCAGGGCACAGccccc                                                                                                                                               <  2:212090/251‑4 (MQ=255)                                                                                                  cgcATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGGCGGTATTCGTATTGAAAATCTCGGCGGCGTAGTTGCAAATATTCCTTTTCAGAACGGTATTCTTGCTTTT‑GATCCTGCCGGGCATGGCACCGTCAAGATTCGCGCCTTGCCATGTAATCTTTTCGGTGGTGGGCTgg                                                                                                                                                                            <  2:255949/217‑1 (MQ=255)                                                                                                  cgcATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGAGAAGAAGTTCAGACAATCTGCCCGCCGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCACATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTACTCTGGTCGGTGGTGAGCTgg                                                                                                                                                                            >  1:255949/1‑217 (MQ=255)                                                                                                    cATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCg                                                                                                                                        <  1:104122/251‑1 (MQ=255)                                                                                                                    ggggTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCTGATGGCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTGGCAGAACGGTATCCTTCCTGGT‑GATCCGTCCGGGCATGGCACCGGCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAAGTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTaaat                                                                                                                        <  2:87332/251‑2 (MQ=255)                                                                                                                                  gTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGtt                                                                                                          >  1:138056/1‑251 (MQ=255)                                                                                                                                       aCCTGGTCGCCAGGAGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGAt                                                                                                     >  2:416535/1‑251 (MQ=255)                                                                                                                                               gCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTaacgca                                                                                              >  2:266694/1‑250 (MQ=255)                                                                                                                                                 cAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTaacgcaacg                                                                                           <  2:183689/251‑1 (MQ=255)                                                                                                                                                 cAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTaacgcaacg                                                                                           >  2:90838/1‑251 (MQ=255)                                                                                                                                                      gaagaaGCTCCGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATATCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCg                                                                                      >  2:293782/1‑251 (MQ=255)                                                                                                                                                      gaagaaATTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAACCTCGGCGGCGGAGTTGCAAATATTCCTTGTCAGAACGTTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCGCCGGCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGTAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGGTTCGTTGTGATTatcgcgtggcttcg                                                                                      <  2:46135/251‑10 (MQ=255)                                                                                                                                                         gaagTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAACATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGc                                                                                    >  2:264161/1‑250 (MQ=255)                                                                                                                                                                  aCAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTAccc                                                                          >  1:311494/1‑251 (MQ=255)                                                                                                                                                                            cGATGTCGGTATTCGTATAGAAAATCTCGGCGGCGGAGCTGCAAATCTTCCTTTTCAGAACGGTATTCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTACTCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCCTTACCGTCATCGTGCGGTAAAAGCAGACGATGCGTTTTGATTAACGAAAAGCTTCGGCACTGTTACCCCTGCTGCTcc                                                                >  2:38624/1‑251 (MQ=255)                                                                                                                                                                                 tCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAgcgc                                                           >  1:66565/1‑251 (MQ=255)                                                                                                                                                                                      aTTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACt                                                      <  1:266694/251‑1 (MQ=255)                                                                                                                                                                                          gTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTtata                                                  <  1:293782/251‑1 (MQ=255)                                                                                                                                                                                                aaaaTCTCGGCGGAGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAgcgc                                            >  2:410882/1‑251 (MQ=255)                                                                                                                                                                                                 aaaTCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCa                                           <  1:166162/251‑1 (MQ=255)                                                                                                                                                                                                       cggcggcggAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCCATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTttc                                     >  1:408375/1‑251 (MQ=255)                                                                                                                                                                                                       cggcggcggAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTttc                                     >  1:305547/1‑251 (MQ=255)                                                                                                                                                                                                            gcggAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCAc                                <  1:416535/251‑1 (MQ=255)                                                                                                                                                                                                                     cAAATATTCCTTTTCAGATGGGTATCCT‑CCCGTTTGACCCTTCCGGGCATTGCACCGTCAATATGAGAGCCTGGCCAGATAATCGGTTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACTTTGCGTTTTGATTAACGCATCGCTTCGTCACTGTTACCCCTGCTGCTCCATCGCATACTTATACAGCGCATTTTTCGTCACGCCGTGAAt                       <  2:404159/251‑1 (MQ=255)                                                                                                                                                                                                                      aaaTATTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAAtt                      >  2:312394/1‑251 (MQ=255)                                                                                                                                                                                                                        atatTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCa                                                                                                                                     >  2:383123/1‑138 (MQ=255)                                                                                                                                                                                                                        atatTCCTTTTCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCa                                                                                                                                     <  1:383123/138‑1 (MQ=255)                                                                                                                                                                                                                                 ttCAGAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGcc             >  1:96949/1‑249 (MQ=255)                                                                                                                                                                                                                                     gAACGGTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGcgcc       >  1:373386/1‑251 (MQ=255)                                                                                                                                                                                                                                          gTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGcgccgccgc  >  1:353095/1‑251 (MQ=255)                                                                                                                                                                                                                                          gTATCCTTCCTGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGcgccgccgc  >  2:122339/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  CACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCCGTT‑GATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGC  >  NC_000913/3292069‑3292552

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: ‑