Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	780,624	C→T	noncoding (71/76 nt)	lysT →	tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	780,624	0	C	T	96.3%	172.1 / NA	55	noncoding (71/76 nt)	lysT	tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (23/29);  minor base G (0/2);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

AACGCGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAG  >  NC_000913/780378‑780869                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       aaCGCGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAAAGAAAAAGTGG‑GTCGTTAGCTCAGTTGGAAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACtcacc                                                                                                                                                                                                                                                    >  1:529785/1‑251 (MQ=255)            aTGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAg‑g                                                                                                                                                                                                                                        <  2:432579/251‑1 (MQ=255)                    cacacGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCg                                                                                                                                                                                                                                >  2:46194/1‑251 (MQ=255)                            cAGAAGTAGCATTATTTATGGTCATGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGGCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAGCAAACCGAAAGCAACGAAAAAGTGG‑TTCGTTAGATCAGTTGGTAGAGCAGTTGACTTTTACTCAATTGATCGCAGGTTCGAATCCTGCACGACTCACCAATCTATCCA‑GTGGAAGCGgtagtca                                                                                                                                                                                                                         >  2:266592/1‑248 (MQ=255)                               aaGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACg                                                                                                                                                                                                                     >  2:106299/1‑251 (MQ=255)                                 gtggatttagttCTGGTCGTGTCGTGCGAATCATAAGTAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCAGACGTTGACACGGGATATCAAACAAACCGAAAGCCACGAAAAAGTGGTGT‑GTTAGCTCAGTTTGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCTAAACCTGCACGACGCACCAATCGCTAAT‑GTGGAGGCGTTAGTAAAAccga                                                                                                                                                                                                                   <  2:23194/241‑4 (MQ=37)                                 gTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTg                                                                                                                                                                                                                   >  1:96401/1‑251 (MQ=255)                                          tttCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAAc                                                                                                                                                                                                          >  1:282836/1‑251 (MQ=255)                                             cTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGtt                                                                                                                                                                                                       >  2:309827/1‑251 (MQ=255)                                                  cgtgtcgtgCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCAt                                                                                                                                                                                                   <  1:118817/250‑1 (MQ=255)                                                       cgtgCGAATCATAAGCAGTTGAGTGATCTACATCGAAATGTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCa                                                                                                                                                                                              <  1:146647/250‑1 (MQ=255)                                                        gtgCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCGCAAGTCTGAAAGCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTGCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGTGCAAc                                                                                                                                                                                            <  2:526201/251‑1 (MQ=255)                                                            gAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGcc                                                                                                                                                                                        <  1:309827/251‑1 (MQ=255)                                                                     gCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCg                                                                                                                                                                               <  1:497296/251‑1 (MQ=255)                                                                               gagcTACAGCGAAATTTTTGTTGCGCTCAATGCTGAAATCAGTAATAGATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGATCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGATCAACGGCCCGAAGGGCGATGCGAAGTc                                                                                                                                                                     <  2:268415/248‑1 (MQ=37)                                                                                    aCATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAg                                                                                                                                                                  <  1:510533/249‑1 (MQ=255)                                                                                     cATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCa                                                                                                                                                                                                    >  2:31463/1‑214 (MQ=255)                                                                                     cATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCa                                                                                                                                                                                                    <  1:31463/214‑1 (MQ=255)                                                                                         gAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCt                                                                                                                                                           <  1:46194/251‑1 (MQ=255)                                                                                         aaaaTTTTTGTTTCTCTCAAATCTGAAATCATTAAGATATGCCGCCTTTGCCACAGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTTTCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGCAGTAAAACGTGAAGGAGAACGTTGCGTGAGCAACGTCCCTAAGGGCGAGACGAAGTCGATTCATCCt                                                                                                                                                           <  2:88962/250‑1 (MQ=25)                                                                                              ttttttgtcgcTAAAGACTGAAATCAGTAATATAGGCCGCCGTTGCCAAGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGGAGAGCAGTTGACGTTGAAGCAATTGGTCGCAGGTTAGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGTATAACTTTGCATGAGCAACGGCCCGAAGGGCGAGACTAAGTCGAGTCATCCTGCAc                                                                                                                                                       <  2:46379/242‑1 (MQ=255)                                                                                                 tgtttctcTCAAGTCTGAAAACAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGGGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGAc                                                                                                                                                    <  2:548225/243‑1 (MQ=255)                                                                                                         tCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACt                                                                                                                                            >  2:50928/1‑250 (MQ=255)                                                                                                             gTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATGGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACAt                                                                                                                                       <  1:229924/251‑1 (MQ=255)                                                                                                                    aTCAGTAATATAGGCCGCTGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACGTTTAATCAATTGGTCGCAGGTTCGAATCCTGCATGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGTATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTATCAtagttagt                                                                                                                                <  1:356812/251‑1 (MQ=25)                                                                                                                              tatGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCAAGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGAATcc                                                                                                                      >  2:535292/1‑251 (MQ=255)                                                                                                                               atGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCTAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATcc                                                                                                                      <  2:28079/250‑1 (MQ=255)                                                                                                                               atGCCGCCGGTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTGTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCCAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCa                                                                                                                     <  1:286887/251‑1 (MQ=255)                                                                                                                                tGCCGCCGTTGCCACGGGATATCCAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCTAAACA‑GAGGACGCGGAAGTTAAACGCGAAGGATAACGTTGCATAAGCCACGGCCCGAAGGGCGAGCGGAACGCGACTCATCCTGCACGACCCACCACTAAAAtattgggtcgtagtatac                                                                                                                      >  1:108783/1‑233 (MQ=18)                                                                                                                                         tGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTCATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGGAGGGCGAGAGGAAGGAGAGTAATCCGGCGCGAACCACCACAAACAGAGATAGTTggagtagcacgcgtagaatc                                                                                                           >  2:390200/1‑237 (MQ=21)                                                                                                                                          gCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCg                                                                                                          <  1:430876/251‑1 (MQ=255)                                                                                                                                           ccACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCGGCACGACTCACCATTGGCTAAGAGTG‑AAGCGGTCGTAAGACGCGGAGGATAACGTTGAATGAGAAAAGCCCCGCAGGGCGAGCGGAAGACGACTCATCCTGCACGACCCACCACTCAAATAGTAAGTTGTAGGATCCAGAgtaggttagg                                                                                                         >  2:455761/1‑245 (MQ=16)                                                                                                                                            cACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCggg                                                                                                        <  1:155991/251‑1 (MQ=255)                                                                                                                                             aCGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGt                                                                                                       <  1:106299/251‑1 (MQ=255)                                                                                                                                                 gATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAGTTTTTAAG‑GCGGAAGCGGATGTTAAACGTGAAGGATAACGTTGCCTCAGCAACGGCACGTCGGGCGAGCCGAAGTCTGGACATCCGTCACGACCCCCCACCACCATAGTCGGCTGCTGCATGCAGCGTAGTATCGGGAGAtt                                                                                                   >  1:246920/1‑251 (MQ=9)                                                                                                                                                    atCAAACCAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGc                                                                                                >  2:28241/1‑251 (MQ=255)                                                                                                                                                    atCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGTTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCTCCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAg                                                                                                 <  1:345704/250‑1 (MQ=255)                                                                                                                                                        aacaaacCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTc                                                                                              >  1:157498/1‑249 (MQ=255)                                                                                                                                                         acaaacCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGc                                                                                           <  1:286445/251‑1 (MQ=255)                                                                                                                                                           aaacCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAGTCGCTACG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGACACCCTGCACCACCCACCTCTAACATAGGGGGTAGTAGTATCTAGCGTAGAATCggaagcgtagcacagggg                                                                                         >  2:132227/1‑235 (MQ=21)                                                                                                                                                           aaacCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAg                                                                                                                                                                                                <  1:150899/148‑1 (MQ=255)                                                                                                                                                           aaacCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAg                                                                                                                                                                                                >  2:150899/1‑148 (MQ=255)                                                                                                                                                                       aCGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGTAGTTGATTTTTAAGCAATTGGTCGCAGGTTCGAATCCTGCAAGACGCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCTAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGCGAGAGCACCTc                                                                             <  2:23604/251‑1 (MQ=37)                                                                                                                                                                        cGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCCATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCACCGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCCCCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTcc                                                                            >  1:513183/1‑251 (MQ=255)                                                                                                                                                                             aaGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTAc                                                                       <  1:203784/251‑1 (MQ=255)                                                                                                                                                                               gTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAGGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGCCACAACACTAACATGGGTGGTTGAAGTATGCAGCGTAGTATCGGGTGACTAGCTCCGCTTGGAGACCACCGCCCTTaaa                                                                      >  2:291661/1‑248 (MQ=18)                                                                                                                                                                                       gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATATTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAgggg                                                              <  1:225425/251‑1 (MQ=255)                                                                                                                                                                                           gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTct                                                          <  2:1520/251‑2 (MQ=255)                                                                                                                                                                                               aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTcggcgg                                                      <  1:302398/251‑1 (MQ=255)                                                                                                                                                                                                  tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAGGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGACCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGCGGGGGACGGCGGTTc                                                   >  2:127067/1‑251 (MQ=255)                                                                                                                                                                                                  tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGtt                                                    <  1:99173/250‑1 (MQ=255)                                                                                                                                                                                                               gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGtcatcc                                      >  2:286789/1‑250 (MQ=255)                                                                                                                                                                                                                      aagcAATTGGTCGCAGGTACGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACAGAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCAGCACccacca                               <  2:520829/248‑1 (MQ=255)                                                                                                                                                                                                                       aTCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCcaccac                              <  1:28241/251‑1 (MQ=255)                                                                                                                                                                                                                              ggTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCCTCACCCACCACCggggcg                       >  1:72468/1‑248 (MQ=255)                                                                                                                                                                                                                                      gTTCGAATCCTGCACGACTCACCAATCGCTAAG‑GTGGACGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTGGTAGTATCCAGCGTAGTATCGGGAGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTTGATCCCCTCATCACCCACCCCCGGGGCGTTAGCTcc               >  2:90279/1‑250 (MQ=37)                                                                                                                                                                                                                                                       cccACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAg  <  1:128266/247‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       AACGCGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGG‑GTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAG‑GTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAG  >  NC_000913/780378‑780869

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