Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,256,722	A→C	E29A (GAA→GCA)	dgkA →	diacylglycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,256,722	0	A	C	99.0%	354.9 / ‑4.0	98	E29A (GAA→GCA)	dgkA	diacylglycerol kinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (45/52);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GAATCATCGGATTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGAAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCAATAGCGCCATCGAAG  >  NC_000913/4256577‑4256868                                                                                                                                                  |                                                                                                                                                   gAATCATCGGATTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCg                                                                                                                                                >  1:511091/1‑150 (MQ=255)  aaTCATCGGATTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGc                                                                                                                                              >  1:482448/1‑151 (MQ=255)     catcGGATTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCAtt                                                                                                                                           <  1:593963/151‑1 (MQ=255)     catcGGATTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCAtt                                                                                                                                           >  1:1288169/1‑151 (MQ=255)     catcGGATTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCAt                                                                                                                                            <  1:32823/150‑1 (MQ=255)      atcGGACTTAACGGTACACTGATATTGACGCTCATCATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAATCTGCTGCCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCACCGGCATTc                                                                                                                                          >  1:1178601/1‑151 (MQ=255)        cGGATTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTcc                                                                                                                                         <  1:766279/150‑1 (MQ=255)        cGGATTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCg                                                                                                                                        <  2:376975/151‑1 (MQ=255)        cGGATTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCg                                                                                                                                        >  2:169152/1‑151 (MQ=255)          gATTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTc                                                                                                                                      >  2:967196/1‑151 (MQ=255)          gATTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTc                                                                                                                                      <  1:1025101/151‑1 (MQ=255)           aTTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCCCCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCCTCGGCATTCCGTc                                                                                                                                      <  1:1299573/150‑1 (MQ=255)           aTTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCa                                                                                                                                     <  2:973174/151‑1 (MQ=255)                       acTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGc                                                                                                                         <  1:779465/151‑1 (MQ=255)                        cTGATATTGACGCTCATAATGTAAACAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAACGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCg                                                                                                                        >  1:1130191/1‑151 (MQ=255)                         tGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGACCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCgg                                                                                                                       >  1:265347/1‑151 (MQ=255)                           atatTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGt                                                                                                                      >  2:1140552/1‑150 (MQ=255)                            tatTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTa                                                                                                                     <  1:533180/150‑1 (MQ=255)                             atTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTAt                                                                                                                    <  1:561541/150‑1 (MQ=255)                                 aCGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTAttgttg                                                                                                               >  2:194604/1‑151 (MQ=255)                                 aCGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTAttgttg                                                                                                               <  2:591904/151‑1 (MQ=255)                                     tCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCgg                                                                                                           >  1:823476/1‑151 (MQ=255)                                        tAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCggtgg                                                                                                        >  1:902307/1‑151 (MQ=255)                                        tAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCggtgg                                                                                                        >  1:106175/1‑151 (MQ=255)                                         aaTGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCggtggt                                                                                                       <  1:245296/151‑1 (MQ=255)                                         aaTGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCggtgg                                                                                                        <  1:1220070/150‑1 (MQ=255)                                            gTAAAAAGGTTCTTTCAATGGCCAATAATTCCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCa                                                                                                     <  1:358852/150‑1 (MQ=255)                                            gTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGAGTAGCGGTATTGTTGGCGGTGGTCAt                                                                                                    <  1:363595/151‑1 (MQ=255)                                              aaaaaGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTAAGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCg                                                                                                  >  1:367832/1‑151 (MQ=255)                                               aaaaGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGc                                                                                                 >  2:759346/1‑151 (MQ=255)                                               aaaaGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGAAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGc                                                                                                 >  1:547103/1‑151 (MQ=255)                                                aaaGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGcc                                                                                                >  2:700246/1‑151 (MQ=255)                                                        tttCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTgctggct                                                                                        >  2:299596/1‑151 (MQ=255)                                                        tttCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTgctggct                                                                                        >  1:1294656/1‑151 (MQ=255)                                                           cAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGctggctgg                                                                                      >  2:315664/1‑150 (MQ=255)                                                            aaTGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGAt                                                                                    >  2:1007695/1‑151 (MQ=255)                                                             aTGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGAtg                                                                                   >  2:52810/1‑151 (MQ=255)                                                               ggCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGAtgtg                                                                                 <  1:708771/151‑1 (MQ=255)                                                                    ataataCCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGAcgc                                                                            >  2:556379/1‑151 (MQ=255)                                                                    ataataCCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGAcgc                                                                            <  1:583798/151‑1 (MQ=255)                                                                      aataCCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGAcgcg                                                                           <  1:198416/150‑1 (MQ=255)                                                                      aataCCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGCGGTCATCGCCTGCTGGCTGGATGTGGAcgcg                                                                           >  2:473552/1‑150 (MQ=255)                                                                        taCCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGAtt                                                                        >  2:1041514/1‑151 (MQ=255)                                                                           cACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTAcc                                                                     >  1:478902/1‑151 (MQ=255)                                                                                gATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCcgcg                                                                 >  1:525323/1‑150 (MQ=255)                                                                                gATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGt                                                                <  1:808526/151‑1 (MQ=255)                                                                                 aTTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGtg                                                               <  2:75649/151‑1 (MQ=255)                                                                                  ttCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGtgc                                                              >  1:975532/1‑151 (MQ=255)                                                                                    cACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGtgctg                                                            <  2:1062109/151‑1 (MQ=255)                                                                                      cccGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGtgctgc                                                           <  2:860965/150‑1 (MQ=255)                                                                                      cccGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTgctgct                                                          <  1:567251/151‑1 (MQ=255)                                                                                      cccGAATTATCAAAGCTGCTGGCAATTCCTGGAAAGGTTTACGCGCGGCATGGATCACCGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGAATTGTTGGCGGGGGCCATCGCCTGCTGGCTGGAGGTGGACGCGATTACCCGCGGgctgct                                                          >  1:1231880/1‑151 (MQ=255)                                                                                       ccGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCtt                                                         >  1:1192825/1‑151 (MQ=255)                                                                                        cGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTa                                                        <  1:888533/151‑1 (MQ=255)                                                                                        cGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGACATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCtt                                                         >  2:1271189/1‑150 (MQ=255)                                                                                         gAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTAt                                                       >  1:922388/1‑151 (MQ=255)                                                                                          aaTTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATc                                                      >  1:111994/1‑151 (MQ=255)                                                                                             tATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAg                                                    <  1:315278/150‑1 (MQ=255)                                                                                              aTCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAg                                                    >  1:166342/1‑149 (MQ=255)                                                                                              aTCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCt                                                  <  2:271562/151‑1 (MQ=255)                                                                                              aTCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCt                                                  >  2:1212927/1‑151 (MQ=255)                                                                                               tCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCt                                                  >  1:940686/1‑150 (MQ=255)                                                                                                cAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTc                                                 <  1:399394/150‑1 (MQ=255)                                                                                                 aaaGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCg                                               >  1:159374/1‑151 (MQ=255)                                                                                                  aaGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCg                                               >  2:641384/1‑150 (MQ=255)                                                                                                   aGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTg                                             <  2:199915/151‑1 (MQ=255)                                                                                                    gctgctGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTg                                             >  2:328387/1‑150 (MQ=255)                                                                                                      tgctgGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATg                                          >  1:718166/1‑151 (MQ=255)                                                                                                       gctgGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGc                                         <  1:206142/151‑1 (MQ=255)                                                                                                       gctgGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGc                                         >  1:783874/1‑151 (MQ=255)                                                                                                       gctgGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTAGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATg                                          <  1:264202/150‑1 (MQ=255)                                                                                                              aTTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGtgat                                  <  2:463229/151‑1 (MQ=255)                                                                                                              aTTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGtg                                    <  2:654334/149‑1 (MQ=255)                                                                                                              aTTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGtg                                    <  1:1258994/149‑1 (MQ=255)                                                                                                              aTTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGtg                                    <  1:937065/149‑1 (MQ=255)                                                                                                              aTTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGtg                                    <  1:117065/149‑1 (MQ=255)                                                                                                                tCCTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGtgatg                                 <  2:461183/150‑1 (MQ=255)                                                                                                                 ccTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGtgatga                                >  1:1071065/1‑150 (MQ=255)                                                                                                                 ccTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGtgatga                                <  1:592065/150‑1 (MQ=255)                                                                                                                 ccTGGAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTgatgat                               >  1:815717/1‑151 (MQ=255)                                                                                                                     gAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATtgtg                           >  1:130929/1‑151 (MQ=255)                                                                                                                     gAAAGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATtgtg                           <  1:541610/151‑1 (MQ=255)                                                                                                                       aaGGTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGa                         <  2:1231242/151‑1 (MQ=255)                                                                                                                         ggTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGaaa                       <  1:689322/151‑1 (MQ=255)                                                                                                                         ggTTTACGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGaaa                       <  1:937625/151‑1 (MQ=255)                                                                                                                              aCGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCt                   <  1:256007/150‑1 (MQ=255)                                                                                                                              aCGCGCTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTc                  >  1:213487/1‑151 (MQ=255)                                                                                                                                gcgcTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCaa                <  2:482461/151‑1 (MQ=255)                                                                                                                                  gcTGCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCAATa              <  1:135022/151‑1 (MQ=255)                                                                                                                                     gCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCAATAgcg           >  2:172328/1‑151 (MQ=255)                                                                                                                                     gCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCAATAgcg           <  2:771692/151‑1 (MQ=255)                                                                                                                                     gCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCAATAgc            <  1:903324/150‑1 (MQ=255)                                                                                                                                     gCATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCAATAg             <  1:1135191/149‑1 (MQ=255)                                                                                                                                      cATGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCAATAgcgc          <  2:63821/151‑1 (MQ=255)                                                                                                                                       aTGGATCAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCAATAGCGcc         <  2:154137/151‑1 (MQ=255)                                                                                                                                             cAACGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCAATAGCGCCATCGaa   <  1:1189921/151‑1 (MQ=255)                                                                                                                                              aaCGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCAATAGCGCCATCGAAg  <  2:839193/151‑1 (MQ=255)                                                                                                                                              aaCGCAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCAATAGCGCCATCGAAg  <  2:1231955/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GAATCATCGGATTTAACGGTACACTGATATTGACGCTCATAATGTAAAAAGGTTCTTTCAATGGCCAATAATACCACTGGATTCACCCGAATTATCAAAGCTGCTGGCTATTCCTGGAAAGGTTTACGCGCTGCATGGATCAACGAAGCGGCATTCCGTCAGGAAGGCGTAGCGGTATTGTTGGCGGTGGTCATCGCCTGCTGGCTGGATGTGGACGCGATTACCCGCGTGCTGCTTATCAGCTCCGTGATGCTGGTGATGATTGTGGAAATCCTCAATAGCGCCATCGAAG  >  NC_000913/4256577‑4256868

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