Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,983,573	C→T	7.0%	L184L (CTG→CTA) ‡	kduI ←	hexuronate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,983,573	0	C	T	7.0%	270.5 / 8.5	86	L184L (CTG→CTA) ‡	kduI	hexuronate isomerase
Reads supporting (aligned to +/- strand):  ref base C (37/43);  new base T (3/3);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

GCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGATACGTC  >  NC_000913/2983431‑2983720                                                                                                                                               |                                                                                                                                                    gCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTcc                                                                                                                                               <  3:856750/149‑1 (MQ=255)   aTCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTcc                                                                                                                                               <  3:368520/147‑1 (MQ=255)   aTCACAATATGACGCGTTTCTTGCGGCTGCCCCATAATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTcc                                                                                                                                               >  4:368520/1‑147 (MQ=255)    tCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTc                                                                                                                                            >  3:17363/1‑149 (MQ=255)      acaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCTGGCTCCgggta                                                                                                                                          <  4:493049/149‑5 (MQ=255)      acaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAg                                                                                                                                          >  3:323944/1‑149 (MQ=255)        aaTATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCTGGCTCCgggt                                                                                                                                           >  2:757446/1‑143 (MQ=255)        aaTATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCTGGCTCCgggt                                                                                                                                           <  1:757446/146‑4 (MQ=255)        aaTATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGcc                                                                                                                                        >  2:323704/1‑149 (MQ=255)          tatGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGcccc                                                                                                                                      <  2:635513/149‑1 (MQ=255)            tGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTTGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCAt                                                                                                                                    <  4:405652/149‑1 (MQ=255)                  gTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCaa                                                                                                                              <  2:387976/149‑1 (MQ=255)                  gTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCaa                                                                                                                              >  2:38175/1‑149 (MQ=255)                  gTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCaa                                                                                                                              >  1:561220/1‑149 (MQ=255)                   tttCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAAt                                                                                                                             <  4:323944/149‑1 (MQ=255)                       ttGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGTGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCCGCCCCATACCCAATTGGc                                                                                                                         <  4:457649/149‑1 (MQ=255)                            gCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGtt                                                                                                                    >  3:25160/1‑149 (MQ=255)                             cTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGttt                                                                                                                   >  1:53264/1‑149 (MQ=255)                              tGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTc                                                                                                                  >  3:860563/1‑149 (MQ=255)                              tGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTc                                                                                                                  <  1:367478/149‑1 (MQ=255)                                ccccATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCa                                                                                                                <  4:193473/149‑1 (MQ=255)                                ccccATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCa                                                                                                                <  3:326048/149‑1 (MQ=255)                                 cccATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAg                                                                                                               <  3:542989/149‑1 (MQ=255)                                 cccATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAg                                                                                                               >  4:289100/1‑149 (MQ=255)                                  ccATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGt                                                                                                              >  3:91786/1‑149 (MQ=255)                                   catcatGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTa                                                                                                             >  2:17688/1‑149 (MQ=255)                                      catGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACAt                                                                                                          <  2:53264/149‑1 (MQ=255)                                        tgtgGAAATCGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATcc                                                                                                        >  3:815014/1‑149 (MQ=255)                                          tgGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCgg                                                                                                      >  2:609842/1‑149 (MQ=255)                                             aaaaCGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGAc                                                                                                   >  3:545551/1‑149 (MQ=255)                                                   cAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAata                                                                                             <  4:545551/149‑1 (MQ=255)                                                    aGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAatat                                                                                            >  4:691461/1‑149 (MQ=255)                                                    aGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAatat                                                                                            >  3:56935/1‑149 (MQ=255)                                                        gTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACaaaa                                                                                               <  4:66180/142‑1 (MQ=255)                                                        gTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACaaaa                                                                                               >  3:66180/1‑142 (MQ=255)                                                        gTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGt                                                                                        >  1:272062/1‑149 (MQ=255)                                                             catcCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATc                                                                                   <  4:636613/149‑1 (MQ=255)                                                              atcCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCg                                                                                  >  2:109543/1‑149 (MQ=255)                                                                ccATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTg                                                                                <  3:330665/149‑1 (MQ=255)                                                                 cATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACA‑TGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCg                                                                              >  4:789653/1‑149 (MQ=255)                                                                   tatTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGa                                                                             >  4:319405/1‑149 (MQ=255)                                                                     ttGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACg                                                                           >  3:239919/1‑149 (MQ=255)                                                                        aaatagaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGtt                                                                        <  2:173037/149‑1 (MQ=255)                                                                          atagaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTAc                                                                      >  1:754586/1‑149 (MQ=255)                                                                          atagaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTAc                                                                      <  2:147934/149‑1 (MQ=255)                                                                            agaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCGGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTg                                                                    >  4:559955/1‑149 (MQ=255)                                                                            agaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTg                                                                    <  2:802219/149‑1 (MQ=255)                                                                             gaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTgg                                                                   <  4:815014/149‑1 (MQ=255)                                                                             gaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGATTACTgg                                                                   <  3:601858/149‑1 (MQ=255)                                                                              gaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTACGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCATGCGACGGTTACTGGt                                                                  >  4:413380/2‑149 (MQ=255)                                                                                ataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGa                                                                <  1:561954/149‑1 (MQ=255)                                                                                 taAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAg                                                               <  3:266366/149‑1 (MQ=255)                                                                                 taAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAg                                                               <  4:17363/149‑1 (MQ=255)                                                                                 taAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAg                                                               <  4:71401/149‑1 (MQ=255)                                                                                      ttCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTAt                                                          >  4:349234/1‑149 (MQ=255)                                                                                      ttCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTAt                                                          >  2:896909/1‑149 (MQ=255)                                                                                           tCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCt                                                     >  3:515949/1‑149 (MQ=255)                                                                                                cgcTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGt                                                <  3:109461/149‑1 (MQ=255)                                                                                                 gcTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACgag                                               <  1:209036/149‑3 (MQ=255)                                                                                                 gcTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTg                                               >  4:240869/1‑149 (MQ=255)                                                                                                  cTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGa                                              <  2:66893/149‑1 (MQ=255)                                                                                                   tCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGAc                                             >  4:152661/1‑149 (MQ=255)                                                                                                   tCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGAc                                             <  2:104492/149‑1 (MQ=255)                                                                                                     gTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTg                                           >  1:504987/1‑149 (MQ=255)                                                                                                     gTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTg                                           <  1:17688/149‑1 (MQ=255)                                                                                                        ggTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGgag                                        >  4:156961/1‑149 (MQ=255)                                                                                                          tgtgACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGAt                                      <  1:171981/149‑1 (MQ=255)                                                                                                               cacGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTc                                 <  3:240869/149‑1 (MQ=255)                                                                                                               acatgCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTc                                 <  3:789653/145‑1 (MQ=255)                                                                                                                   gCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAAGTGGCAGGTTTCCAGTACATCCGGGGCAAAATATTTGTTAATCGTGCTACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTcc                             >  2:697717/1‑149 (MQ=255)                                                                                                                      tGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGgt                          <  4:91786/149‑1 (MQ=255)                                                                                                                        gtgtTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGgtgt                        <  1:109543/149‑1 (MQ=255)                                                                                                                              cacaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTgg                                                                   >  2:10867/1‑100 (MQ=255)                                                                                                                              cacaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTgg                                                                   <  1:10867/100‑1 (MQ=255)                                                                                                                              cacaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTg                                                                 <  4:133834/102‑1 (MQ=255)                                                                                                                              cacaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTg                                                                 >  3:133834/1‑102 (MQ=255)                                                                                                                              cacaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACttt                  >  3:484794/1‑149 (MQ=255)                                                                                                                                      tgacGGAGCTGGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTtggtgg          >  3:700145/4‑149 (MQ=255)                                                                                                                                      tgacGGAGCTGGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTGTTTtggtgg          >  1:15710/4‑149 (MQ=255)                                                                                                                                      tgacGGAGCTGGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTAGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTtggtgg          <  4:700145/146‑1 (MQ=255)                                                                                                                                      aCCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTtggtgg          <  3:319405/149‑1 (MQ=255)                                                                                                                                      aCCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTtggtgg          <  1:323704/149‑1 (MQ=255)                                                                                                                                       cccGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTggg         <  4:432745/149‑1 (MQ=255)                                                                                                                                       cccGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTggg         <  3:321762/149‑1 (MQ=255)                                                                                                                                       cccGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTggg         <  1:327322/149‑1 (MQ=255)                                                                                                                                           gAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGATAc     <  4:141510/149‑1 (MQ=255)                                                                                                                                              ccAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGATAcgtc  <  1:362133/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    GCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGATACGTC  >  NC_000913/2983431‑2983720

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑