Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,910,295	A→T	6.1%	D7V (GAC→GTC)	yebQ →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,910,295	0	A	T	6.1%	336.1 / 10.4	98	D7V (GAC→GTC)	yebQ	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (48/44);  new base T (2/4);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

TATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTA  >  NC_000913/1910148‑1910443                                                                                                                                                    |                                                                                                                                                     tataATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGAc                                                                                                                                                     >  3:172333/1‑149 (MQ=255)    aaTTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGc                                                                                                                                                  <  2:351157/149‑1 (MQ=255)      tttCGCACCAGGTTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGACTCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCt                                                                                                                                                <  2:1059830/149‑1 (MQ=255)      tttCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCt                                                                                                                                                >  2:979611/1‑149 (MQ=255)      tttCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCt                                                                                                                                                >  4:30250/1‑149 (MQ=255)      tttCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCt                                                                                                                                                >  2:108726/1‑149 (MQ=255)        tCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGc                                                                                                                                              <  4:1079069/149‑1 (MQ=255)         cGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGcc                                                                                                                                             >  1:822978/1‑149 (MQ=255)         cGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGcc                                                                                                                                             <  3:968191/149‑1 (MQ=255)          gCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaac                                                                                                                                            <  4:866802/149‑4 (MQ=255)          gCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaac                                                                                                                                            <  2:12096/149‑4 (MQ=255)              cAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaactttt                                                                                                                                        <  3:1001739/149‑8 (MQ=255)              cAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGc                                                                                                                                        >  1:858009/1‑149 (MQ=255)              cAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTGCAGGCCGTCGGCCTgaactttt                                                                                                                                        <  4:78154/149‑8 (MQ=255)                gggTGCTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGccc                                                                                                                                      <  3:993557/149‑1 (MQ=255)                 ggtggtCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGcccc                                                                                                                                     <  3:30250/149‑1 (MQ=255)                        gCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATa                                                                                                                              <  1:104570/149‑1 (MQ=255)                        gCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATa                                                                                                                              <  1:108726/149‑1 (MQ=255)                         cAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATAc                                                                                                                             >  3:689350/1‑149 (MQ=255)                         cAATCCATCTTTTGCCGGTTAGTTACAAGTCTGCGACATCCACCGTGAAGATCAGTGCTAGAATCATACCCCTGTTGATTATGCACCAAAGAGATAAAATTCCTATCCCAAAAGATCAGCCCGACGGCCTGCCATTGCCGCGGCGGTAc                                                                                                                             >  3:193641/1‑149 (MQ=255)                            tCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGt                                                                                                                          <  3:339458/149‑1 (MQ=255)                               aTCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCg                                                                                                                       >  4:439630/1‑149 (MQ=255)                                  ttttGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGa                                                                                                                      <  3:107278/147‑1 (MQ=255)                                  ttttGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGa                                                                                                                      >  4:107278/1‑147 (MQ=255)                                      gCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTaa                                                                                                                <  1:880450/149‑1 (MQ=255)                                       ccGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAAc                                                                                                               <  3:151719/149‑1 (MQ=255)                                        cGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAAcc                                                                                                              >  4:226601/1‑149 (MQ=255)                                         ggTTAGTTACATTTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCa                                                                                                             >  1:818344/1‑149 (MQ=255)                                           ttagttaCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCAtt                                                                                                           <  1:303514/149‑1 (MQ=255)                                           ttagttaCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCAtt                                                                                                           <  3:1058060/149‑1 (MQ=255)                                              gttaCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCt                                                                                                                                                >  3:183083/1‑109 (MQ=255)                                              gttaCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCt                                                                                                                                                <  4:183083/109‑1 (MQ=255)                                               ttaCAATTCTGCGACAGCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGa                                                                                                       <  3:920945/149‑1 (MQ=255)                                                taCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGAt                                                                                                      <  3:688522/149‑1 (MQ=255)                                                taCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGAt                                                                                                      <  3:773286/149‑1 (MQ=255)                                                   aaTTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATtgt                                                                                                         >  3:95160/1‑143 (MQ=255)                                                   aaTTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATtgt                                                                                                         <  4:95160/143‑1 (MQ=255)                                                   aaTTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTgg                                                                                                   >  3:991518/1‑149 (MQ=255)                                                   aaTTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTgg                                                                                                   >  3:31631/1‑149 (MQ=255)                                                   aaTTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCAGACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTgg                                                                                                   <  4:156027/149‑1 (MQ=255)                                                     ttCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTa                                                                                                 <  2:208596/149‑1 (MQ=255)                                                      tCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTAt                                                                                                <  3:766234/149‑1 (MQ=255)                                                              aTCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATgg                                                                                        <  4:790973/149‑1 (MQ=255)                                                                ccACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGcc                                                                                      <  1:1075632/149‑1 (MQ=255)                                                                   ccGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTc                                                                                   <  3:867548/149‑1 (MQ=255)                                                                     gTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCt                                                                                 >  2:681816/1‑149 (MQ=255)                                                                       gAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTg                                                                               <  1:174629/149‑1 (MQ=255)                                                                            tCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGgc                                                                          <  1:974050/149‑1 (MQ=255)                                                                              aGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGgcgc                                                                        >  3:1021449/1‑149 (MQ=255)                                                                                   tAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCg                                                                   >  4:513165/1‑149 (MQ=255)                                                                                   tAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCg                                                                   >  3:228104/1‑149 (MQ=255)                                                                                   tAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCg                                                                   >  1:323307/1‑149 (MQ=255)                                                                                    aGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGc                                                                  >  1:696868/1‑149 (MQ=255)                                                                                      aaTCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCa                                                                >  4:641878/1‑149 (MQ=255)                                                                                             ccccTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGc                                                         <  3:792487/149‑1 (MQ=255)                                                                                              cccTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGAt                                                                                                      >  2:209445/1‑103 (MQ=255)                                                                                              cccTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGAt                                                                                                      <  1:209445/103‑1 (MQ=255)                                                                                              cccTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTg                                                                                                    >  3:218530/1‑105 (MQ=255)                                                                                              cccTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTg                                                                                                    <  4:218530/105‑1 (MQ=255)                                                                                              cccTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGcc                                                        >  2:201214/1‑149 (MQ=255)                                                                                                 tGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTg                                                     >  3:882317/1‑149 (MQ=255)                                                                                                   ttGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGcc                                                   >  2:953938/1‑149 (MQ=255)                                                                                                        tattCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCcaacaa                                              >  3:788709/1‑149 (MQ=255)                                                                                                            cACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGc                                          <  3:985487/149‑1 (MQ=255)                                                                                                             aCCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGcc                                         <  3:513165/149‑1 (MQ=255)                                                                                                             aCCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGcc                                         >  3:350893/1‑149 (MQ=255)                                                                                                              ccAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCGATCGCCAACGTCGCCCTGCCAACAATCGCCa                                        >  1:1070917/1‑149 (MQ=255)                                                                                                              ccAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCa                                        >  4:431416/1‑149 (MQ=255)                                                                                                              ccAAAGAGATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCa                                        >  2:951775/1‑149 (MQ=255)                                                                                                               cAAAGATATAACATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATc                                            >  3:613494/1‑144 (MQ=255)                                                                                                               cAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATc                                            <  4:613494/144‑1 (MQ=255)                                                                                                               cAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCAc                                       <  4:1021449/149‑1 (MQ=255)                                                                                                                aaaGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACg                                      >  4:981938/1‑149 (MQ=255)                                                                                                                aaaGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACg                                      >  1:488040/1‑149 (MQ=255)                                                                                                                aaaGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACg                                      >  3:258330/1‑149 (MQ=255)                                                                                                                aaaGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCACTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACg                                      >  4:571791/1‑149 (MQ=255)                                                                                                                       taAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTc                               >  2:521485/1‑149 (MQ=255)                                                                                                                        aaaaTTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCa                              >  3:173580/1‑149 (MQ=255)                                                                                                                          aactCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCAtc                            <  4:171524/146‑2 (MQ=255)                                                                                                                               cTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGa                                    >  3:805712/1‑136 (MQ=255)                                                                                                                               cTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGa                                    <  4:805712/136‑1 (MQ=255)                                                                                                                               cTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTc                               >  2:314111/1‑141 (MQ=255)                                                                                                                               cTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTc                               <  1:314111/141‑1 (MQ=255)                                                                                                                               cTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACg                       >  1:646736/1‑149 (MQ=255)                                                                                                                                tATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgc                      <  3:259610/149‑1 (MQ=255)                                                                                                                                 aTGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgcc                     >  2:124314/1‑149 (MQ=255)                                                                                                                                    ccAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCAt                             <  4:408137/138‑1 (MQ=255)                                                                                                                                    ccAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCAt                             >  3:408137/1‑138 (MQ=255)                                                                                                                                    ccAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgccagc                  >  4:288927/1‑149 (MQ=255)                                                                                                                                    ccAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgccagc                  >  1:557666/1‑149 (MQ=255)                                                                                                                                     cAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgccagcc                 >  2:912093/1‑149 (MQ=255)                                                                                                                                      aaaaGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgccagcca                >  3:421785/1‑149 (MQ=255)                                                                                                                                         aGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGtt             <  3:566082/149‑1 (MQ=255)                                                                                                                                           ttCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTcc           <  3:981938/149‑1 (MQ=255)                                                                                                                                           ttCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTcc           <  4:350893/149‑1 (MQ=255)                                                                                                                                              aGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATc        <  4:882317/149‑1 (MQ=255)                                                                                                                                                gCCGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTg      >  1:949493/1‑149 (MQ=255)                                                                                                                                                 ccGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTgg     <  2:71391/149‑1 (MQ=255)                                                                                                                                                    aCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTa  <  1:912093/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTA  >  NC_000913/1910148‑1910443

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