Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,983,574	A→G	5.9%	L184P (CTG→CCG) ‡	kduI ←	hexuronate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,983,574	0	A	G	5.9%	270.3 / 6.0	84	L184P (CTG→CCG) ‡	kduI	hexuronate isomerase
Reads supporting (aligned to +/- strand):  ref base A (34/45);  new base G (3/2);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

GCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGATACGTC  >  NC_000913/2983431‑2983720                                                                                                                                                |                                                                                                                                                   gCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTcc                                                                                                                                               <  4:114917/149‑1 (MQ=255)   aTCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGt                                                                                                                                             >  2:520354/1‑149 (MQ=255)   aTCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGt                                                                                                                                             >  3:151135/1‑149 (MQ=255)     cacaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCTGGCTCCgggt                                                                                                                                           <  3:1879/149‑4 (MQ=255)        aaTATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGcc                                                                                                                                        <  4:526015/149‑1 (MQ=255)         atatGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCATCAGccc                                                                                                                                       <  1:315468/149‑1 (MQ=255)           atGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCa                                                                                                                                     <  2:698391/149‑1 (MQ=255)           atGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCAACCATATTGAAATAGAAATAAACTCCCCTCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCa                                                                                                                                     <  2:741842/149‑1 (MQ=255)             gACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATa                                                                                                                                   >  2:671249/1‑149 (MQ=255)                gcgTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTc                                                                                                                                >  1:561100/1‑149 (MQ=255)                gcgTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTc                                                                                                                                >  1:741345/1‑149 (MQ=255)                 cgTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCa                                                                                                                               >  1:690860/1‑149 (MQ=255)                   tttCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCTGGCTCCgg                                                                                                                                             <  3:529680/133‑2 (MQ=255)                   tttCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCTGGCTCCgg                                                                                                                                             >  4:529680/1‑132 (MQ=255)                   tttCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAAt                                                                                                                             >  2:30206/1‑149 (MQ=255)                    ttCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAAtt                                                                                                                            <  4:149343/149‑1 (MQ=255)                    ttCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAAtt                                                                                                                            <  1:671249/149‑1 (MQ=255)                         gCGGCTGCCGCATCATGTGGAAAACGCATGCGTCATCATCCATATTGAACTAGAAATAAACTTCCATCCGGCGCTTGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAg                                                                                                                       <  3:475375/149‑1 (MQ=255)                                 cccATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACGCGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAg                                                                                                               >  2:488338/1‑149 (MQ=255)                                  ccATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGt                                                                                                              >  1:228971/1‑149 (MQ=255)                                    atcatGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTAc                                                                                                            >  2:337186/1‑149 (MQ=255)                                         gtgGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCg                                                                                                       <  1:386269/149‑1 (MQ=255)                                          tgGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCgg                                                                                                      <  2:228971/149‑1 (MQ=255)                                              aaaCGCAGGCGTCATCATCCATATTGAAATAGAAATACACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACa                                                                                                  >  2:69823/1‑149 (MQ=255)                                              aaaCGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACa                                                                                                  >  4:329054/1‑149 (MQ=255)                                                 cGCAGGCGTCATCATCCATATGGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACaaaa                                                                                               <  2:138778/149‑1 (MQ=255)                                                  gCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAat                                                                                              <  2:561100/149‑1 (MQ=255)                                                      gCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACAGCCGGGACAAAATAttt                                                                                          <  4:67802/149‑1 (MQ=255)                                                        gTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGt                                                                                        <  1:281974/149‑1 (MQ=255)                                                         tcatcaTCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGtt                                                                                       >  4:408568/1‑149 (MQ=255)                                                            tcatcCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAAt                                                                                    <  2:741345/149‑1 (MQ=255)                                                             catcCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATc                                                                                   <  1:318594/149‑1 (MQ=255)                                                             catcCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATc                                                                                   <  3:27050/149‑1 (MQ=255)                                                             catcCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATc                                                                                   <  2:690860/149‑1 (MQ=255)                                                             catcCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATc                                                                                   <  3:757951/149‑1 (MQ=255)                                                              atcCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATAttt                                                                                          >  1:341340/1‑141 (MQ=255)                                                              atcCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATAttt                                                                                          <  2:341340/141‑1 (MQ=255)                                                              atcCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCg                                                                                  >  1:678519/1‑149 (MQ=255)                                                                 cATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGc                                                                               <  3:173030/149‑1 (MQ=255)                                                                    atTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGAc                                                                            <  4:736528/149‑1 (MQ=255)                                                                         aatagaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTa                                                                       <  4:151135/149‑1 (MQ=255)                                                                           tagaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACt                                                                     <  1:472221/149‑1 (MQ=255)                                                                            agaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTg                                                                    <  3:312355/149‑1 (MQ=255)                                                                               aataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTg                                                                 <  1:548986/149‑1 (MQ=255)                                                                                ataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCGTGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGa                                                                <  4:226094/149‑1 (MQ=255)                                                                                ataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTg                                                                 >  2:388829/1‑148 (MQ=255)                                                                                ataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTg                                                                 <  1:388829/148‑1 (MQ=255)                                                                                ataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGa                                                                <  4:566437/149‑1 (MQ=255)                                                                                  aaaCTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAgg                                                              <  3:97961/149‑1 (MQ=255)                                                                                   aaCTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGt                                                                         <  1:165864/137‑1 (MQ=255)                                                                                   aaCTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGt                                                                         >  2:165864/1‑137 (MQ=255)                                                                                   aaCTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGt                                                             >  2:37364/1‑149 (MQ=255)                                                                                    aCTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATAttt                                                                                          >  4:572097/1‑119 (MQ=255)                                                                                    aCTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATAttt                                                                                          <  3:572097/119‑1 (MQ=255)                                                                                    aCTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGtt                                                            <  3:372867/149‑1 (MQ=255)                                                                                      ttCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTAt                                                          >  1:99217/1‑149 (MQ=255)                                                                                      ttCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTAt                                                          >  2:709973/1‑149 (MQ=255)                                                                                      ttCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTAt                                                          >  4:40159/1‑149 (MQ=255)                                                                                          aTCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAgg                                                              <  1:722218/141‑1 (MQ=255)                                                                                          aTCCGGCGCTCGTGGGTGTGACACGGCATCGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAgg                                                              >  2:722218/1‑141 (MQ=255)                                                                                            ccGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTa                                                    >  4:204957/1‑149 (MQ=255)                                                                                              ggCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAAc                                                  <  2:396422/149‑1 (MQ=255)                                                                                               gcgcTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTc                                                                                                                                                <  1:410979/54‑1 (MQ=255)                                                                                               gcgcTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTc                                                                                                                                                >  2:410979/1‑54 (MQ=255)                                                                                                    cGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCCGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACt                                            <  4:487234/149‑1 (MQ=255)                                                                                                      tGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTgg                                          >  3:436854/1‑149 (MQ=255)                                                                                                       gggTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGga                                         <  4:586876/149‑1 (MQ=255)                                                                                                          tgtgACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGAt                                      >  2:522549/1‑149 (MQ=255)                                                                                                            tgACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATAc                                    <  3:329054/149‑1 (MQ=255)                                                                                                                     aTGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGCGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCgg                           >  4:504820/1‑149 (MQ=255)                                                                                                                      tGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGgt                          <  1:128908/149‑1 (MQ=255)                                                                                                                          gtTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGa                      <  1:734822/149‑1 (MQ=255)                                                                                                                              cacaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACttt                  >  4:48454/1‑149 (MQ=255)                                                                                                                               acaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACtttt                 >  2:588071/1‑149 (MQ=255)                                                                                                                               acaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACtttt                 >  4:267712/1‑149 (MQ=255)                                                                                                                                   gTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTtgg             <  3:40159/149‑1 (MQ=255)                                                                                                                                    ttACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTtggt            <  1:262191/149‑1 (MQ=255)                                                                                                                                       gacGGAGCTGGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTggg         >  4:553787/3‑149 (MQ=255)                                                                                                                                       gacGGAGCTGGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTggg         >  2:187228/3‑149 (MQ=255)                                                                                                                                        ccGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGa        <  1:210991/149‑1 (MQ=255)                                                                                                                                              ccAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACt                                            >  2:335068/1‑107 (MQ=255)                                                                                                                                              ccAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACt                                            <  1:335068/107‑1 (MQ=255)                                                                                                                                              ccAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGATAcgtc  >  4:702892/1‑149 (MQ=255)                                                                                                                                              ccAGCTCCGTAAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGAAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGATAcgtc  >  4:493417/1‑149 (MQ=255)                                                                                                                                                aGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTaa                                                   <  4:49975/98‑1 (MQ=255)                                                                                                                                                aGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTaa                                                   >  3:49975/1‑98 (MQ=255)                                                                                                                                                |                                                                                                                                                   GCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGATACGTC  >  NC_000913/2983431‑2983720

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