Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,747,168:1	+C	6.5%	coding (85/1497 nt)	lyxK →	L‑xylulose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,747,168	1	.	C	6.5%	369.6 / 21.5	92	coding (85/1497 nt)	lyxK	L‑xylulose kinase
Reads supporting (aligned to +/- strand):  ref base . (47/39);  new base C (4/2);  total (51/41)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

TGTTTGTTAAAGATGTGCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATC  >  NC_000913/3747019‑3747311                                                                                                                                                       |                                                                                                                                                tgtttgttAAAGATGTGCCGCAGGGAGCTGATCTGCTGAAAGCAGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCg                                                                                                                                                  >  2:965301/1‑150 (MQ=255)     tgttAAAGATGTGCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCa                                                                                                                                             <  1:347899/150‑1 (MQ=255)      gttAAAGATGTGCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAg                                                                                                                                            <  2:679273/150‑1 (MQ=255)             aTGTGCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGCTGCA‑CGCCTGc                                                                                                                                     <  2:951742/150‑1 (MQ=255)              tgtgCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCGCTGCA‑CGCCTGcc                                                                                                                                    <  1:1220960/150‑1 (MQ=255)              tgtgCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGcc                                                                                                                                    <  1:967916/150‑1 (MQ=255)               gtgCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGcca                                                                                                                                   <  2:753467/150‑2 (MQ=255)                 gCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCt                                                                                                                                 >  1:596117/1‑150 (MQ=255)                 gCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCt                                                                                                                                 >  1:1032553/1‑150 (MQ=255)                 gCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCt                                                                                                                                 >  2:700911/1‑150 (MQ=255)                             tGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATCCTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAg                                                                                                                     <  1:1576635/150‑1 (MQ=255)                                 ctgctgAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCg                                                                                                                 <  1:1327633/150‑1 (MQ=255)                                   gctgAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCa                                                                                                               >  2:132584/1‑150 (MQ=255)                                      gAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGcc                                                                                                            >  2:808572/1‑150 (MQ=255)                                      gAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGcc                                                                                                            >  2:961492/1‑150 (MQ=255)                                        aaGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAg                                                                                                          <  2:884524/150‑1 (MQ=255)                                        aaGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAg                                                                                                          >  1:1299818/1‑150 (MQ=255)                                          gCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGc                                                                                                                <  1:222744/142‑1 (MQ=255)                                          gCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGc                                                                                                                >  2:222744/1‑142 (MQ=255)                                             gTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTgg                                                                                                     <  2:468701/150‑1 (MQ=255)                                               ggATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGc                                                                                                   >  2:1010155/1‑150 (MQ=255)                                                gATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGct                                                                                                  <  1:100523/150‑2 (MQ=255)                                                  tGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCaga                                                                                                >  1:1444124/1‑150 (MQ=255)                                                       tGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGcgcg                                                                                           >  1:1467631/1‑150 (MQ=255)                                                       tGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGcgcg                                                                                           >  2:995890/1‑150 (MQ=255)                                                       tGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGcgcg                                                                                           <  2:152808/150‑1 (MQ=255)                                                       tGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGcgcg                                                                                           <  1:1101011/150‑1 (MQ=255)                                                       tGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGcgcg                                                                                           <  2:969893/150‑1 (MQ=255)                                                        gCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGa                                                                                          <  1:1256073/150‑1 (MQ=255)                                                         cAATAAACCATGAGGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGat                                                                                         >  2:632436/1‑150 (MQ=255)                                                         cAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGat                                                                                         >  2:59612/1‑150 (MQ=255)                                                            tAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAg                                                                                                                     >  2:594109/1‑119 (MQ=255)                                                            tAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAg                                                                                                                     <  1:594109/119‑1 (MQ=255)                                                              aaCCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGCGCGCGATATGGc                                                                                    >  2:1051788/1‑150 (MQ=255)                                                              aaCCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGc                                                                                    <  2:1299818/150‑1 (MQ=255)                                                               aCCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGACCGCAGCCAGGCTGGGCAGAGCGCGATATGGCa                                                                                   >  1:174142/1‑150 (MQ=255)                                                                ccATGACGCAATACTGGCGGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAg                                                                                  <  2:256282/150‑1 (MQ=255)                                                                 cATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGCGCTGAATGACCGCGAAGGCCGGGATGCATGCG‑TGCAGCGCCTGCCGCAGTGCGCATTTAGCCAGCCGCAAGGCTGGGCAGAGCGCGATATGGCAGa                                                                                 >  1:1088811/1‑150 (MQ=255)                                                                      cGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACtgt                                                                            >  1:450571/1‑150 (MQ=255)                                                                       gCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACtgtg                                                                           <  1:724714/150‑1 (MQ=255)                                                                               ggctggGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAAtgct                                                                   >  1:1162289/1‑150 (MQ=255)                                                                               ggctggGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAAtgct                                                                   >  1:1206424/1‑150 (MQ=255)                                                                                gctggGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAAtgctg                                                                  <  2:116938/150‑1 (MQ=255)                                                                                 ctggGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAAtgctgc                                                                 <  2:174142/150‑1 (MQ=255)                                                                                  tggGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCa                                                                <  2:1032553/150‑1 (MQ=255)                                                                                   ggggTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCAt                                                               >  2:431360/1‑150 (MQ=255)                                                                                    gggTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATg                                                              >  1:506978/1‑150 (MQ=255)                                                                                     ggTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATgg                                                             >  2:414892/1‑150 (MQ=255)                                                                                      gTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGc                                                            >  1:407727/1‑150 (MQ=255)                                                                                      gTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGc                                                            <  1:132584/150‑1 (MQ=255)                                                                                      gTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGc                                                            >  1:401787/1‑150 (MQ=255)                                                                                         aGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACtgtg                                                                           <  1:1022169/132‑1 (MQ=255)                                                                                         aGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACtgtg                                                                           >  2:1022169/1‑132 (MQ=255)                                                                                          gATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTc                                                        >  1:1444725/1‑150 (MQ=255)                                                                                             tgtgGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCAtt                                                     >  1:1303259/1‑150 (MQ=255)                                                                                               tgGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTcg                                                   >  1:586978/1‑150 (MQ=255)                                                                                                ggCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTcgc                                                  >  1:214212/1‑150 (MQ=255)                                                                                                 gCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTcgcg                                                 <  2:1293078/150‑1 (MQ=255)                                                                                                    gTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGccc                                              <  1:510291/150‑1 (MQ=255)                                                                                                    gTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGccc                                              <  1:59612/150‑1 (MQ=255)                                                                                                       gctggctgAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGc                                           <  1:808572/150‑1 (MQ=255)                                                                                                           gctgAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGCGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTAc                                       >  2:624246/1‑150 (MQ=255)                                                                                                            ctgAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACt                                      >  1:858869/1‑150 (MQ=255)                                                                                                              gAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCa                                    >  1:1500805/1‑150 (MQ=255)                                                                                                              gAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCa                                    <  2:549479/150‑1 (MQ=255)                                                                                                                aaGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCAt                                                               >  2:1031552/1‑121 (MQ=255)                                                                                                                aaGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCAt                                                               <  1:1031552/121‑1 (MQ=255)                                                                                                                     gggCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGgt                             <  1:1303704/150‑1 (MQ=255)                                                                                                                      ggCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGgtg                            <  1:859403/150‑1 (MQ=255)                                                                                                                       gCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTTCGCATTAAGCCCGCAGCCAGTGTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGgtgt                           <  2:983147/150‑1 (MQ=255)                                                                                                                       gCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGgtgt                           <  1:700911/150‑1 (MQ=255)                                                                                                                         tGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTa                         >  2:1213939/1‑150 (MQ=255)                                                                                                                          gTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGc                                                                                    >  2:803076/1‑90 (MQ=255)                                                                                                                          gTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGc                                                                                    <  1:803076/90‑1 (MQ=255)                                                                                                                            aTGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCg                      >  2:1495865/1‑150 (MQ=255)                                                                                                                             tGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCgg                     >  2:1013547/1‑150 (MQ=255)                                                                                                                              gACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCggg                    <  1:181952/150‑1 (MQ=255)                                                                                                                               aCCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCgggg                   <  2:214212/150‑1 (MQ=255)                                                                                                                                ccGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGa                  >  1:687247/1‑150 (MQ=255)                                                                                                                                ccGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGa                  >  1:681857/1‑150 (MQ=255)                                                                                                                                  gcgAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAAc                >  2:738144/1‑150 (MQ=255)                                                                                                                                     aaGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACaaa             >  1:1538052/1‑150 (MQ=255)                                                                                                                                      aGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAAt            >  1:1113952/1‑150 (MQ=255)                                                                                                                                       ggCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAAtt           >  2:502793/1‑150 (MQ=255)                                                                                                                                         ccGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGt         >  1:394691/1‑150 (MQ=255)                                                                                                                                          cGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTc        <  2:596117/150‑1 (MQ=255)                                                                                                                                           gggAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCg       <  2:1467631/150‑1 (MQ=255)                                                                                                                                           cagcggcaggcGCTGCA‑CGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCg       >  1:1145996/5‑150 (MQ=255)                                                                                                                                            ggAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCgg      <  2:1206424/150‑1 (MQ=255)                                                                                                                                               ggcaggcGCTGCA‑CGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTAt   >  1:599875/1‑150 (MQ=255)                                                                                                                                               ggcaggcGCTGCA‑CGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTAt   >  1:983921/1‑150 (MQ=255)                                                                                                                                               ggcaggcGCTGCA‑CGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTAt   >  1:594304/1‑150 (MQ=255)                                                                                                                                                gcaggcG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCCGCCAGGCTGGGCAGAGCGCGATATGGCAGCACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATc  >  1:817758/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                TGTTTGTTAAAGATGTGCCGCAGGGAGCTGATCTGCTGAAAGCCGTGGATGAGGTGCAATAAACCATGACGCAATACTGGCTGGGGTTAGATTGTGGCGGTAGCTGGCTGAAAGCCGGGCTGTATGACCGCGAAGGCCGGGAGGCAGGCG‑TGCAGCGCCTGCCGCTGTGCGCATTAAGCCCGCAGCCAGGCTGGGCAGAGCGCGATATGGCAGAACTGTGGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATC  >  NC_000913/3747019‑3747311

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