Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,213,765	0	A	C	26.8%	132.3 / 13.7	71	V120G (GTG→GGG)	yjaB	peptidyl‑lysine N‑acetyltransferase YjaB
Reads supporting (aligned to +/- strand):  ref base A (42/10);  new base C (0/19);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 2.31e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGT  >  NC_000913/4213516‑4214015                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            gTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTACCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCAc                                                                                                                                                                                                                                                            >  1:397322/1‑251 (MQ=255)   cGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATAGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACct                                                                                                                                                                                                                                                          >  2:519727/1‑251 (MQ=255)    gCTGGGGCAACTGTAAGGAAATATCGCCAGTTCCACAGTCATGCGCAAGCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGACCCGTAACCTTAAAACCCAccgt                                                                                                                                                                                                                                                         <  2:389838/251‑3 (MQ=255)    gCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACctt                                                                                                                                                                                                                                                         <  1:337976/251‑1 (MQ=255)       gggtcAACGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATATTGCGGTCATGGCAATACTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCACCGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACctt‑ctt                                                                                                                                                                                                                                                     <  2:179622/247‑1 (MQ=255)        gggCAACGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTTGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTta                                                                                                                                                                                                                                                    <  1:468123/251‑1 (MQ=255)          gCAACGTTAATGAAATATCGCCAGTTCCACACCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTtat                                                                                                                                                                                                                                                   >  1:702554/1‑250 (MQ=255)             aCGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGAACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGa                                                                                                                                                                                                                                                >  2:707764/1‑250 (MQ=255)              cGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAAc                                                                                                                                                                                                                                              >  1:114006/1‑251 (MQ=255)                         tatCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTg                                                                                                                                                                                                                                   >  2:114010/1‑251 (MQ=255)                            cGCCATTTCCACAACCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACTCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCCCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCt                                                                                                                                                                                                                                 <  1:334068/250‑1 (MQ=255)                               cAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCAtt                                                                                                                                                                                                                             >  1:620878/1‑251 (MQ=255)                                   gccACTTCCATGCGCAAACAGCGGTTCTCAGTGATAGGGCGGGTATGGCAATGCTAAAGCAGAAAAAATCGGGTCACCATGGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAAGCGCCTCCCCCCCCCCATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGCCCACCTCAGAGCGTCCCGTACCCTAAAACCCCCCCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTg                                                                                                                                                                                                                          <  2:81715/249‑1 (MQ=255)                                        aTCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGCACCCAACCGCCTGCTCATTTTGCTCAtt                                                                                                                                                                                                                    >  2:542562/1‑251 (MQ=255)                                          ccATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCAGGGCAATGCTTAAGCAGAAATAATCGTGCCACCTTTGGTGGTTACTAAACCTGATGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCCCATCCGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCCCCTCAGAGCGTCCCGTAACCTTAAACCCCCCCTT‑CTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTaa                                                                                                                                                                                                                  <  1:581498/251‑1 (MQ=255)                                                  aagcAGCGTTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTTTCACCATTGGTGGGTACTAAACCGGAAGTTCACCCCACCGGGATGAGAAAAAACCGCCTACGCCCCCACATACGCCAGATTCAGCAGCGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCCCCTT‑CTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACgttggtt                                                                                                                                                                                                          <  1:232432/248‑1 (MQ=255)                                                   aTCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACgttggttg                                                                                                                                                                                                         >  1:269635/1‑251 (MQ=255)                                                      aGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTTCTCATTAACGTTGGATGTc                                                                                                                                                                                                       >  2:535568/1‑250 (MQ=255)                                                            aCTCAGAGATAGTGCGGTCATGGCAATGCTTAAGCAGCAATAATCGTGTCACCATTGGTGGGTACTAAACCTGCAGTTACGCCCACCGGGCTGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCAAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGATCCg                                                                                                                                                                                                >  1:391982/1‑251 (MQ=255)                                                               cAGTGATAGTGCGGTCAAGGCAATGCTTAAGCAGAAATAATCGTTCCCCCATTGGTGGGTAATAAACTTGAAGTTCAGCCCCCCGGGATGAGAAAAAATCGCCTTCCCCCCCCCATACCCCAGATTAAGCAACGGATCCGTTTTCCCAAAACCGTCCACCTCAGAGCGTCCCGTA‑CCTTTAAACCCCCCTTTCTAAGAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCgg                                                                                                                                                                                               <  1:325420/249‑1 (MQ=255)                                                                aGTGATAGTGCGGTCATGGCAATGCTTAAGAAGAAATAATCGTGCCCCCATGGGTGGGTACTAACCCTGAAGTTCAGCCCCCCGGGATGAGAAAAAATCGCCTACGCCCCCACATCCGCCAGATTCAGCAACGGATACGTTTCCCCCAAATCGTCCCCCTCAGAGCGTCCCGTAACCTTAAACCCCCCCTT‑CTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGc                                                                                                                                                                                            <  1:114010/251‑1 (MQ=255)                                                                   gATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCAACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCAt                                                                                                                                                                                         >  1:670195/1‑251 (MQ=255)                                                                    agaGTGCGGTCATGGCAATGCTTAAGCAGAAAAAACCGTGTCCCCATTGGTGGGTACTAACCCTGAAGTTCAGCCCCCGGGGATGAGAAAAAACCGCCTACGCCCCCACATCCGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTACCCTTAAAACCCACCTT‑CTTATAGAACCCAGCCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATc                                                                                                                                                                                        <  1:4935/249‑1 (MQ=255)                                                                           cgcgtttccgatctgaaaGCAGAAATAATCGTGTCACAATGGGTGGGTACTAACCCTGAAGTTCAGCCCACCGGGAGGAGAAAAAATCGCCTACGCCCCCACAGCCGCCAGATTCAGCAACGGATCCGTTTCCCCCAAATCGTCCACCCCAGAGCGTCCCGTACCCTTAAATCCCCCCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAgcg                                                                                                                                                                                 <  2:264061/235‑1 (MQ=255)                                                                                 ggCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCt                                                                                                                                                                           >  1:378016/1‑251 (MQ=255)                                                                                       gCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTccacc                                                                                                                                                                      >  1:261681/1‑250 (MQ=255)                                                                                       gCTTAAGCAAAAAAAATCGTGTCCCCATGGGGGGGTACTAACCCGGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTCCGCCCCCACATCCGCCAGATTCAGCAACGGATACGTTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCCCCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTccacca                                                                                                                                                                     <  1:206399/251‑1 (MQ=255)                                                                                        cTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGGAGTTCAGCCCACCGGGCTGAGAAAAAATCGCCTACGCCCACACATACGCCAGATTCAGCAACGGATACGGTTTCCCAAAATAGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACACAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTccaac                                                                                                                                                                      >  1:645722/1‑247 (MQ=255)                                                                                         ttAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAgcg                                                                                                                                                                                 >  1:264061/1‑237 (MQ=255)                                                                                           aaGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCAc                                                                                                                                                                 >  1:563970/1‑251 (MQ=255)                                                                                                  aataatCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCt                                                                                                                                                          >  1:173074/1‑251 (MQ=255)                                                                                                     aagcgggtCACCATTGGTGGGTACTAAACCATATATTCAGCCCGCCGGGAGGAGAAAAAACCGCCTCCCCCCCCACATACCCCAGATTCAGCAACGGATACGGTTCCCCCAAATCTCCCACCTCAGGGTGCCCCGTACCCTTAAAACCCCCCTT‑TTTATAGAACCCAGCCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTAc                                                                                                                                                        <  2:106467/244‑1 (MQ=255)                                                                                                         gtgtCACCATTGGGGGGTACTAAAACTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGAAACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGCGAGCGCATGCTCCACCAGCACCCGACCTAcgccgc                                                                                                                                                   >  2:530928/1‑251 (MQ=255)                                                                                                              aCCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCcg                                                                                                                                              >  1:252528/1‑251 (MQ=255)                                                                                                                 aTTGGTGGGTACTAAACCGGAAGTCCAGCCCCCCGGGATGAAAAAAAACCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTCCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAACCCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCgcgc                                                                                                                                           <  1:362366/251‑1 (MQ=255)                                                                                                                   tggggggTACGAAACCTAAAGGTCCGCCCACCGGGAAAAAAAAAAACCCCCGCCCCCCCCAAATACGCCAGATTCCCCACCGGACACTTTCTCCCCAAATCGCCCCCCCCGAAGCGCCCCGCACCCTTAAACCCCCCCTT‑CTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACATTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGca                                                                                                                                          <  2:107712/246‑1 (MQ=255)                                                                                                                    ggggggTACTAAACCTAAAGTCCACCCCACCGGGAGGAAAAAAAACCCCCTCCGCCCCCACATACGCCAATTTCAGCAACGGATACTGTTTCCCAAAATCGTCCACCTCAGCGCGTCCCTGAACCTTAACCCCCCCCTT‑CTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTCCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCtcgccc                                                                                                                                         <  2:549624/247‑3 (MQ=255)                                                                                                                      ggggtccaaaaCCTGAAGTTCACCCAACGGGGAGGAGAAAAAACCGCCTCCGCCCCCACATCCGCCAGATTCAGCAACGAATACGGTTCCCCCAAACCGTCCACCTCAGATCGCCCCGTACCCTAAAACCCCCCCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACTTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATc                                                                                                                                      <  2:620878/243‑1 (MQ=255)                                                                                                                        ggTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAACCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCCCCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAg                                                                                                                                    <  1:546025/251‑1 (MQ=255)                                                                                                                                      aGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACTGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAgc                                                                                                                      <  2:707213/251‑1 (MQ=255)                                                                                                                                       gTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAgcg                                                                                                                     >  1:174146/1‑251 (MQ=255)                                                                                                                                         tCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAgcgc                                                                                                                    >  1:435756/1‑250 (MQ=255)                                                                                                                                             cccACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCCGCCGCGCCCCTCAGGCTCGATACACACCGCATccc                                                                                                               >  2:11859/1‑250 (MQ=255)                                                                                                                                                        gagaAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATCCGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTACCCTTAAAACCCCCCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGccc                                                                                                     <  1:147281/250‑1 (MQ=255)                                                                                                                                                            aaaacgcGCCTACGCCCCCCGATACGGCAGTTTCATCAGTGGATACGGTTCCCCCAAATCGTCCCCCTCAGCGCGTCCCGTAACCTTAAACCCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCTCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTa                                                                                                <  1:451341/245‑1 (MQ=255)                                                                                                                                                                       ccGCCCCCACATACGCCAGATTCAGCAACGGATACGTTTTCCCCAAATCGTCCACCTCAGAGCCTCCCGTAACCTTAAAACCCCCCTT‑CTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATaaa                                                                                      <  1:585078/249‑1 (MQ=255)                                                                                                                                                                        cGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATc                                                                                    >  1:544720/1‑251 (MQ=255)                                                                                                                                                                           ccccccATACGCCAGATTCATCAACGGATCCGGTTTCCCCAAATCGTCCCCCTCAGGGCGCCCCGTACCCTTAAACCCCCCCTT‑CTTATGGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAgcgc                                                                                                                                                                                <  1:466308/156‑1 (MQ=255)                                                                                                                                                                           ccccACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAgcgc                                                                                                                                                                                >  2:466308/1‑156 (MQ=255)                                                                                                                                                                               acaTACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCgg                                                                             >  1:348264/1‑251 (MQ=255)                                                                                                                                                                                  tACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACCGCGCATCCATATGCNGCCCACTTAGCAACATAAATCCAACCGGCTg                                                                          >  2:324634/1‑251 (MQ=255)                                                                                                                                                                                       cAGATTCAGCAACGGATACGGTTTCACCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCAAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGATGaaacc                                                                     >  1:375727/1‑247 (MQ=255)                                                                                                                                                                                        aGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCg                                                                    >  1:156805/1‑251 (MQ=255)                                                                                                                                                                                        aGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCg                                                                    >  2:522016/1‑251 (MQ=255)                                                                                                                                                                                                        aCGGTTCCCCCAAAGCGCCCCCCCCAGAGCGTCCCGTAACCTAAAAACCCACCTT‑CTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGAc                                                     <  2:269635/250‑1 (MQ=255)                                                                                                                                                                                                         cGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATACACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGAccc                                                   >  1:130158/1‑251 (MQ=255)                                                                                                                                                                                                         cGGTTTCCCCAAACCGTCCCCCTCAGAGCGTCCCGTAACCTTAAAACCCCCCTG‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGAccc                                                   <  1:542562/251‑1 (MQ=255)                                                                                                                                                                                                          ggTTTCCCCAAATCGTCCCCCTCAGAGCGTCCCGTACCCTTAAAACCCCCCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGcc                                                                                                                                               <  1:465749/158‑1 (MQ=255)                                                                                                                                                                                                          ggTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGcc                                                                                                                                               >  2:465749/1‑158 (MQ=255)                                                                                                                                                                                                          ggTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCat                                                                                                              >  1:554445/1‑191 (MQ=255)                                                                                                                                                                                                          ggTTCCCCCAAATCGTCCCCCTCAGAGCGTCCCGTACCCTAAAACCCCCCCTT‑TTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCat                                                                                                              <  2:554445/191‑1 (MQ=255)                                                                                                                                                                                                                  cAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACgccg                                          >  2:71164/1‑248 (MQ=255)                                                                                                                                                                                                                             cccccacgcgccccttacccTTAAAACCCCCCTC‑CTTATAGAACCCACCCGCCTTCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCAGCGAGAGCGCATGCTCCCCCAGCACCCGACCTACGCCGCAGCCGCGCACCTCAGGATCGATAAACGGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACACCGGCGCTTCCTGCagg                               <  2:413116/233‑1 (MQ=255)                                                                                                                                                                                                                                           gTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACt                   >  1:141029/1‑249 (MQ=255)                                                                                                                                                                                                                                              aCCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCCGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTAATTAACAGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGtc              >  1:266143/1‑251 (MQ=255)                                                                                                                                                                                                                                               ccTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTCACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGTCGCTTCCGTCAGGAAAGACGGAACTAGGtcc             >  2:483411/1‑251 (MQ=255)                                                                                                                                                                                                                                                cTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGtcct            >  1:509931/1‑251 (MQ=255)                                                                                                                                                                                                                                                   aaaaCCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGtcctcc          >  1:260536/1‑250 (MQ=255)                                                                                                                                                                                                                                                   aaaaCCCACCTT‑CTAATAGCACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGCCCTACGCCGCAGCCGCGCACATCAGGAGCGATAAACAGCGCATCCTTATGCCGCCCACTTAGAAACATAAATCCAACCCGCTGATCCAGCTCATTAACCGCGACCCACAACTCCGCTTACCGCAGGAAGGAACACACTACGtcctcc          >  2:684712/1‑250 (MQ=255)                                                                                                                                                                                                                                                     aaCCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGc       >  1:457102/1‑251 (MQ=255)                                                                                                                                                                                                                                                      aCCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCt      >  1:616139/1‑251 (MQ=255)                                                                                                                                                                                                                                                          aCCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTCTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCCTCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGACCGAACTAGGTCCTCCAGCTCggc  >  1:425218/1‑250 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            GTCGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTT‑CTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGT  >  NC_000913/4213516‑4214015

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