Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,163,318	C→T	P20L (CCG→CTG)	thiK →	thiamine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,163,318	0	C	T	98.8%	326.7 / ‑4.6	86	P20L (CCG→CTG)	thiK	thiamine kinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (39/45);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AGCGCCTCTGGCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGC  >  NC_000913/1163175‑1163459                                                                                                                                                |                                                                                                                                               aGCGCCTCTGGCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTat                                                                                                                                         <  4:641941/151‑1 (MQ=255)     ccTCTGGCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTatcatc                                                                                                                                     >  2:904356/1‑151 (MQ=255)     ccTCTGGCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTatcatc                                                                                                                                     >  4:1683412/1‑151 (MQ=255)      ctctGGCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATcc                                                                                                                                    >  2:933110/1‑151 (MQ=255)       tctGGCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATccc                                                                                                                                   >  4:1448593/1‑151 (MQ=255)        ctGGCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATcc                                                                                                                                    <  4:422998/149‑1 (MQ=255)         tGGCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGt                                                                                                                                 >  4:693523/1‑151 (MQ=255)          ggCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCg                                                                                                                                  <  2:795717/149‑1 (MQ=255)          ggCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGt                                                                                                                                 <  4:111251/150‑1 (MQ=255)          ggCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTc                                                                                                                                <  4:35293/151‑1 (MQ=255)           gCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCa                                                                                                                               <  4:587957/151‑1 (MQ=255)            cAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCac                                                                                                                              <  2:1623248/151‑1 (MQ=255)             aaCGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCacg                                                                                                                             <  2:973326/151‑1 (MQ=255)                 ttAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGt                                                                                                                         >  2:529361/1‑151 (MQ=255)                   aaCGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGttt                                                                                                                       <  3:266468/151‑1 (MQ=255)                   aaCGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGttt                                                                                                                       <  2:1185993/151‑1 (MQ=255)                    aCGCTCCGACCCTACAAAGGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGtt                                                                                                                        <  1:419691/149‑1 (MQ=255)                     cGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTaa                                                                                                                     <  3:265099/151‑1 (MQ=255)                      gCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAAt                                                                                                                    <  4:931590/151‑1 (MQ=255)                       cTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAAt                                                                                                                    >  3:252955/1‑150 (MQ=255)                       cTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAAt                                                                                                                    <  4:1011605/150‑1 (MQ=255)                        tCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAAt                                                                                                                    >  2:1334956/1‑149 (MQ=255)                        tCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAAt                                                                                                                    >  2:45864/1‑149 (MQ=255)                         ccGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAg                                                                                                                  <  2:980861/150‑1 (MQ=255)                           gACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTgg                                                                                                               >  2:1072242/1‑151 (MQ=255)                           gACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTgg                                                                                                               >  1:177578/1‑151 (MQ=255)                           gACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTgg                                                                                                               <  1:231106/151‑1 (MQ=255)                           gACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTgg                                                                                                               >  3:321393/1‑151 (MQ=255)                              ccTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCt                                                                                                            >  2:1565442/1‑151 (MQ=255)                                tACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCtt                                                                                                           >  4:109220/1‑150 (MQ=255)                                 aCAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCtt                                                                                                           >  1:121633/1‑149 (MQ=255)                                 aCAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCtt                                                                                                           >  1:252608/1‑149 (MQ=255)                                  gaaaTGCAGATGATGCTGGTGGGGAAGAGCGAAATTATCTGGTTAGGTAAAGGTGGGGTTTCGCAGCAATAATAAAATAAAGAGAGAGGAATTGCTGTCGCGATTTTTACGGCAGTATCATATCGTCACGACGTTTAATAGTGGGCATAg                                                                                                         >  2:488877/2‑150 (MQ=255)                                  cAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGt                                                                                                        >  4:495068/1‑151 (MQ=255)                                    aaTGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTgg                                                                                                      <  4:1131132/151‑1 (MQ=255)                                    aaTGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTgg                                                                                                      <  2:236698/151‑1 (MQ=255)                                    aaTGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTgg                                                                                                      <  1:296566/151‑1 (MQ=255)                                     aTGCAGCTGATGCTGGTGCGGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGc                                                                                                     >  2:1256781/1‑151 (MQ=255)                                     aTGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGc                                                                                                     <  4:1042356/151‑1 (MQ=255)                                     aTGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGc                                                                                                     <  2:1287827/151‑1 (MQ=255)                                     aTGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGc                                                                                                     <  1:427724/151‑1 (MQ=255)                                       gCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑gg                                                                                                  >  2:1511124/1‑151 (MQ=255)                                        cAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑ggg                                                                                                 >  4:1056091/1‑151 (MQ=255)                                           cTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑ggg                                                                                                 >  3:271588/1‑148 (MQ=255)                                            tGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑ggggg                                                                                               >  2:1599516/1‑147 (MQ=255)                                             gATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGtt                                                                                             <  2:765038/150‑1 (MQ=255)                                             gATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGtt                                                                                             <  2:721835/150‑1 (MQ=255)                                             gATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGt                                                                                              <  4:1050353/149‑1 (MQ=255)                                                  tGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGgttttt                                                                                         >  4:1653190/1‑147 (MQ=255)                                                   ggTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCAt                                                                                      <  3:475625/151‑1 (MQ=255)                                                     tGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCCttt                                                                                    >  2:1372558/1‑150 (MQ=255)                                                     tGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGGTTTTTCCtt                                                                                     >  3:31959/1‑150 (MQ=255)                                                        aGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGCGgttttttcttt                                                                                    >  3:247413/1‑141 (MQ=255)                                                        aGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAAc                                                                                 <  4:694317/151‑1 (MQ=255)                                                          aCGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACAt                                                                               <  2:678917/151‑1 (MQ=255)                                                               cGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCTTTTTACATc                                                                              >  2:1292414/1‑147 (MQ=255)                                                                 aaaTTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGcc                                                                        <  2:447173/151‑1 (MQ=255)                                                                     tATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGc                                                                     <  2:52972/150‑1 (MQ=255)                                                                      aTCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCCTTTTACATCACGGCCcgcg                                                                    >  4:1313843/1‑150 (MQ=255)                                                                          ggTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGt                                                                   <  4:175919/147‑1 (MQ=255)                                                                              aGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTtgttgt                                                           <  2:1529906/151‑1 (MQ=255)                                                                               ggTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTtgtg                                                          <  4:768503/151‑1 (MQ=255)                                                                                gTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGTGCCAGCGTTTTGTTGTGc                                                         <  4:557168/151‑1 (MQ=255)                                                                                gTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGGTTTTTTTTTTTACATCAAGGcccacggtttttt                                                             >  4:730743/1‑136 (MQ=255)                                                                                    tggTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTGTTTTCTTTTTAACTCCAGGCCCGCGTtttttttttggtcg                                                     >  2:101738/2‑143 (MQ=255)                                                                                       tGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTTATTTAACATCTGGGCCCGCGGTTTTTTTTgtgtgcgcc                                                  >  2:1219870/1‑146 (MQ=255)                                                                                        gCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGcc                                                  <  2:231140/150‑1 (MQ=255)                                                                                          cGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTTTTTTAAAATCAAGGCCAGCGGTTTTTtgtgggggagac                                                  >  3:25443/1‑140 (MQ=255)                                                                                          cGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCa                                               <  1:89220/151‑1 (MQ=255)                                                                                           gTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCAc                                              <  2:1515210/151‑1 (MQ=255)                                                                                            tttCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCTTTTTACATAAAGGCCAGCGGTTTTTTTTGTGTCAGcc                                                  >  4:3237/1‑146 (MQ=255)                                                                                             ttCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACg                                             <  2:1266278/150‑1 (MQ=255)                                                                                                  agcaATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCCTTTTACATCCAGGCCAGCGGTTTTTTGTTGGTCAGcc                                                  >  1:509/1‑140 (MQ=255)                                                                                                   gcaATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGa                                      <  4:545493/151‑1 (MQ=255)                                                                                                      ataatCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGTTTTTTCTTTTTAAATTAAGGCCAGCGTTTTTTTGTGGGGTAGCCGcccccacccgatg                                    >  2:486973/1‑138 (MQ=255)                                                                                                         atCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTTTTTTAACATCAAGGCCAGCGGTTTTTTTTGGGTGGGGCGCCCGATCCttgtgtg                                  >  1:447818/1‑143 (MQ=255)                                                                                                         atCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCATTTAACAAAAGGGGCAGCGGCTTGTTGTGGGTCAGCCgggcggacctcat                                     >  2:233352/1‑134 (MQ=255)                                                                                                                cACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTcc                         <  4:1390968/151‑1 (MQ=255)                                                                                                                    cgcgACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑TTTTGTTTTCTCATTGACCATCAGGGCCAGCTTTTTTTTGTGCCTCACCCTCTCGACCCTGATGCGCCTCACTCcgcc                      <  2:481404/150‑2 (MQ=255)                                                                                                                      cgacgaATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGGGGTTTTTTCTTTTAA‑AACAGCGCCCGCGTTTTTTTGTTGGTTAGCCCCCCCAACCTGATGGGCCCCCGGCcgc                       >  4:1004047/1‑147 (MQ=255)                                                                                                                       gacgaATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCt                  <  2:755433/151‑1 (MQ=255)                                                                                                                          gaATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTTATTTTACAACAAGGGCAGCGGTTTTTTTTGGGTCAGCCGCCCCATCCCGCTGCGCCCCAGGCCGCGGTCtt                 >  4:551778/1‑149 (MQ=255)                                                                                                                          gaATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTgc               <  2:274886/151‑1 (MQ=255)                                                                                                                            aTTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTgcgc             <  4:14993/151‑1 (MQ=255)                                                                                                                                    cgcgCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCg     <  1:289503/151‑1 (MQ=255)                                                                                                                                       gcTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGc  <  4:372607/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                               AGCGCCTCTGGCAACGTTAACGCTCCGACCCTACAAATGCAGCTGATGCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGC  >  NC_000913/1163175‑1163459

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