Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,362,553	G→T	54.1%	noncoding (74/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,553	0	G	T	54.1%	4.5 / 65.3	68	noncoding (74/76 nt)	pheU	tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base G (16/15);  new base T (25/12);  total (41/27)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

CGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGA  >  NC_000913/4362408‑4362700                                                                                                                                                  |                                                                                                                                                    cGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGc                                                                                                                                                  >  1:762829/1‑149 (MQ=255) cGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGc                                                                                                                                                  >  2:629065/1‑149 (MQ=255)    aTTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCg                                                                                                                                               <  2:893486/149‑1 (MQ=255)    aTTGGGATATGCCCGCAAAATGCTTGTGATGTGGGGAGAATCTGGATGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTAAGTGGAGGTCAGTTTTTATGAATTTGGTggcag                                                                                                                                               >  1:522440/1‑145 (MQ=255)          atatGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAACCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCG                                                                                                                                         <  1:550774‑M1/149‑1 (MQ=255)          atatGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCG                                                                                                                                         <  2:413250‑M1/149‑1 (MQ=255)          atatGCCCGCAAATTGCTGCTGATGTGGGGAGAATGTGGTAGAGTTCGGGATAATTGATTTGGAGGCAGAACGCTAAAATCGTGCCGTCCTGAAACGAAAAACGCACCTCCCTGGACGGCCGTTTATATGAATTTGGTGCCCGGACTCg                                                                                                                                         <  2:325838/149‑1 (MQ=255)              gCCCGCAAATTGCTGGTGATGTGGGCAGAATCTGGTGGAGGTCGGTAGAATTGATTTGGAGGCAGCACGCTTAAATCGTGGCGTCCTGAATCGAAAAACGGACCTCGGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAAt                                                                                                                                     <  1:322869/149‑1 (MQ=255)                  gCAAATTGCTGGGGATGTGGGGAGAAGCTGGTTGCGTCCGGTAGAATTGATAGGGAGGCAGAACGCGTAACTCGTGGCGGCCTGAAACGCAAACCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGaa                                                                                                                                 <  2:203006/149‑1 (MQ=255)                             gTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACg                                                                                                                      >  2:973841/1‑149 (MQ=255)                               caggagcggACAATCAGGTTGAGAGCGCTCGAATTGATTTGGAGGCACAACGCTCAAATCGAGGCGTCCTGACACGAAAAACGGACCTCCGGGGAGGTCCGTTTATATGAATTTCGTGCGCGGACTCGCAAACGAACCAACGACACggg                                                                                                                    <  2:98731/142‑1 (MQ=255)                                     gggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGCACCAAGGACACGGGGAtttt                                                                                                              >  2:414430/1‑149 (MQ=255)                                       gagaATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCa                                                                                                            >  1:217189/1‑149 (MQ=255)                                       gagaATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGAGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCa                                                                                                            >  1:71554/1‑149 (MQ=255)                                            tCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATccc                                                                                                       >  1:329806/1‑149 (MQ=255)                                            tCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGTAAACCGGGATTTTCAATccc                                                                                                       >  2:617392/1‑149 (MQ=255)                                              gtgTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCt                                                                                                     <  1:997840/147‑1 (MQ=255)                                               ggTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTg                                                                                                    >  2:183367/1‑149 (MQ=255)                                               ggTTGAGTTCGGTAGAATAGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTCCCCGGACACGGAATCGAACCAAGGACACGGGGATATTCAATAACCTg                                                                                                    >  1:740458/1‑149 (MQ=39)                                                gcggAGTTCGGTAGTATTGATTTGGAGGAAGAACGCTTAAATCGTGGCGTCCTGATACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGc                                                                                                   <  1:227755/146‑1 (MQ=255)                                                gTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGc                                                                                                   <  1:223996/149‑1 (MQ=255)                                                   gAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGCACCAAGGACACGGGGATTTTCAATCCCCTGctct                                                                                                >  1:643683/1‑149 (MQ=255)                                                          gTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGTAATCGAACCAAGGACACGGGGATTTTCAATCCACTGCTCTACCGACt                                                                                         >  1:194216/1‑149 (MQ=255)                                                            aGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa                                                                                       <  2:942554/149‑1 (MQ=255)                                                            aGAATTGATTAGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCTGGGAGGTCCGTATATATGAATTTGTTGCCCGGACTCGGAATCGAACCACGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa                                                                                       <  2:890347/149‑1 (MQ=37)                                                             gAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTATATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACTGGGATTTTCAATCCACAGCTCTACCGACTGAg                                                                                      >  1:479925/1‑149 (MQ=38)                                                              aaTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGc                                                                                     >  2:940512/1‑149 (MQ=255)                                                              aaTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGCCACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGc                                                                                     >  1:992845/1‑149 (MQ=255)                                                              aaTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGc                                                                                     >  1:725039/1‑149 (MQ=255)                                                                 tGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTAt                                                                                  <  2:245328/149‑1 (MQ=255)                                                                   aTTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAACCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATcc                                                                                <  2:398460/149‑1 (MQ=255)                                                                      tGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCggg                                                                             <  1:528374/149‑1 (MQ=255)                                                                           gCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGAGTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACg                                                                        >  1:660815/1‑149 (MQ=255)                                                                             aGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGCCTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGTGCGATCCGGGCCACggg                                                                      >  2:855908/1‑149 (MQ=25)                                                                              gAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACgggg                                                                     >  2:217090/1‑149 (MQ=255)                                                                                 agCTTAAATCGTGGCGTCCTGACACGAAAAACGGACCGCCGGGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAAGCAACGACACGGGGATTTTCAATCCCCTGATCTACCGACTGAGCTATCCGGGCAACGGGgcgc                                                                  <  1:469203/148‑1 (MQ=21)                                                                                   ctgaaatcgTGGCGTCCTGAAACGAAACACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGaaccaaggactcggggattttcaatcccctgctctaccgtctgcgctatccgggcaacggggcgcattaaacc                                                                <  1:52721‑M1/146‑74 (MQ=255)                                                                                    ttAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGCGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAtt                                                               <  1:629065/149‑1 (MQ=255)                                                                                       aaTCGTGGCGTCCTGACACGAAAAACGGACCTCCGTGGAGGTCCGTTTATTTGAATTTGTTGCCCGGACTCGGAATCGATCCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCGATCCGGGCAACGGGGGGCAGTaaa                                                            >  2:577865/1‑149 (MQ=21)                                                                                       aaTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTaaa                                                            <  1:183367/149‑1 (MQ=255)                                                                                                   cTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACAAATGACACGGGGATTTTCAATCACCTGCTATACCGAAGGAGCTAACCGGGCACCGGGGCGCATTTAACCGAAATCTGca                                                >  1:268621/1‑149 (MQ=21)                                                                                                   cTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTATATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGca                                                <  2:2266/149‑1 (MQ=255)                                                                                                    tGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGCCTCGGAATCGAACCAAGGACACGGGGATATTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGcac                                               <  1:343977/149‑1 (MQ=255)                                                                                                     gAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGcaca                                              >  1:377320/1‑149 (MQ=255)                                                                                                        aCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACAtct                                           >  2:924360/1‑149 (MQ=255)                                                                                                          tataaaCGGACCTCCACGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAAGCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGCCTGAGCTAGCCGGGCAACGGGGCGCATTAATCCGTAATCTGCACATCTCg                                         >  1:312000/4‑149 (MQ=21)                                                                                                          gAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCGACCGACTGAGGTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCg                                         >  2:859102/1‑149 (MQ=255)                                                                                                          gAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGaaccaaggacacggggattttcaatcccctgctctaccgactgagctatccgggcaacggggcgcattaaaccgtaatctgcacatctcgtcaacc                                         >  2:378999‑M1/1‑53 (MQ=255)                                                                                                          gAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCg                                         >  2:910089/1‑149 (MQ=255)                                                                                                           ataaaCGGACCTCCACGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGt                                        >  2:282845/3‑149 (MQ=255)                                                                                                           ataaaCGGACCTCCACGGAGGTCCGGTTATTTGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGCCGGAGCTATCCGGGCAGCGGGGCGCATTTCACAGTAAGCTGCACATCTCGt                                        >  1:1003526/3‑149 (MQ=17)                                                                                                             aaaCGGACCTCCGTTGAGGTCCTTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTAGCGACTGAGATATCCGGGCTACGGGGCGCATTAATCCGTATTCTGCACATGTCGTCa                                      >  2:765260/1‑149 (MQ=18)                                                                                                                  gACCTCCACGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTa                                 >  2:786711/1‑149 (MQ=255)                                                                                                                    ccTCCGGGGATGTCCGTTTATATGAATGTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCGGCTCTACCGACTTAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAt                               <  2:630891/149‑1 (MQ=18)                                                                                                                     cTCCGTGGAGGTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGCGATTTTCAAACCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAAACTGCACCTCTCGTCAACCTTAtt                              >  1:850763/1‑149 (MQ=21)                                                                                                                     cTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAtt                              <  2:424988/149‑1 (MQ=255)                                                                                                                     cTCCGTGGAGGTCCGTGTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAtt                              <  1:199316/149‑1 (MQ=255)                                                                                                                               gTCCGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAg                    >  2:739960/1‑149 (MQ=255)                                                                                                                               gTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGGAATCTGCACATATCGTCAACCTAATTTCAGTAAAAg                    >  2:914709/1‑149 (MQ=255)                                                                                                                               gTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACGGCGCTATCCGGGCAACGGGGCGCATTAAACCGTTATCTGCACATCTGGTCATCGTAATTTCAGTAAAAg                    >  2:829605/1‑149 (MQ=18)                                                                                                                                  cGTTTATATGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGTGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGa                 <  1:973841/149‑1 (MQ=255)                                                                                                                                       agatGAAATTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCAATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTAGTCAACCTAATTTCAGTAAAAGCGATTCaa            >  1:328295/3‑149 (MQ=255)                                                                                                                                         atGAATTTGGTGCCCGGACTCGaaccaaggactcggggattttcaatcccctgctctaccgactgagctatccgggcaacggggcgcattaaaccgtaatctgcacatctcgtcaacctaatttcagtaaaagcgattcaactgcttaa          <  1:378999‑M1/149‑128 (MQ=255)                                                                                                                                          tGAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGCCTGAGCTAGCCGGGCTACGGGGAGCATTAATCCGTGATCTGCACATGGCGTCACCCTAATTTCAGTAAAAGCGATTCAACTg         >  1:309798/1‑149 (MQ=18)                                                                                                                                           gAATTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAGTAAACCGTAATCTGCGCATCTCGTCAACCGAATTTCAGTAAAAGCGATTCAACTgg        >  1:389702/1‑148 (MQ=255)                                                                                                                                            aaTTTGTTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCt       >  1:357658/1‑149 (MQ=255)                                                                                                                                            aaTTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCt       >  2:122361/1‑149 (MQ=255)                                                                                                                                               ttGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTaa    <  2:71554/149‑1 (MQ=255)                                                                                                                                                 ggTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGATCTACCGCCTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTATTTTCAGTAAAAGCGATTCAAGTggtgtagg  >  2:187293/1‑142 (MQ=18)                                                                                                                                                  |                                                                                                                                                    CGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGA  >  NC_000913/4362408‑4362700

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