Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,217,943	T→C	5.2%	Y279H (TAT→CAT)	aceA →	isocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,217,943	0	T	C	5.2%	235.6 / 2.7	76	Y279H (TAT→CAT)	aceA	isocitrate lyase
Reads supporting (aligned to +/- strand):  ref base T (42/30);  new base C (2/2);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAA  >  NC_000913/4217795‑4218086                                                                                                                                                     |                                                                                                                                                ggTTGAGCGTACCGATGCTGATGCGGCGGATCTTATGACCTCCGATTGCGACTCGTATGACATCGAATTTATTACCGGCGAGCGTACCATTGAAGTCTTGTTGAGTTCTCATGCGTGCATTTATCAAGCGATCAGCCTTGTCCTGGCGTa                                                                                                                                                <  2:133920/150‑1 (MQ=255)  gTTGCCCGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTAt                                                                                                                                               >  1:1007838/1‑150 (MQ=255)       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:872345/150‑1 (MQ=255)                                  gATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAAc                                                                                                               >  2:1043787/1‑150 (MQ=255)                                       ccTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGGGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCa                                                                                                          >  1:347692/1‑150 (MQ=255)                                          ccGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACctggtctgg                                                                                                                       <  1:1551294/134‑1 (MQ=255)                                          ccGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACctggtctgg                                                                                                                       >  2:1551294/1‑134 (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                                                                                     >  2:685210/1‑150 (MQ=255)                                                            aCAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGgcgc                                                                                     >  1:52031/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                                                                            <  1:354247/127‑1 (MQ=255)                                                                                            gAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTg                                                                            >  2:354247/1‑127 (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:1175305/1‑150 (MQ=255)                                                                                                                     gCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACt                            >  1:1117682/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:1299654/3‑150 (MQ=255)                                                                                                                                              tatGGCGCATACGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGaa   >  1:63897/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)                                                                                                                                                     |                                                                                                                                                GGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAA  >  NC_000913/4217795‑4218086

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