Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,362,578	G→T	50.4%	noncoding (49/76 nt)	pheU ←	tRNA‑Phe

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,362,578	0	G	T	50.4%	11.7 / 86.0	84	noncoding (49/76 nt)	pheU	tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base G (25/17);  new base T (18/24);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTG  >  NC_000913/4362430‑4362724                                                                                                                                                     |                                                                                                                                                   ttGCGGGTGATGTGGGGAGAATCTGGTTGAGGTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGGCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTAAATGAATTTGGTGCCCGGACTCGGAATCGAACCAAg                                                                                                                                                    <  2:380059/149‑1 (MQ=255)   gCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGTCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCCAGGa                                                                                                                                                  >  1:619179/1‑149 (MQ=255)     tGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGaca                                                                                                                                                >  2:325800/1‑149 (MQ=255)                ggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTc                                                                                                                                     <  1:1151694/149‑1 (MQ=255)                ggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTc                                                                                                                                     <  1:36931/149‑1 (MQ=255)                    aaTCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAATCGAAAAACGGACCTCCGTGGAGGTCGGTTTATATGAATTTGGTGCCCGGACTCGGAATCGACCCAAGGACACGCGGATTTTCAAta                                                                                                                                 >  1:676686/1‑148 (MQ=255)                    aaTCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATc                                                                                                                                 >  1:118150/1‑149 (MQ=255)                         ggTTGAGTTCGGTAGAATTGATTTGGAGTCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGAGGAGGTCCGTTTATATGAATTTGGTGCCCGTACTCTGAATAGAACCAATGACACTGGGATTTTCAATCCCCTg                                                                                                                            >  2:778221/1‑149 (MQ=255)                          gTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGGGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGc                                                                                                                           >  2:803004/1‑149 (MQ=255)                               gTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTAc                                                                                                                      <  2:1014194/149‑1 (MQ=255)                                  cGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACACGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGa                                                                                                                   >  2:546526/1‑149 (MQ=255)                                   ggTAGACTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGAc                                                                                                                  <  1:1021647/149‑1 (MQ=255)                                   ggTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGGCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGAc                                                                                                                  <  1:1026353/149‑1 (MQ=255)                                     tAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTg                                                                                                                >  1:56946/1‑149 (MQ=255)                                      aGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa                                                                                                               <  1:853283/149‑1 (MQ=255)                                       gAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAg                                                                                                              <  1:138836/149‑1 (MQ=255)                                       cttcagATCTGGAGGCAGAACGCTTAAATCGTGGAGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGTGATTTTCAATCCCCTGCTCTACCGACTGAg                                                                                                              <  1:97054/144‑1 (MQ=255)                                          ttGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATAAGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCGTCCGACTGAGCTa                                                                                                           >  1:1177294/1‑149 (MQ=255)                                              tttGGAGGCATAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCg                                                                                                       >  2:1130446/1‑149 (MQ=255)                                               ttGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCgg                                                                                                      <  2:90908/149‑1 (MQ=255)                                                     gCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACg                                                                                                <  1:457616/149‑1 (MQ=255)                                                      cAGAACGCTTAAATCGTGGCGTCCTGAAACTAAAAACGGACCTCCGTGGAGGTCCGTTTATATGTATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACgg                                                                                               <  2:85801/149‑1 (MQ=255)                                                         aaagCTTAACTCGTGGCGTCCTGAAACGAAAAACGGACCTACGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGgc                                                                                            <  2:191300/146‑1 (MQ=255)                                                         aaCGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGgc                                                                                            >  1:97449/1‑149 (MQ=255)                                                                 aaTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTaaa                                                                                    <  2:633546/149‑1 (MQ=255)                                                                 aaTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCGCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTaaa                                                                                    <  1:660579/149‑1 (MQ=255)                                                                  aTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAGTAAAc                                                                                   <  1:614799/149‑1 (MQ=255)                                                                    cGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCg                                                                                 <  1:235628/149‑1 (MQ=255)                                                                          gTCCTGAAACGAAAAACCGACCGCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGTACCAAGGCCACGGGGATTTTAAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCt                                                                           <  1:677389/149‑1 (MQ=21)                                                                            ccTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGc                                                                         <  1:430443/149‑1 (MQ=255)                                                                             cTGAAACGAAAAACGGACCTCCGTGGAGCTCCGTTTATAGGAATTAGGTGCCCCGACTCGGAATCGACCGACGTACACGTGGCTTTTCAATCCCCTGCTATACCGACTGACCTATCCGGGCATCGCGGCGCATCAAACCGTAATCAGca                                                                        >  2:81109/1‑149 (MQ=255)                                                                              tGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACAGAGCTATCCGGGCAACGGGGCGCATTAAGCCGTAATCTGcac                                                                       >  2:223542/1‑149 (MQ=255)                                                                                aaaCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGCCACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACAt                                                                     >  1:597111/1‑149 (MQ=255)                                                                                aaaCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACCCGGGGATGTTCAATCCCCTGCTCTCCCGACTGCGCTATCCGGGCAACGGGGCTCATTAAACCGTAATCTGCACAt                                                                     >  1:670349/1‑149 (MQ=21)                                                                                 aaCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGCGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAAGCGTAATCTGCACAtg                                                                    >  2:436956/1‑148 (MQ=255)                                                                                  aCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCCCAtct                                                                   >  2:129483/1‑149 (MQ=255)                                                                                  aCGAAAAACGGACCTCCGTGGAGGTCCGGTTATATGAATATGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGAAATCTGCACAtct                                                                   >  1:374450/1‑149 (MQ=255)                                                                                    gAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCg                                                                 <  1:789791/149‑1 (MQ=255)                                                                                    gAAAAACGGACCTCCGAGGAGGTCAGTTTATATGAATTTGGTGACCGGACTCGGAATCGAACCAATGACACGGGGAATTTCAATCCCCTGCTCTACCGACTGAGCTATCCCGCCAACGGTGCGCATTAAACCGTAATCTGAACATCTCg                                                                 >  1:1031752/1‑149 (MQ=21)                                                                                     aaaaaCGGACCTCCGTGGAGGTCCGTTTATATGATTTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCAGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAGCCGGAATCTGCAGATGTCGt                                                                >  2:584945/1‑149 (MQ=18)                                                                                      taaaCGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTc                                                               >  1:404721/2‑149 (MQ=255)                                                                                      aaaaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTc                                                               <  1:803004/149‑1 (MQ=255)                                                                                       aaaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCa                                                              >  2:1005430/1‑149 (MQ=255)                                                                                         aCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCacc                                                            >  1:199611/1‑147 (MQ=255)                                                                                           ggACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCt                                                          >  1:969188/1‑149 (MQ=255)                                                                                              ccTCCGTGGAGGTCCGGTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCACCTGCTCTACCGACTGAGCTATCCGGGCAACGGGTCGCATTAAACCGTAATCTGCACATCTCGGCAACCTAAt                                                       >  2:666937/1‑149 (MQ=21)                                                                                                tCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAttt                                                     <  1:546526/149‑1 (MQ=255)                                                                                                  cGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGCCTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCa                                                   <  2:48334/149‑1 (MQ=255)                                                                                                  cGTGGAGGTCCGTTTATAGGAATTTGGTGCCCTGACTCGGAATCGAACCAATGACACGGTGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCa                                                   <  1:785646/149‑1 (MQ=255)                                                                                                      gAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTaa                                               >  2:95868/1‑149 (MQ=255)                                                                                                      gAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAAACCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCTCATAAAACCGTAATCTGCACATCTCGTCAACCTAATTTCATTaa                                               >  2:1069779/1‑149 (MQ=21)                                                                                                        ggTCGGTTTATATGAAATTGGTGCCCGGTCTGGGAATCTAACCAATGACACGGGGTTGTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTaaaa                                             <  2:1137817/149‑1 (MQ=21)                                                                                                           ccGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCg                                          >  2:16424/1‑149 (MQ=255)                                                                                                            cGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGa                                         >  2:488928/1‑149 (MQ=255)                                                                                                            cGTTTATATGAATTTGGTGCACGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGa                                         >  1:580456/1‑149 (MQ=255)                                                                                                             gTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGAt                                        >  1:325452/1‑149 (MQ=255)                                                                                                              tttATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGAtt                                       <  1:481209/149‑1 (MQ=255)                                                                                                              tttATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAATAAACCGTAAACAGCACATCACCTCAACCTAATTTAAGAAAAAGCGagt                                       >  1:766876/1‑147 (MQ=18)                                                                                                               ttATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCACCCTAATTTCAGTAAAAGCGATTc                                      >  2:404981/1‑149 (MQ=255)                                                                                                                  tatGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAAc                                   >  2:973821/1‑149 (MQ=255)                                                                                                                      aaTTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCGATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCt                               <  2:943806/149‑1 (MQ=255)                                                                                                                       aTTTGGTGCCCTGAGTCTGAATCTAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGAATGAGCTATCCGGGCAACGGGGCGCTTTAAACCGTAATCTGCTCATCTCTTCAAGCTAATTTCAGTAAAAGCGATTCAACTGCtt                              <  2:44051/149‑1 (MQ=18)                                                                                                                        tttGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTa                             <  1:325800/149‑1 (MQ=255)                                                                                                                         ttGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTaa                            <  1:45599/149‑1 (MQ=255)                                                                                                                           ggTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGa                          <  2:179471/149‑1 (MQ=255)                                                                                                                             tGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGAtt                        >  1:570472/1‑149 (MQ=255)                                                                                                                              gCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTg                       >  1:442646/1‑149 (MQ=255)                                                                                                                               cccGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGc                      <  2:969188/149‑1 (MQ=255)                                                                                                                                ccGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGTCTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCg                     <  1:95868/149‑1 (MQ=255)                                                                                                                                ccGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCg                     <  1:1069779/149‑1 (MQ=255)                                                                                                                                ccGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCGGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCAGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCg                     <  2:274988/149‑1 (MQ=255)                                                                                                                                      tCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAAt               >  1:793493/1‑149 (MQ=255)                                                                                                                                        ggAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGCGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCACCCTAATTTCAGTAAACGCGATTCAACTGCTTAACATTGCGGCAAATCg             >  2:994770/1‑149 (MQ=255)                                                                                                                                        ggAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAGCCGTAATCTGCACATCTCGTCGACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCTCATCg             >  1:839781/1‑149 (MQ=255)                                                                                                                                         gAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCACCCTAATTTCTGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGc            <  1:404981/149‑1 (MQ=255)                                                                                                                                          aaTCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGct           <  2:544668/149‑1 (MQ=255)                                                                                                                                               aaCCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCttt      >  2:592208/1‑149 (MQ=255)                                                                                                                                                aCCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCtttt     <  2:619179/149‑1 (MQ=255)                                                                                                                                                 ccAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTGTTTTc    >  2:591775/1‑149 (MQ=255)                                                                                                                                                  cAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAACGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCt   <  2:1068220/149‑1 (MQ=255)                                                                                                                                                  cAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCt   <  2:379770/149‑1 (MQ=255)                                                                                                                                                  cAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAATTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCt   <  2:1072738/149‑1 (MQ=255)                                                                                                                                                   aaGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCACCCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCtt  >  2:275801/1‑148 (MQ=255)                                                                                                                                                   aaGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCtg  <  1:485463/149‑1 (MQ=255)                                                                                                                                                   aaGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCtg  >  1:76976/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTG  >  NC_000913/4362430‑4362724

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