Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,102,752	C→G	5.7%	noncoding (35/76 nt)	PP_t21 →	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,102,752	0	C	G	5.7%	220.9 / 3.5	70	noncoding (35/76 nt)	PP_t21	tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base C (29/37);  new base G (2/2);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

GCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTC  >  NC_002947/2102604‑2102897                                                                                                                                                     |                                                                                                                                                  gCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                                                                                                                                   >  5:242638/1‑149 (MQ=255)     ggAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGTGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAg                                                                                                                                               <  3:166241/149‑1 (MQ=255)           aaaGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgt                                                                                                                                          <  2:279101/148‑2 (MQ=255)           aaaGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                          <  3:48931/148‑1 (MQ=255)           aaaGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                          <  6:242638/148‑1 (MQ=255)           aaaGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACACCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                         <  1:325939/149‑1 (MQ=255)            aaGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                         <  7:74700/148‑1 (MQ=255)               caggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcg                                                                                                                                     >  8:54170/1‑149 (MQ=255)                     gTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGtt                                                                                                                               <  6:343137/149‑1 (MQ=255)                        aTGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                            <  2:111033/149‑1 (MQ=255)                           ccAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                                                                                                                          <  2:354835/148‑1 (MQ=40)                                  cGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                                                                  <  2:90715/149‑1 (MQ=34)                                   gtgtCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                                                                                                 <  3:288836/149‑1 (MQ=33)                                     gtCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCg                                                                                                               >  6:185771/1‑149 (MQ=33)                                          cTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                                                            >  5:125516/1‑147 (MQ=33)                                               tCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                            <  6:202996/126‑1 (MQ=34)                                               tCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                            >  5:202996/1‑126 (MQ=34)                                                   gTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTtctc                                                                                                   <  5:185771/147‑1 (MQ=35)                                                   gTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCg                                                                                                 >  7:337781/1‑149 (MQ=38)                                                   gTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGCGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCg                                                                                                 >  6:340027/1‑149 (MQ=38)                                                      gCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCa                                                                                              >  2:188018/1‑149 (MQ=40)                                                        cTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACt                                                                                            <  5:243381/149‑1 (MQ=255)                                                         tAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTg                                                                                           >  5:6706/1‑149 (MQ=255)                                                          aGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGc                                                                                          >  1:202051/1‑149 (MQ=255)                                                                 cGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGt                                                                                                                                       <  7:242957/97‑1 (MQ=35)                                                                 cGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGt                                                                                                                                       >  8:242957/1‑97 (MQ=35)                                                                 cGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                                    >  7:151475/1‑100 (MQ=35)                                                                 cGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                                    <  8:151475/100‑1 (MQ=35)                                                                 cGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAgcgg                                                                                   >  7:333939/1‑149 (MQ=255)                                                                  gCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAgcggc                                                                                  >  8:133072/1‑149 (MQ=255)                                                                     cTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCggg                                                                               <  3:213574/149‑1 (MQ=255)                                                                      tttCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGACgggg                                                                              >  3:108394/1‑149 (MQ=255)                                                                        tCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCgggg                                                                              <  6:125516/147‑1 (MQ=255)                                                                          aCGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAg                                                                           >  6:89510/1‑148 (MQ=255)                                                                          aCGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAg                                                                           <  7:295165/148‑1 (MQ=255)                                                                             gcggTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTtct                                                                                                    >  2:312424/1‑120 (MQ=37)                                                                              cggTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaaa                                                                      <  7:133072/149‑1 (MQ=255)                                                                                gTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACTACCTTGCCAAGTTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaaaga                                                                    <  5:340027/149‑1 (MQ=255)                                                                                 tAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaaagaa                                                                   <  7:54170/149‑1 (MQ=255)                                                                                   aCAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaaagaaaa                                                                 >  5:372883/1‑149 (MQ=255)                                                                                   aCAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaaagaaaa                                                                 <  6:372883/149‑1 (MQ=255)                                                                                    cAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAAc                                                                >  5:93457/1‑149 (MQ=255)                                                                                          ttCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGc                                                          <  5:89510/149‑1 (MQ=255)                                                                                             gAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTg                                                        >  2:55968/1‑148 (MQ=255)                                                                                              aGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGACGGGGCAGGAAAAGAAAACCCAGGCTGaa                                                      <  4:108394/149‑1 (MQ=255)                                                                                                  cccTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCgggg                                                                              <  3:249562/121‑1 (MQ=255)                                                                                                  cccTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCgggg                                                                              >  4:249562/1‑121 (MQ=255)                                                                                                    cTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGa                                                <  8:337781/149‑1 (MQ=255)                                                                                                       ggggACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAgg                                              <  2:202051/148‑1 (MQ=255)                                                                                                        gggACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGtt                                            >  7:280301/1‑149 (MQ=255)                                                                                                        gggACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGt                                             <  8:333939/148‑1 (MQ=255)                                                                                                         ggACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGAGAGGtt                                            >  6:222899/1‑148 (MQ=255)                                                                                                             gCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCa                                          <  4:220771/146‑1 (MQ=255)                                                                                                             gCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCa                                          >  3:220771/1‑146 (MQ=255)                                                                                                             gCCAATGCGGGAATAGCTCACTTGGTAGAGCACGACCTTGCCAAGGTCGGCGTCGCGAGATCGAGTATCGTTTCCCGCACCAGTATATCCGGCAATGCTTAGCGGCGGGGAACGAAAAGAAAACCCAGCCTGAATCGTGAGGTTCAtct                                       >  2:308612/1‑147 (MQ=21)                                                                                                                        aaTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACt                            <  6:6706/149‑1 (MQ=255)                                                                                                                            gCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTgg                                  <  4:45892/139‑1 (MQ=255)                                                                                                                            gCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTgg                                  >  3:45892/1‑139 (MQ=255)                                                                                                                                  ttGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGTAACTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgtt                  <  1:308612/149‑1 (MQ=255)                                                                                                                                  ttGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgtt                  <  1:188018/149‑1 (MQ=255)                                                                                                                                     gTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCTTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTTAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGtgttgt                <  1:55968/148‑1 (MQ=255)                                                                                                                                      tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTtgtgt              <  8:280301/149‑1 (MQ=255)                                                                                                                                         agCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTc           >  3:298222/1‑149 (MQ=255)                                                                                                                                            ccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCa                              >  3:193512/2‑127 (MQ=255)                                                                                                                                            ccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCa                              <  4:193512/126‑1 (MQ=255)                                                                                                                                            aCGACCGTGCCTAGCTCGGGGTCGCGAGTTCGAGTATCGTTTCCAGCTCCAGTTTCTCCTGCACTCCTTAGCCGCGGGGCAGGAAAAGAAAACCGAGGGTGAATCGGGAGCTACATGTTCTGGTGCAATGAAAAGTCTtctggccgca         >  2:245309/1‑142 (MQ=21)                                                                                                                                              gACCTTGCCAAGGTCGGGGTCGCTAGTTCGACTCTAGTTTAACGCTACATTTTCTCAGGCACTGCTTAGCGGCGGGGCAGGAAAATAAAACAAAGTCTGAATCGTGAGGTTCATGTTCTg                                   >  6:378673/1‑120 (MQ=25)                                                                                                                                              gACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTg                                   <  5:378673/120‑1 (MQ=255)                                                                                                                                                ccTTGGCAAGGTCGTGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCg    >  7:254944/1‑149 (MQ=255)                                                                                                                                                 cTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGt   <  6:93457/149‑1 (MQ=255)                                                                                                                                                  ttGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTc  >  3:9367/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCATGTTCTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTC  >  NC_002947/2102604‑2102897

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