Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,209,135	A→G	5.1%	pseudogene (36/2001 nt)	yehQ →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,209,135	0	A	G	5.1%	238.0 / 2.8	78	pseudogene (36/2001 nt)	yehQ	pseudogene
Reads supporting (aligned to +/- strand):  ref base A (39/35);  new base G (2/2);  total (41/37)
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.01e-01

CTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTGGCGAACG  >  NC_000913/2208994‑2209283                                                                                                                                              |                                                                                                                                                     cTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAgg                                                                                                                                               >  7:271872/1‑149 (MQ=255)         gACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACg                                                                                                                                       >  4:286954/1‑149 (MQ=255)          aCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACgg                                                                                                                                      >  4:68298/1‑149 (MQ=255)             gcgATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGt                                                                                                                                   <  8:113421/149‑1 (MQ=255)                aTACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTaa                                                                                                                                <  7:55983/149‑1 (MQ=255)                 tACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAg                                                                                                                               <  6:40888/149‑1 (MQ=255)                      cccAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCt                                                                                                                                           <  6:14702/132‑1 (MQ=255)                        cAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGcc                                                                                                                        >  8:24049/1‑149 (MQ=255)                         aGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCg                                                                                                                       <  6:357972/149‑1 (MQ=255)                         aGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCg                                                                                                                       <  7:186278/149‑1 (MQ=255)                         aGGCGCTGGTTAATGTAGGCGCACAAATAGCCGCAATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAAGTATTATAACTGACACCGCAGGCCCGGAAGGCCTTAAGCACTGCCg                                                                                                                       >  4:127797/1‑149 (MQ=255)                            cgcTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGtt                                                                                                                    <  7:273301/149‑1 (MQ=255)                              cTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGtttt                                                                                                                  <  4:323404/149‑1 (MQ=255)                                 gTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGt                                                                                                                >  7:24405/1‑148 (MQ=255)                                        tCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAgcgc                                                                                                         <  3:274504/148‑1 (MQ=255)                                        tCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCa                                                                                                        <  7:89769/149‑1 (MQ=255)                                        tCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGCTTTGTTAAgcg                                                                                                          >  6:297696/1‑147 (MQ=255)                                             gcACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCt                                                                                                    >  2:309435/1‑148 (MQ=255)                                              cacaAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCtt                                                                                                   >  8:123261/1‑148 (MQ=255)                                                caAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAg                                                                                                <  1:309435/149‑1 (MQ=255)                                                 aaaTAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCa                                                                                                        <  4:8092/140‑1 (MQ=255)                                                 aaaTAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAg                                                                                                >  1:243740/1‑148 (MQ=255)                                                 aaaTAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAg                                                                                                >  8:391519/1‑148 (MQ=255)                                                 aaaTAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGACATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCa                                                                                                        >  3:8092/1‑140 (MQ=255)                                                   aTAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGaa                                                                                             >  3:67123/1‑149 (MQ=255)                                                       ccgccATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTgg                                                                                         <  3:68298/149‑1 (MQ=255)                                                           cATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGaaaa                                                                                     >  8:211156/1‑149 (MQ=255)                                                               aCGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATgg                                                                                  <  7:185961/148‑1 (MQ=255)                                                                  ccGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAAc                                                                              <  7:211156/149‑1 (MQ=255)                                                                   cGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCCCTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCCACg                                                                             >  4:274565/1‑149 (MQ=255)                                                                        gAGCTGGCAGCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCg                                                                        <  1:253749/149‑1 (MQ=255)                                                                              gCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAAc                                                                              >  6:109260/1‑137 (MQ=255)                                                                              gCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAAc                                                                              <  5:109260/137‑1 (MQ=255)                                                                               catcatGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCa                                                                 <  8:271872/149‑1 (MQ=255)                                                                                atcatGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAg                                                                >  5:50661/1‑149 (MQ=255)                                                                                    tGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCCGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCAt                                                            <  8:310231/149‑1 (MQ=255)                                                                                    tGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCAt                                                            <  7:391519/149‑1 (MQ=255)                                                                                      gCTTGCGGAGAATCTTCATTCATGAATTCACTACGTCGGGAATTATTAGAACTGACAGCGCAGGCCCTGACGGCGTTAAGCATTGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGTAAAATGGCAACGTCCCGGAGATCAGCCATga                                                          <  5:243534/149‑1 (MQ=255)                                                                                      gCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATg                                                           <  7:401501/148‑1 (MQ=255)                                                                                            ggAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAacga                                                    >  6:353456/1‑149 (MQ=255)                                                                                            ggAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAacga                                                    >  5:355154/1‑149 (MQ=255)                                                                                             gagaATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAacgac                                                   >  4:36769/1‑149 (MQ=255)                                                                                                aaTCGTCAGTCATGAAATGACAACGTCCGGAAATATTAGAACAGACACCCCAGGCCCAGACGGCCTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCATAAGGAACAGGAAAAAGGCAACGGCCCGGAGATCAGCCAATAGAACGACGca                                                >  2:332263/1‑148 (MQ=255)                                                                                                  tCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAAc                                                                              <  8:343491/117‑1 (MQ=255)                                                                                                  tCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAAc                                                                              >  7:343491/1‑117 (MQ=255)                                                                                                     tCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTaa                                                                                                 <  6:358690/95‑1 (MQ=255)                                                                                                     tCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTaa                                                                                                 >  5:358690/1‑95 (MQ=255)                                                                                                      cAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTaa                                            <  1:332263/147‑1 (MQ=255)                                                                                                       aGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATc                                          <  4:67123/148‑1 (MQ=255)                                                                                                       aGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATc                                          <  8:24405/148‑1 (MQ=255)                                                                                                         tCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGc                                        >  1:34686/1‑148 (MQ=255)                                                                                                           aTGAATTCACTACGTCCGGAATTATTAGAACTGACGCCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCAc                                     >  5:22644/1‑149 (MQ=255)                                                                                                           aTGAATTCACTACGTCCGGAATTATTAGAACTGACGCCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCAc                                     >  5:17612/1‑149 (MQ=255)                                                                                                           aTGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCAc                                     >  1:388943/1‑149 (MQ=255)                                                                                                           aTGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCAc                                     >  1:187021/1‑149 (MQ=255)                                                                                                           aTGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCAc                                     >  5:352413/1‑149 (MQ=255)                                                                                                           aTGAATTCACTACCTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAACCAATGCCGGTTTTGTTAAGCGCAGACTTAAGGAACTGGAAAATGGCAACGTCCCCGAGATCAGCCATGAGAACGACGCTTTAATCGCCAc                                     >  2:19993/1‑149 (MQ=255)                                                                                                             gAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCt                                   >  2:68158/1‑149 (MQ=255)                                                                                                              aaTTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCtt                                  <  3:286954/149‑1 (MQ=255)                                                                                                                ttCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGGTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCGATGAGAACGACGCTTTAATCTCCACCtt                                  >  2:258791/1‑147 (MQ=255)                                                                                                                   aCTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCa                                                                 >  3:260002/1‑113 (MQ=255)                                                                                                                   aCTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCa                                                                 <  4:260002/113‑1 (MQ=255)                                                                                                                   aCTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGt                              >  3:377289/1‑148 (MQ=255)                                                                                                                   aCTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGt                              >  7:28985/1‑148 (MQ=255)                                                                                                                    cTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGCTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGa                            <  5:297696/149‑1 (MQ=255)                                                                                                                    cTACGGCCGGAATTATTATAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGa                            <  4:269834/149‑1 (MQ=255)                                                                                                                       cGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACgg                         >  6:246366/1‑149 (MQ=255)                                                                                                                        gTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCttt                                              <  4:173307/127‑1 (MQ=255)                                                                                                                        gTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCttt                                              >  3:173307/1‑127 (MQ=255)                                                                                                                          ccGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGgtgt                      >  3:379201/1‑149 (MQ=255)                                                                                                                           cGGAATTATTAGAACTGACGCCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGgtgt                      <  6:22644/148‑1 (MQ=255)                                                                                                                           cGGAATTATTAGAACTGACGCCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGgtgt                      <  6:17612/148‑1 (MQ=255)                                                                                                                           cGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGgtgt                      <  7:24049/148‑1 (MQ=255)                                                                                                                           cGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAAAGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTc                     >  4:179091/1‑149 (MQ=255)                                                                                                                                   ttaGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCa                                <  7:173519/130‑1 (MQ=255)                                                                                                                                   ttaGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCa                                >  8:173519/1‑130 (MQ=255)                                                                                                                                   ttaGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTAc                >  3:150365/1‑146 (MQ=255)                                                                                                                                   ttaGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTAc                <  4:150365/146‑1 (MQ=255)                                                                                                                                            acacCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTGGCGaa    >  3:189634/1‑149 (MQ=255)                                                                                                                                             cacCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTGGCGAAc   <  5:353456/149‑1 (MQ=255)                                                                                                                                              acCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCATTCTTTAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTGGCGATCg  <  8:128182/149‑1 (MQ=255)                                                                                                                                              acCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTTGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCAGTACCCAGCTGGCGAACg  >  7:313643/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                     CTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTGGCGAACG  >  NC_000913/2208994‑2209283

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