Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	228,997	T→C	6.1%	noncoding (70/77 nt)	aspU →	tRNA‑Asp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	228,997	0	T	C	6.1%	180.8 / 3.3	65	noncoding (70/77 nt)	aspU	tRNA‑Asp
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base C (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

GGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTCTGAAGCCCATAATCACCTCAGTTA  >  NC_000913/228860‑229138                                                                                                                                          |                                                                                                                                              gggAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGCTCCGCCACttat                                                                                                                                   >  1:295029/1‑150 (MQ=17)  ggAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACttatt                                                                                                                                  >  1:717534/1‑150 (MQ=33)     aCTGCCAGTCATCAACTAACCCCCGAAGCCGCGCCTAAACCCGGCGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGGCCGTTCCGCCACTTATTaag                                                                                                                               <  2:80755/150‑1 (MQ=11)      cTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCACGTCAGTTCCGCCACTTATTaaga                                                                                                                              >  2:302458/1‑150 (MQ=255)        gCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTaagaag                                                                                                                            >  1:1285843/1‑150 (MQ=38)          cAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGACGca                                                                                                                          >  1:1180715/1‑149 (MQ=37)                tCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGt                                                                                                                    <  2:600153/150‑1 (MQ=255)                tCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGGCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGt                                                                                                                    <  2:230400/150‑1 (MQ=255)                 cAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACATGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGtt                                                                                                                   >  2:437944/1‑150 (MQ=255)                     ttAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACg                                                                                                               <  2:170658/150‑1 (MQ=255)                       aaGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCt                                                                                                             <  1:668685/150‑1 (MQ=255)                         gCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCg                                                                                                           >  2:1132643/1‑150 (MQ=255)                           aGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTGGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAg                                                                                                         >  2:753943/1‑150 (MQ=255)                           aGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAg                                                                                                         >  2:596532/1‑150 (MQ=255)                            gTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGt                                                                                                                    >  2:887072/1‑138 (MQ=255)                            gTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGt                                                                                                                    <  1:887072/138‑1 (MQ=255)                            gTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAgg                                                                                                        <  2:336992/150‑1 (MQ=255)                                   ggTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGttttttt                                                                                                 >  1:332405/1‑150 (MQ=255)                                      cATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCg                                                                                              >  2:912114/1‑150 (MQ=255)                                      cATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGCTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCg                                                                                              >  2:976239/1‑150 (MQ=255)                                          aaaCCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTg                                                                                          >  1:1086656/1‑150 (MQ=255)                                             ccGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGtat                                                                                       <  1:1259074/150‑1 (MQ=255)                                             ccGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGATCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGtat                                                                                       >  2:293737/1‑150 (MQ=255)                                                  ggtTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCta                                                                                  >  1:535209/1‑150 (MQ=255)                                                      gTAAAAGTATTCGGTGGAGCGGTAGTTCAGTCGGTTAGATTTCCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTTTTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCtattat                                                                              <  1:400972/150‑1 (MQ=255)                                                      gTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCtattat                                                                              >  1:106686/1‑150 (MQ=255)                                                       tAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTattatt                                                                             >  1:140951/1‑150 (MQ=255)                                                       tAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTattatt                                                                             >  2:1288218/1‑150 (MQ=255)                                                         aaaGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGc                                                                           >  1:742396/1‑150 (MQ=255)                                                          aaGATTTCGGTGGAGCTGTTGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGGCGCGGTTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCAGTATATCTATTATTGcc                                                                          <  2:173527/150‑1 (MQ=25)                                                          aaGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGcc                                                                          <  1:1173287/150‑1 (MQ=255)                                                          aaGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGCTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTTcc                                                                          <  2:295029/150‑1 (MQ=255)                                                           aGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTACGCCACATATTAAGAAGCCTAGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATATattatttaca                                                                         >  1:591869/1‑146 (MQ=21)                                                            gAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAg                                                                        <  1:441211/150‑1 (MQ=255)                                                                   gTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCa                                                                 <  2:932369/150‑1 (MQ=255)                                                                   gTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGCTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCa                                                                 <  1:976239/150‑1 (MQ=255)                                                                        gCGGTAGTTCAGTCTGTTTGAATACCTGCCTGACACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAAt                                                            <  2:979460/150‑1 (MQ=255)                                                                          ggTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATcc                                                          >  2:1205318/1‑150 (MQ=255)                                                                              gTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTg                                                      <  2:803246/150‑1 (MQ=255)                                                                                 cAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCAt                                                   <  2:676533/150‑1 (MQ=255)                                                                                      ggTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGATGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTAc                                              <  2:751030/150‑1 (MQ=255)                                                                                      ggTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTaa                                              >  1:460348/1‑149 (MQ=255)                                                                                      ggTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTAc                                              <  1:1257572/150‑1 (MQ=255)                                                                                         tAGAATACCTGCCTGTCACGCAGGGGGTCACGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTATCGCTCGAGGTTTTTTTTCGTGTGGATATCAATTATTGCAAGAATCCCTAAAAACCTCTGCTTTTTAAGca                                           >  1:233520/1‑149 (MQ=25)                                                                                          aGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGctc                                          <  2:5754/150‑1 (MQ=255)                                                                                            aaTACCTGCCTGTCACGCAGGGGGTCGCGTGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTTTTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCtt                                        <  1:857514/150‑1 (MQ=255)                                                                                             aTACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCttt                                       >  2:619887/1‑150 (MQ=255)                                                                                                cctgcctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTc                                    >  1:866758/1‑150 (MQ=255)                                                                                                  tgcctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCt                                  >  1:1132218/1‑150 (MQ=255)                                                                                                   gcctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTc                                 >  2:1256255/1‑150 (MQ=255)                                                                                                    cctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCa                                <  1:586642/150‑1 (MQ=255)                                                                                                     ctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGGCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTcta                               >  1:1120580/1‑148 (MQ=255)                                                                                                     ctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCaa                               <  1:700196/150‑1 (MQ=255)                                                                                                      tgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAAc                              >  1:1183600/1‑150 (MQ=255)                                                                                                        tCAGGCAGGGGGTGGCGGGTTCGAGTCCCGTCCGTTCCGCCTCTTATTTAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAg                            <  1:590170/150‑1 (MQ=255)                                                                                                         ctcGCAGGGGGTCGCGGGTTCGAGTCCCGTACGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGt                           <  1:856288/148‑1 (MQ=255)                                                                                                         cACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGt                           >  2:1049300/1‑150 (MQ=255)                                                                                                            gCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTCTg                        <  1:625033/150‑1 (MQ=255)                                                                                                             cAGGGGGTCGCTGGTTCGAGTCCCGTCAGTTACGCCAGTTATTAATAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTCTGa                       <  1:42678/150‑1 (MQ=255)                                                                                                                 gggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTCTGAAGcc                   <  2:42223/150‑1 (MQ=255)                                                                                                                  ggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTCTGAAGccc                  <  2:717534/150‑1 (MQ=255)                                                                                                                   gTCGCGGGTTCGAGTACAGTCCGTTCCGCCACTTATTAAGTAGCCTAGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGGATAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTCTGAAGCCCa                 <  1:397567/150‑1 (MQ=255)                                                                                                                       cgGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTCTGAAGCCCATAAt             <  2:332405/150‑1 (MQ=255)                                                                                                                                gTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTCTGAAGCCCATAATCACCTCAGt    >  1:1233194/1‑150 (MQ=255)                                                                                                                                  cccGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAGTCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTCTGAAGCCCATAATCACCTCAGTTa  >  1:644656/1‑150 (MQ=255)                                                                                                                                          |                                                                                                                                              GGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTTATTAAGAAGCCTCGAGTTAACGCTCGAGGTTTTTTTTCGTCTGTATATCTATTATTGCCAGAATCGCAAAAATCCTCTGCATTTTACGCTCTTTTTCCTCAACAGTCTGAAGCCCATAATCACCTCAGTTA  >  NC_000913/228860‑229138

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