Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,579,154	A→C	9.3%	A163A (GCA→GCC)	talA →	transaldolase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,579,154	0	A	C	9.3%	262.5 / 19.7	97	A163A (GCA→GCC)	talA	transaldolase A
Reads supporting (aligned to +/- strand):  ref base A (50/38);  new base C (4/5);  total (54/43)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

CTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGACTACTATA  >  NC_000913/2579005‑2579302                                                                                                                                                      |                                                                                                                                                     cTGGTGGACTTGTATCAGCAACAATGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCa                                                                                                                                                      >  2:1207903/1‑150 (MQ=255)   ggtggACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACg                                                                                                                                                    >  2:1175784/1‑150 (MQ=255)       gACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGcc                                                                                                                                                <  1:1082599/150‑1 (MQ=255)       gACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGcc                                                                                                                                                <  1:967255/150‑1 (MQ=255)         cTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCtg                                                                                                                                              >  2:348711/1‑150 (MQ=255)              aTCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTTCACAGGCACGGGCCTGTGCgg                                                                                                                                         >  1:651446/1‑150 (MQ=255)              aTCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaa                                                                                                                                         <  2:1340096/150‑3 (MQ=255)              aTCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCgg                                                                                                                                         >  1:644624/1‑150 (MQ=255)              aTCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCgg                                                                                                                                         >  1:597618/1‑150 (MQ=255)              aTCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCgg                                                                                                                                         >  2:338727/1‑150 (MQ=255)               tCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaa                                                                                                                                        <  1:1482055/150‑4 (MQ=255)               tCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaa                                                                                                                                        <  1:410581/150‑4 (MQ=255)                 aGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaaag                                                                                                                                      <  1:1252444/150‑6 (MQ=255)                 aGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGATGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAg                                                                                                                                      >  2:223848/1‑150 (MQ=255)                     acaaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAgg                                                                                                                                  <  2:1418439/150‑1 (MQ=255)                      caaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGc                                                                                                                                 >  2:1482835/1‑150 (MQ=255)                       aaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCg                                                                                                                                <  2:223800/150‑1 (MQ=255)                           cGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGtttt                                                                                                                            >  1:808361/1‑150 (MQ=255)                              tGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCt                                                                                                                         >  2:269060/1‑150 (MQ=255)                              tGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCt                                                                                                                         >  2:1218578/1‑150 (MQ=255)                               gagaAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTg                                                                                                                        <  1:1175784/150‑1 (MQ=255)                               gagaAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTg                                                                                                                        <  2:187228/150‑1 (MQ=255)                                 gaAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGAt                                                                                                                      >  2:1522509/1‑150 (MQ=255)                                  aaaTCACGCATTCTGATCATGCTGGCTTCGACCTGGTAAGGAATTCGCGCGTCAGAAGAGCTGCAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGAtt                                                                                                                     <  1:89435/150‑1 (MQ=255)                                  aaaTCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGAtt                                                                                                                     <  2:482967/150‑1 (MQ=255)                                      cACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGTAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGc                                                                                                                 <  1:1256236/150‑1 (MQ=255)                                         gCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGt                                                                                                              <  1:1182376/150‑1 (MQ=255)                                            ttCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTg                                                                                                           >  1:1443526/1‑150 (MQ=255)                                              cTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTc                                                                                                         >  1:481802/1‑150 (MQ=255)                                                gATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCgg                                                                                                       >  1:266692/1‑150 (MQ=255)                                                   cAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCg                                                                                                    >  2:1354425/1‑150 (MQ=255)                                                   cAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCg                                                                                                    <  2:1187589/150‑1 (MQ=255)                                                    aaGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGt                                                                                                   >  2:1438071/1‑150 (MQ=255)                                                    aaGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGAGTTTTTCTGATTTCGCCGTTTGTCGGGCGt                                                                                                   >  2:1432366/1‑150 (MQ=255)                                                        tGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTAttt                                                                                               <  1:633338/150‑1 (MQ=255)                                                         ggCTTCGACCTGGGAAGGAATTCGCGCGGCAGAACAGCTGGAAAAAGAACGTATTAACTGCAACCTGAGGCTGCTGTTTTCTTTTGGACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTa                                                                                              <  2:678883/150‑1 (MQ=255)                                                           cTTCGACCTGGGAAGGAATTCGCTGGGCAGAAGAGCTTTAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGTAAGCAGGCGTTTTTCTGATTTCGGCGTTTGTCGGGCGTATTTATg                                                                                            <  2:608864/150‑1 (MQ=255)                                                           cTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATg                                                                                            >  2:1580439/1‑150 (MQ=255)                                                           cTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTAGTTATg                                                                                            >  1:1011320/1‑150 (MQ=255)                                                               gACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTg                                                                                        <  1:1207903/150‑1 (MQ=255)                                                                aCCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTgg                                                                                       <  1:112509/150‑1 (MQ=255)                                                                  cTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTa                                                                                     >  1:210004/1‑150 (MQ=255)                                                                     ggaaggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCa                                                                                  <  1:670035/150‑1 (MQ=255)                                                                          gaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCAc                                                                             >  1:125029/1‑150 (MQ=255)                                                                            aTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGc                                                                           <  1:1204703/150‑1 (MQ=255)                                                                                     gCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATg                                                                  <  2:997089/150‑1 (MQ=255)                                                                                     gCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATg                                                                  <  1:150009/150‑1 (MQ=255)                                                                                       agaagaGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGa                                                                >  1:496399/1‑150 (MQ=255)                                                                                        gaagaGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCAGCTGTTTTATTTTTCACAGGCACGGCCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTGTGTCGGGCGTATTTATGACTGGTATCCGGCCCGCAAGCCGATGGaa                                                               >  1:156878/1‑149 (MQ=255)                                                                                           gagCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCg                                                            <  2:1340790/150‑1 (MQ=255)                                                                                           gagCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCg                                                            <  2:1335249/150‑1 (MQ=255)                                                                                            agCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGt                                                           >  1:1396160/1‑150 (MQ=255)                                                                                             gCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTa                                                          <  1:894373/150‑1 (MQ=255)                                                                                             gCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTa                                                          >  2:1478660/1‑150 (MQ=255)                                                                                             gCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTa                                                          <  1:348711/150‑1 (MQ=255)                                                                                             gCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTa                                                          >  2:1185675/1‑150 (MQ=255)                                                                                               tGGAAACAGATGGTATTAACTGCAACCTGTCGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGGAGGCGGTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACTCAAGCCGATGGACCCGTAtg                                                        <  2:1445216/150‑1 (MQ=255)                                                                                                ggAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTAtgt                                                       >  2:1119465/1‑150 (MQ=255)                                                                                                ggAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGGTTTTCTGATTTCTCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTAtgt                                                       >  2:74981/1‑150 (MQ=255)                                                                                                   aaaaGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATgtggt                                                    >  2:939614/1‑150 (MQ=255)                                                                                                   aaaaGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATgtggt                                                    <  1:338727/150‑1 (MQ=255)                                                                                                    aaaGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATgtggtg                                                   <  1:1354425/150‑1 (MQ=255)                                                                                                     aagaagGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGtggtgg                                                  <  2:937940/150‑1 (MQ=255)                                                                                                      agaagGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGga                                                 <  1:1218578/150‑1 (MQ=255)                                                                                                       gaagGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaa                                                <  2:1202424/150‑1 (MQ=255)                                                                                                          ggTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaagaa                                             >  2:1495117/1‑150 (MQ=255)                                                                                                          ggTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaagaa                                             >  2:70300/1‑150 (MQ=255)                                                                                                           gTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATg                                                                  <  2:837047/128‑1 (MQ=255)                                                                                                           gTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATg                                                                  >  1:837047/1‑128 (MQ=255)                                                                                                            tATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAagaaga                                           <  2:81640/150‑1 (MQ=255)                                                                                                            tATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAagaaga                                           >  1:1470601/1‑150 (MQ=255)                                                                                                             aTTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGAt                                          >  1:154141/1‑150 (MQ=255)                                                                                                                    gCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCg                                   >  1:1013941/1‑150 (MQ=255)                                                                                                                       aCCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTa                                <  1:1151892/150‑1 (MQ=255)                                                                                                                       aCCTGACGCTGCTGGTTTGTTTTGCACAGGCACGGGCCTGTGCGGAAGGAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTa                                <  1:340301/150‑1 (MQ=255)                                                                                                                        ccTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAACCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTaa                               >  2:1136902/1‑150 (MQ=255)                                                                                                                             ctctgctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAATCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCTTATGTGGTGGAAGAAGATCCGGGCGTTAAATCgg                          <  1:702442/148‑1 (MQ=255)                                                                                                                               ctgctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTg                        >  2:794564/1‑150 (MQ=255)                                                                                                                                 gctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTgcg                      >  2:1615089/1‑150 (MQ=255)                                                                                                                                 gctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTgcg                      >  2:1004374/1‑150 (MQ=255)                                                                                                                                   tgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCa                    <  1:1457940/150‑1 (MQ=255)                                                                                                                                    gTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCaa                   >  1:68007/1‑150 (MQ=255)                                                                                                                                    gTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCaa                   >  1:780540/1‑150 (MQ=255)                                                                                                                                     ttttCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGTGTTAAATCGGTGCGCAat                  >  2:1580019/1‑150 (MQ=255)                                                                                                                                      tttCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAata                 <  2:644624/150‑1 (MQ=255)                                                                                                                                       tgcttccgCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAatat                >  1:733995/8‑150 (MQ=255)                                                                                                                                         cttccgCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCt              >  1:466636/6‑150 (MQ=255)                                                                                                                                         cTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCt              >  1:572718/1‑150 (MQ=255)                                                                                                                                         cTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCt              >  1:1098528/1‑150 (MQ=255)                                                                                                                                         cTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCt              >  1:382345/1‑150 (MQ=255)                                                                                                                                           tttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTAc            <  2:266692/150‑1 (MQ=255)                                                                                                                                             tGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACgc          >  2:1212869/1‑149 (MQ=255)                                                                                                                                             tGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGa          >  2:1422969/1‑150 (MQ=255)                                                                                                                                             cgCACAGGCCCGTGACTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATAAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCTGGCGTTAAATCGGTGCGAAAAATCTACGa          >  2:666065/2‑150 (MQ=255)                                                                                                                                              gCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGac         >  1:63105/1‑150 (MQ=255)                                                                                                                                              gCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACgcc         >  2:453174/1‑148 (MQ=255)                                                                                                                                                    gCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGACTACtat   >  1:1616591/1‑150 (MQ=255)                                                                                                                                                     cACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGACTACtata  <  1:1006102/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGACTACTATA  >  NC_000913/2579005‑2579302

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