Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,423,622	C→T	8.1%	noncoding (34/76 nt)	thrV ←	tRNA‑Thr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,622	0	C	T	8.1%	248.5 / 9.4	87	noncoding (34/76 nt)	thrV	tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base C (50/29);  new base T (3/4);  total (53/34)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.99e-01

ACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCA  >  NC_000913/3423478‑3423766                                                                                                                                                 |                                                                                                                                                 aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAggg                                                                                                                                             <  2:261859/150‑1 (MQ=255)   ggCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTg                                                                                                                                           <  1:1327537/150‑1 (MQ=255)    gCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTga                                                                                                                                          <  2:1644402/150‑2 (MQ=255)    gCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgc                                                                                                                                          >  1:648252/1‑150 (MQ=255)     cGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgag                                                                                                                                         <  2:1079675/150‑3 (MQ=255)      gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgagg                                                                                                                                        <  1:944512/150‑4 (MQ=255)      gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTgcgc                                                                                                                                        >  2:1233500/1‑150 (MQ=255)       tttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGct                                                                                                                                       <  1:572497/150‑1 (MQ=255)        ttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgaggtc                                                                                                                                      <  2:1788375/150‑6 (MQ=255)         tCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGctct                                                                                                                                     >  1:1030031/1‑150 (MQ=255)         tCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGctct                                                                                                                                     >  1:1134868/1‑150 (MQ=255)         tCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGctct                                                                                                                                     <  1:614598/150‑1 (MQ=255)            cTTCTGAGTTCGGCATGGGGGCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAcc                                                                                                                                  <  2:960019/150‑1 (MQ=255)                   gTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAg                                                                                                                           >  2:1653491/1‑150 (MQ=255)                   gTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAg                                                                                                                           <  2:1583407/150‑1 (MQ=255)                   gTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAg                                                                                                                           >  1:1067964/1‑150 (MQ=255)                     tCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGcc                                                                                                                         <  1:1556558/150‑1 (MQ=255)                      cGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCa                                                                                                                        <  1:1733660/150‑1 (MQ=255)                           tGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCa                                                                                                                   <  1:1654306/150‑1 (MQ=255)                               gTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCAc                                                                                                               <  2:1067964/150‑1 (MQ=255)                                   ggTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTGTAATCTGCCGAACTTTAGGGGTAAAAGTGGTGCTGATACCCAGAGTCGAGCTGGGGGCGTCACGCTTACCAGGGGTGCGCTGTACCAACTGCGCCATAGCAGCACGGTa                                                                                                           >  2:1579937/1‑150 (MQ=255)                                   ggTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTa                                                                                                           >  1:1279555/1‑150 (MQ=255)                                         accaccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTg                                                                                                     <  2:869393/150‑1 (MQ=255)                                           caccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGAt                                                                                                   <  1:1899921/150‑1 (MQ=255)                                                cgcTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTg                                                                                              <  2:738416/150‑1 (MQ=255)                                                cgcTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTg                                                                                              <  2:626315/150‑1 (MQ=255)                                                       gccgccATACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTc                                                                                       <  2:852228/150‑1 (MQ=255)                                                       gccgccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTc                                                                                       >  1:227952/1‑150 (MQ=255)                                                       gccgccAGACAAATTCTTTTCTAATCTGCCGAACATTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGATCTGGGGACCTCACCCTTACCAAGGGTGCGCTGTACCAACTGAGGCATATCAGCAGGCTAATGTTGATGCCTGGCAGttg                                                                                       >  2:1173579/1‑149 (MQ=37)                                                        ccgccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTcc                                                                                      >  1:1153032/1‑150 (MQ=255)                                                        ccgccAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTcc                                                                                      >  1:935146/1‑150 (MQ=255)                                                              gACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTg                                                                                              <  1:1903944/136‑1 (MQ=255)                                                              gACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTg                                                                                              >  2:1903944/1‑136 (MQ=255)                                                               aCAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTActct                                                                               >  2:681585/1‑150 (MQ=255)                                                               aCAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTActct                                                                               >  1:1985650/1‑150 (MQ=255)                                                                   aTTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCa                                                                           >  1:990115/1‑150 (MQ=255)                                                                      cTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCGCCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTATCGCATgg                                                                        >  2:1637380/1‑150 (MQ=255)                                                                      cTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCGCCCTTACCAAGGGTGCGCTCGACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATgg                                                                        >  2:1644415/1‑150 (MQ=255)                                                                        tttCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATgggg                                                                      <  2:1134868/150‑1 (MQ=255)                                                                         ttCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGga                                                                     >  2:791244/1‑150 (MQ=255)                                                                          tCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGgag                                                                    <  2:1030031/150‑1 (MQ=255)                                                                           cTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGgagc                                                                   >  2:1271808/1‑149 (MQ=255)                                                                           cTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCAAGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGgaga                                                                   >  1:1808151/1‑150 (MQ=255)                                                                            tAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTATACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAc                                                                  <  2:1758982/150‑1 (MQ=255)                                                                               tCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAcccc                                                               <  1:1233500/150‑1 (MQ=255)                                                                                cTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCCCGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCca                                                              >  2:1278203/1‑150 (MQ=255)                                                                                 tGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcac                                                             <  1:802773/150‑1 (MQ=255)                                                                                  gCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcaca                                                            >  2:1637533/1‑150 (MQ=255)                                                                                  gCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcaca                                                            >  2:1279235/1‑150 (MQ=255)                                                                                     gAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTa                                                         >  1:1445520/1‑150 (MQ=255)                                                                                      aaCTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTAc                                                        <  1:1653491/150‑1 (MQ=255)                                                                                                 aaaaaGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTa                                             >  2:1412395/1‑150 (MQ=255)                                                                                                   aaaGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACg                                           >  1:279427/1‑150 (MQ=255)                                                                                                         gtgCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                     >  1:127038/1‑150 (MQ=255)                                                                                                              gATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGATCCATATCAGCACGCTAAATTTGATTCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCACACACTACCATCGGCGCTACGGCGTTTCACtt                                >  2:453170/1‑150 (MQ=255)                                                                                                              gATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACtt                                >  2:1700457/1‑150 (MQ=255)                                                                                                              gATACCCAGAGTCGAACTGGGGACCACACCCTTACCAAGGGTGCGCTCTACCAACAGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTATCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACtt                                >  2:601612/1‑150 (MQ=255)                                                                                                                tACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTACCTACTCTCGCATGGGGAGACCCCACCCTACCATCGGCGCTACGGCGTTTCACGTCt                              >  1:1687903/1‑150 (MQ=255)                                                                                                                 aCCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                             <  2:279427/150‑1 (MQ=255)                                                                                                                 aCCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                             >  1:151672/1‑150 (MQ=255)                                                                                                                  cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACTCTAAATTTGATTCCTGGCAGTTCCCTACTCTCGCATGGGGATACCCCACACTACCATCTTCGCTACGGCGTTTCACTTCTGa                            <  2:648252/150‑1 (MQ=25)                                                                                                                  cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                            >  2:988675/1‑150 (MQ=255)                                                                                                                  cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                            >  1:1040379/1‑150 (MQ=255)                                                                                                                  cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                            >  1:1801131/1‑150 (MQ=255)                                                                                                                  cccAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACGGAGCCATAGCAGCACGCTAAATTTGATGCCGGGCAGTTCCCTACTCTCGCATGGGGAGCCCCCACACGACCATCGGCGCGACGGCGTTTCACGTCTGa                            >  1:529420/1‑150 (MQ=25)                                                                                                                          cGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCa                    >  2:1612017/1‑150 (MQ=255)                                                                                                                          cGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACACCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCa                    >  2:1453373/1‑150 (MQ=255)                                                                                                                           gAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCAt                   >  1:926182/1‑150 (MQ=255)                                                                                                                                ggggACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGCCCCCACACTACCATCGGCGCTACGGCGTTTCCCTTCTGAGTTCGGCATggggg              >  1:472585/1‑149 (MQ=255)                                                                                                                                ggggACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGt              >  2:1290011/1‑150 (MQ=255)                                                                                                                                cgagACCTAACACTTACCAAAGATGCGCTATACCAACAGAGCCATATCAGCACCCTAAATTTGATAAATGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGt              >  2:191678/4‑150 (MQ=17)                                                                                                                                 gggACCTCACCCTTACAAACCGTTCGGTCTACCAACTGAGCCTTTTCATGTGGTTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTc             <  1:1551282/150‑1 (MQ=14)                                                                                                                                   gagcgcACCCTTGGTAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg           >  1:762650/6‑150 (MQ=255)                                                                                                                                      cTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg        >  2:219619/1‑150 (MQ=255)                                                                                                                                       gcACCCTTGGTAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTgg       >  1:1057512/2‑150 (MQ=255)                                                                                                                                       gcACCCTTGGTAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTgg       >  2:1187722/2‑150 (MQ=255)                                                                                                                                        cACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg      <  2:127038/150‑1 (MQ=255)                                                                                                                                        cACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCACACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg      <  2:990115/150‑1 (MQ=255)                                                                                                                                        cACCATTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg      >  1:1121301/1‑150 (MQ=255)                                                                                                                                         aCCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGa     <  2:1279555/150‑1 (MQ=255)                                                                                                                                          cccTTACCACGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAATTTTGATGCCTGGCAGTTCCCTACTCTCGCCTGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCCTGGGGTCAGGCgggcc    <  2:529420/150‑3 (MQ=18)                                                                                                                                          cccTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac    >  2:397504/1‑150 (MQ=255)                                                                                                                                            cTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGacca  >  1:626936/1‑150 (MQ=255)                                                                                                                                            cTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGacca  >  2:51467/1‑150 (MQ=255)                                                                                                                                            cTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCACACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGacca  >  2:1549826/1‑150 (MQ=255)                                                                                                                                              tACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCAc                                  >  1:260419/1‑116 (MQ=255)                                                                                                                                              tACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCAc                                  <  2:260419/116‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                 ACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCA  >  NC_000913/3423478‑3423766

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