Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,061,199	A→G	6.0%	Q12R (CAA→CGA)	yihM →	putative sugar phosphate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,061,199	0	A	G	6.0%	197.9 / 2.4	67	Q12R (CAA→CGA)	yihM	putative sugar phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (25/37);  new base G (2/2);  total (27/40)
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.00e-01

AGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCG  >  NC_000913/4061052‑4061342                                                                                                                                                    |                                                                                                                                                aGCGATCTTAATTGAAACGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGGCGTAGAGCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACaa                                                                                                                                                >  4:44433/1‑149 (MQ=255) aGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACaa                                                                                                                                                >  7:28421/1‑149 (MQ=255) aGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACaa                                                                                                                                                >  5:199738/1‑149 (MQ=255) aGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACaa                                                                                                                                                >  4:84756/1‑149 (MQ=255)  gCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATTACGCAAGAAAGATTCTACAAc                                                                                                                                               <  2:226123/149‑1 (MQ=255)     aTCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACgtg                                                                                                                                            >  8:10691/1‑149 (MQ=255)       cTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTc                                                                                                                                          >  3:1611/1‑149 (MQ=255)        ttAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTc                                                                                                                                          >  4:62279/1‑148 (MQ=255)        ttAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTc                                                                                                                                          >  4:13895/1‑148 (MQ=255)         tAATTGAAAAGCATCGCGCCGACCAGCAAGGCATTATTGTCCAAATAGATATTTATTATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACATGTAATAACGCAAGAAAGATTCGTCAACGTGTCGa                                                                                                                                        >  5:76044/1‑149 (MQ=255)           aTTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGac                                                                                                                                       <  1:161709/148‑1 (MQ=255)            ttGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGacac                                                                                                                                     >  2:89223/1‑149 (MQ=255)            ttGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGacac                                                                                                                                     >  8:13738/1‑149 (MQ=255)               aaaaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACActct                                                                                                                                  >  4:54462/1‑149 (MQ=255)               aaaaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACActct                                                                                                                                  <  4:134300/149‑1 (MQ=255)                aaaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCtt                                                                                                                                 >  8:170408/1‑149 (MQ=255)                 aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTc                                                                                                                                >  8:215935/1‑149 (MQ=255)                 aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTc                                                                                                                                <  3:229844/149‑1 (MQ=255)                    cATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCctc                                                                                                                             <  3:101081/149‑1 (MQ=255)                     aTCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATATCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCctct                                                                                                                            >  1:54512/1‑149 (MQ=255)                     aTCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCctct                                                                                                                            <  7:215935/149‑1 (MQ=255)                        gcgcCGACAAGCACGGCAATATTGTCCAACTAGATATTGAATATTGGCGAGTTGAGGCCGTAGACATCATCATTAATCTGTAGCTGTTTTATGGTGACAATAAATACCGCAAGAAAGATTCTCCAACGTGGCGACACTCTTCCTCttt                                                                                                                          <  5:59534/148‑1 (MQ=255)                         cgcCGACCAGCAACGCAATACTGTCCAAATAGATATTTAATATTGGCGATTTGATGCCGTAGACCTCTTTATTTTTATGTAGTTGTTTTATGGTGACTATAAATAACGTAATTAAGATTGTATATCTTGTCGACACTCTTCCTCtttat                                                                                                                        <  4:31287/149‑1 (MQ=255)                              aCCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACa                                                                                                                   <  2:227910/149‑1 (MQ=255)                                cAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACGACGTGTCGACACTCTTCCTCTTTATTTACATg                                                                                                                 >  7:150309/1‑149 (MQ=255)                                cAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACGACGTGTCGACACTCTTCCTCTTTATTTACATg                                                                                                                 >  7:170535/1‑149 (MQ=255)                                   cAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCt                                                                                                               >  1:62943/1‑148 (MQ=255)                                    aaGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTa                                                                                                             >  2:225070/1‑149 (MQ=255)                                    aaGGCAATATTCTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGGGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTa                                                                                                             <  6:222412/149‑1 (MQ=255)                                     tggCAATATTGTCCAAATAGAAATTGAATATTGCCTACTTGAGGCCGTAGACATCATCATTAATCAGTAGGTGTTTTATGCTGACAATAAATAACGCAAGAAAGATTCTACAACGTCTCGACACTCTTCCTCTTTATTTACATGCTTAt                                                                                                            <  6:76044/148‑1 (MQ=255)                                      ggCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATg                                                                                                           >  2:27814/1‑149 (MQ=255)                                       gCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGc                                                                                                          <  2:1507/149‑1 (MQ=255)                                             ttGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTc                                                                                                     <  8:72609/148‑1 (MQ=255)                                                   aaaTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTaaa                                                                                              >  7:64185/1‑149 (MQ=255)                                                   aaaTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTaaa                                                                                              <  8:141833/149‑1 (MQ=255)                                                         atatTGAATATTGGCGATTTTAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAAAGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTATTTAATTACATGCTTATGcc                                                                                                         >  8:173134/1‑132 (MQ=255)                                                            ttGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGc                                                                                                          <  7:87686/128‑1 (MQ=255)                                                            ttGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGc                                                                                                          >  8:87686/1‑128 (MQ=255)                                                                      gctaTTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAgt                                                                            <  3:13895/145‑1 (MQ=255)                                                                        gATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGtt                                                                         <  1:57111/149‑1 (MQ=255)                                                                        gATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGtt                                                                         <  7:170408/149‑1 (MQ=255)                                                                        gATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGtt                                                                         <  4:1611/149‑1 (MQ=255)                                                                            tGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGcc                                                                      <  2:62943/148‑1 (MQ=255)                                                                            tGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGcc                                                                      <  7:71790/148‑1 (MQ=255)                                                                                 ccGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGGGTCGACACTCTTCCTCTTTATTTACATGCTTAAGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGa                                                                <  4:93141/149‑1 (MQ=255)                                                                                       aCCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACt                                                          <  7:13738/149‑1 (MQ=255)                                                                                        ccTCATCATTAATCTGTAGGTGTTTTATGGTTACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACt                                                          <  3:62279/148‑1 (MQ=255)                                                                                        ccTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACtt                                                         <  6:199738/149‑1 (MQ=255)                                                                                               atTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCg                                                  >  7:171775/1‑149 (MQ=255)                                                                                                ttAATCTGTATGGTTTATATGGTGACAAAAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGc                                                 <  6:91360/149‑1 (MQ=255)                                                                                                ttAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCg                                                  >  8:210775/1‑148 (MQ=255)                                                                                                ttAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATACTACCACGTGGCGACACTCTTCCTCTTTATTTACATGCATATGCCTTTCATTTAAATATGCGGCTGGACAGAGTGTTGCCTGCTGATTTACTTGATATCGc                                                 <  6:119975/149‑1 (MQ=255)                                                                                                   aTCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCaa                                               <  1:225070/148‑1 (MQ=255)                                                                                                        tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGaaa                                         <  3:84756/149‑1 (MQ=255)                                                                                                        tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGaaa                                         <  3:44433/149‑1 (MQ=255)                                                                                                              ttttATGGTGACAATAAATAACGCAAGAAAGATTCTACGACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAaataat                                    <  8:150309/148‑1 (MQ=255)                                                                                                              ttttATGGTGACAATAAATAACGCAAGAAAGATTCTACGACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAaataat                                    <  8:170535/148‑1 (MQ=255)                                                                                                                ttATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTa                                 >  1:119343/1‑149 (MQ=255)                                                                                                                   tGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGt                              <  6:67337/149‑1 (MQ=255)                                                                                                                         cAATAAATAACGCAAGAAACATTCTACAACGTATCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAATTATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTc                        <  2:178155/149‑1 (MQ=255)                                                                                                                            taaataaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCtt                                                                                                       <  5:214010/67‑1 (MQ=255)                                                                                                                            taaataaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCtt                                                                                                       >  6:214010/1‑67 (MQ=255)                                                                                                                              aataaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCaaaa                    <  3:219104/148‑1 (MQ=255)                                                                                                                               ataaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATATACATGCTTATCCCTTTCATTTAAATATGCGGCTGCAAAGAGTCTTGCCTGCTGATATACTTGATATCGCAAGTGAAAATAATCTACTTGGCGTCAAAAt                   >  3:230304/1‑148 (MQ=255)                                                                                                                                  aCGCCAGAAAGATTCTACAACGTGCCGACACTCTTCCTCTTTATTAACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCa                <  8:28421/148‑1 (MQ=255)                                                                                                                                    gCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAaataat                                    <  5:102740/126‑1 (MQ=255)                                                                                                                                    gCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAaataat                                    >  6:102740/1‑126 (MQ=255)                                                                                                                                         aaagaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCATTTCATTTAAATATGCGGCTGGAAAGAGAGTTGACTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTtctc         >  2:96597/1‑148 (MQ=255)                                                                                                                                            gaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATg     <  2:119343/149‑1 (MQ=255)                                                                                                                                            gaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATg     >  6:9209/1‑149 (MQ=255)                                                                                                                                                cTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCg  <  8:64185/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                AGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCG  >  NC_000913/4061052‑4061342

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