Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,217,946	G→A	5.2%	A280T (GCG→ACG)	aceA →	isocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,217,946	0	G	A	5.2%	241.7 / 2.0	76	A280T (GCG→ACG)	aceA	isocitrate lyase
Reads supporting (aligned to +/- strand):  ref base G (43/29);  new base A (2/2);  total (45/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

CCGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGAC  >  NC_000913/4217801‑4218095                                                                                                                                                  |                                                                                                                                                      ccGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGcc                                                                                                                                                   <  1:391062/150‑1 (MQ=255)   ccaCCGATTCTTATCCGTCGGATCTGATTTCCCCCCAATGCGCCCCTTCTGACAGGGAATTTATTCCCTGCCCGCTTTGTTTTGTTTTCATCTTCCTTAGTCCTGCGGGAATTTAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCat                                                                                                                                                 <  2:554768/148‑1 (MQ=255)    tACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCata                                                                                                                                                <  2:909879/150‑1 (MQ=255)    tACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGCATACGCCagg                                                                                                                                                <  2:397060/150‑3 (MQ=255)     aCCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGCATACGCCaggc                                                                                                                                               <  2:911810/150‑4 (MQ=255)      ccGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCCTGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATg                                                                                                                                              >  1:879011/1‑150 (MQ=255)        gATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCt                                                                                                                                            >  2:1592817/1‑150 (MQ=255)             tGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGCTCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACct                                                                                                                                       >  2:1588685/1‑150 (MQ=255)                 gcggcgGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACctggtc                                                                                                                                   >  1:212831/1‑150 (MQ=255)                  cggcggATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACctggtct                                                                                                                                  >  1:105186/1‑150 (MQ=255)                     cggATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACctggtctggt                                                                                                                               >  2:556214/1‑150 (MQ=255)                      ggATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGgtg                                                                                                                              <  2:1360635/150‑1 (MQ=255)                      ggATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGgtg                                                                                                                              <  1:730098/150‑1 (MQ=255)                           tGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGaaa                                                                                                                         >  2:1620082/1‑150 (MQ=255)                            gATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAAc                                                                                                                        >  2:1043787/1‑150 (MQ=255)                            gATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAAc                                                                                                                        <  2:872345/150‑1 (MQ=255)                                 ccTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGGGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCa                                                                                                                   >  1:347692/1‑150 (MQ=255)                                    ccGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACctggtctgg                                                                                                                                >  2:1551294/1‑134 (MQ=255)                                    ccGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACctggtctgg                                                                                                                                <  1:1551294/134‑1 (MQ=255)                                          gCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATc                                                                                                          <  2:1007838/150‑1 (MQ=255)                                              cccGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGa                                                                                                      <  2:212831/150‑1 (MQ=255)                                              cccGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGa                                                                                                      <  1:546093/150‑1 (MQ=255)                                               ccGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGaa                                                                                                     >  1:149174/1‑150 (MQ=255)                                                cGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAAc                                                                                                    <  2:1290144/150‑1 (MQ=255)                                                 gTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACt                                                                                                   <  1:1627748/150‑1 (MQ=255)                                                      aCAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGgcgc                                                                                              >  1:52031/1‑150 (MQ=255)                                                      aCAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGgcgc                                                                                              >  2:685210/1‑150 (MQ=255)                                                       cAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGcgcg                                                                                             <  2:645101/150‑1 (MQ=255)                                                       cAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCAAGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGcgcg                                                                                             >  2:1442150/1‑150 (MQ=255)                                                          cGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCg                                                                                          >  1:1173169/1‑150 (MQ=255)                                                          agAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTAATCATGCGGGCATTGAGCAAGCGATCATCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCg                                                                                          >  2:425381/2‑150 (MQ=255)                                                               ttattaCCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTg                                                                                     <  1:1592817/150‑1 (MQ=255)                                                                tattaCCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCAGATCTGGAACTGGCGCGTCGCTTTGc                                                                                    <  1:77230/150‑1 (MQ=255)                                                                 attaCCGGCGAGCGTACCATTGAAGGCTTCTTACGTACTAATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCAATATTCTGACCTGATCACGAGTGAAACCTCCACGCCGGATCTGGAACTCGCGCGTCCCTTAGca                                                                                   >  2:753048/1‑150 (MQ=255)                                                                 attaCCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGca                                                                                   <  1:1620082/150‑1 (MQ=255)                                                                   taCCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGGGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGcaca                                                                                 <  2:347692/150‑1 (MQ=255)                                                                   taCCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGcaca                                                                                 >  2:1106692/1‑150 (MQ=255)                                                                      cGGCGAGCGTACCAGGGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGc                                                                              <  2:462524/150‑1 (MQ=255)                                                                           aGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATcc                                                                         <  2:149174/150‑1 (MQ=255)                                                                              gTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCAcg                                                                      >  2:1339167/1‑150 (MQ=255)                                                                               tACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCAcgc                                                                     >  1:255411/1‑150 (MQ=255)                                                                                aCCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAATTGGCGCGTCGCTTTGCACAAGCTATCCAcgcg                                                                    >  1:159073/1‑150 (MQ=255)                                                                                  cAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGaa                                                                  >  1:825529/1‑150 (MQ=255)                                                                                  cAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGaa                                                                  >  1:451522/1‑150 (MQ=255)                                                                                   aGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGaaa                                                                 >  1:567115/1‑150 (MQ=255)                                                                                    gTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAat                                                                >  1:1053450/1‑150 (MQ=255)                                                                                      gAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTg                                                                                     >  2:354247/1‑127 (MQ=255)                                                                                      gAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTg                                                                                     <  1:354247/127‑1 (MQ=255)                                                                                         ggCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATGTCCg                                                           >  2:866158/1‑150 (MQ=255)                                                                                          gCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCATATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCgg                                                          <  1:1436362/150‑1 (MQ=255)                                                                                              cttcCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCaa                                                      >  1:413061/1‑150 (MQ=255)                                                                                              cttcCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCaa                                                      >  1:1340483/1‑150 (MQ=255)                                                                                                  cGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAActg                                                  <  2:1053450/150‑1 (MQ=255)                                                                                                   gTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAActgc                                                 <  2:413061/150‑1 (MQ=255)                                                                                                      cTCATGCGGGCATTGAGCAACCGATAAGCCGAGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTgg                                              >  2:1546230/1‑150 (MQ=255)                                                                                                      cTCATGCAGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTgg                                              >  1:1206405/1‑150 (MQ=255)                                                                                                        cATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCt                                            >  2:1160064/1‑150 (MQ=255)                                                                                                        cATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCt                                            <  1:556214/150‑1 (MQ=255)                                                                                                          tGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAAATGATGGATta                                          >  2:1036975/1‑150 (MQ=255)                                                                                                            cGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTtata                                        <  2:305272/150‑1 (MQ=255)                                                                                                               gCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACt                                     >  1:1117682/1‑150 (MQ=255)                                                                                                               gCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACt                                     >  1:1175305/1‑150 (MQ=255)                                                                                                                cATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGca                                                                                   >  2:1532315/1‑103 (MQ=255)                                                                                                                cATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGca                                                                                   <  1:1532315/103‑1 (MQ=255)                                                                                                                 aTTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGc                                   >  1:281734/1‑150 (MQ=255)                                                                                                                          gCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCcgtcg                          <  1:1546379/150‑1 (MQ=255)                                                                                                                           cGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCcgtcgt                         >  2:394356/1‑150 (MQ=255)                                                                                                                           cGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCcgtcgt                         >  2:833320/1‑150 (MQ=255)                                                                                                                                       gCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGa             >  1:1178461/1‑150 (MQ=255)                                                                                                                                        tatGGCGCATACGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGaa            >  1:63897/3‑150 (MQ=255)                                                                                                                                        tatGGCGCATACGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGaa            >  1:1299654/3‑150 (MQ=255)                                                                                                                                        ccTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGaa            >  2:1591926/1‑150 (MQ=255)                                                                                                                                         cTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGaaa           <  2:1003584/150‑1 (MQ=255)                                                                                                                                         cTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGaaa           <  2:1009536/150‑1 (MQ=255)                                                                                                                                                aTGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTcg    >  1:794301/1‑150 (MQ=255)                                                                                                                                                aTGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTcg    >  1:142758/1‑150 (MQ=255)                                                                                                                                                  gcgcCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTcgac  <  2:567115/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CCGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCTCGAC  >  NC_000913/4217801‑4218095

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: ‑