Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,102,711	C→T	6.5%	noncoding (72/76 nt)	PP_t20 →	tRNA‑Glu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,102,711	0	C	T	6.5%	180.3 / 2.2	61	noncoding (72/76 nt)	PP_t20	tRNA‑Glu
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base T (2/2);  total (29/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 = 7.94e-01

TGCGCGCCACTTGCAGCGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCA  >  NC_002947/2102565‑2102857                                                                                                                                                   |                                                                                                                                                   tGCGCGCCACTTGCAGCGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGc                                                                                                                                                  >  1:117275/1‑149 (MQ=255) tGCGCGCCACTTGCAGCGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGc                                                                                                                                                  <  6:271576/149‑1 (MQ=255) tGCGCGCCACTTGCAGCGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGc                                                                                                                                                  >  5:309965/1‑149 (MQ=255)       ccACTTGCAGCGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGc                                                                                                                                            <  7:61428/149‑1 (MQ=255)       ccACTTGCAGCGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGc                                                                                                                                            >  7:118959/1‑149 (MQ=255)         aCTTGCAGCGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCgg                                                                                                                                          >  5:135805/1‑149 (MQ=255)         aCTTGCAGCGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCgg                                                                                                                                          <  7:248726/149‑1 (MQ=255)         aCTTGCAGCGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCgg                                                                                                                                          <  5:260179/149‑1 (MQ=255)               aGCGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAg                                                                                                                                    <  4:295323/149‑1 (MQ=255)                 cGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGa                                                                                                                                        <  6:20887/143‑1 (MQ=255)                 cGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGa                                                                                                                                        >  5:20887/1‑143 (MQ=255)                 cGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCt                                                                                                                                  <  6:309965/149‑1 (MQ=255)                           acaGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTgg                                                                                                                           >  2:151179/1‑146 (MQ=255)                           acaGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTgg                                                                                                                           >  4:204090/1‑146 (MQ=255)                              gCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGc                                                                                                                     <  3:204090/149‑1 (MQ=255)                                        gCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                                                                                           <  8:118959/149‑1 (MQ=255)                                         cAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                                                                                           >  7:118834/1‑148 (MQ=255)                                            ggAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAg                                                                                                       >  4:163683/1‑149 (MQ=255)                                                tGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATa                                                                                                                                     <  2:147475/115‑1 (MQ=255)                                                tGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATa                                                                                                                                     >  1:147475/1‑115 (MQ=255)                                                  aaaGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgt                                                                                                  <  8:96808/148‑2 (MQ=255)                                                    aGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                                                                                                <  2:117275/148‑1 (MQ=255)                                                     gcaggcagTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgtgctc                                                                                              <  4:133413/149‑6 (MQ=255)                                                             tAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                                      >  4:197049/1‑149 (MQ=255)                                                               aTGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGTATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                    <  1:76663/149‑1 (MQ=255)                                                               aTGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                    <  3:163683/149‑1 (MQ=255)                                                                    agcagcGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                               <  8:118834/149‑1 (MQ=38)                                                                      cagcGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGt                                                                             >  3:177822/1‑149 (MQ=37)                                                                               cccTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCt                                                                     >  3:79849/1‑148 (MQ=33)                                                                                 cTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                    >  5:182773/1‑147 (MQ=33)                                                                                 cTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCa                                                                  <  3:103132/149‑1 (MQ=33)                                                                                      tCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGt                                                                >  7:162617/1‑146 (MQ=33)                                                                                      tCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGt                                                                <  8:162617/146‑1 (MQ=33)                                                                                      tCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTtc                                                             >  1:125830/1‑149 (MQ=34)                                                                                          gTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCg                                                         <  4:177822/149‑1 (MQ=38)                                                                                              ccTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGATGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                 >  7:168356/1‑105 (MQ=25)                                                                                              ccTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGATGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                 <  8:163916/105‑1 (MQ=25)                                                                                              ccTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGATGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                 <  8:168356/105‑1 (MQ=25)                                                                                              ccTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGATGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAAGTCggg                                                                                                 >  7:163916/1‑105 (MQ=255)                                                                                                  ggACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATTCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGtt                                                                            <  5:93439/122‑1 (MQ=21)                                                                                                  ggACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGtt                                                                            >  6:93439/1‑122 (MQ=34)                                                                                                     cacCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                          <  2:72458/121‑1 (MQ=34)                                                                                                     cacCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                          >  1:72458/1‑121 (MQ=34)                                                                                                      acCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTtctc                                                           >  2:276058/1‑135 (MQ=37)                                                                                                      acCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTtctc                                                           <  1:276058/135‑1 (MQ=37)                                                                                                         gCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAgcggc                                          >  1:271629/1‑149 (MQ=255)                                                                                                         gCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAgcggc                                          <  2:125830/149‑1 (MQ=255)                                                                                                           ccTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGcggcgg                                        <  6:182773/149‑1 (MQ=255)                                                                                                           ccTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGcggcgg                                        <  6:135805/149‑1 (MQ=255)                                                                                                           ccTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGcggcgg                                        <  5:288332/149‑1 (MQ=255)                                                                                                           ccTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGcggcgg                                        <  3:197049/149‑1 (MQ=255)                                                                                                             tttCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCgggg                                      >  2:286732/1‑149 (MQ=255)                                                                                                                 aCGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGGTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAg                                   >  2:230852/1‑148 (MQ=255)                                                                                                                          aCAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGgaaaagaaaa                         <  1:180815/149‑1 (MQ=255)                                                                                                                                    gcgTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGa               <  2:243299/147‑1 (MQ=255)                                                                                                                                     aGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGaa              >  4:196642/1‑149 (MQ=255)                                                                                                                                      gTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGaa              >  4:21054/1‑148 (MQ=255)                                                                                                                                       tCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATc            <  1:151179/149‑1 (MQ=255)                                                                                                                                          ccTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGt          >  3:107387/1‑148 (MQ=255)                                                                                                                                                 gACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCa  >  7:13625/1‑149 (MQ=255)                                                                                                                                                 gACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCa  >  7:161717/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TGCGCGCCACTTGCAGCGTAAAGCACACAGCGAAACGCGGCAGGGAGTGAAAGCAGGCAGTAATGCCAGCAGCGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTCTCCGGCACTGCTTAGCGGCGGGGCAGGAAAAGAAAACCCAGGCTGAATCGTGAGGTTCA  >  NC_002947/2102565‑2102857

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