Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,116,273	G→A	9.8%	T520I (ACC→ATC)	valS ←	valine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,116,273	0	G	A	9.8%	165.5 / 7.2	61	T520I (ACC→ATC)	valS	valine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (25/30);  new base A (3/3);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

CCGTCGACGATGTCCAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCCACAGGCCCGAGCTGAACCAGGTGT  >  NC_002947/1116129‑1116420                                                                                                                                                 |                                                                                                                                                     ccGTCGACGATGTCCAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCa                                                                                                                                                   <  2:249089/148‑1 (MQ=255)      gacgaTGTCCAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATg                                                                                                                                              <  6:313499/148‑1 (MQ=255)      gacgaTGTCCAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGa                                                                                                                                             >  2:64034/1‑149 (MQ=255)         gaTGTCCAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCa                                                                                                                                          <  4:166079/149‑1 (MQ=255)         gaTGTCCAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCa                                                                                                                                          >  4:185487/1‑149 (MQ=255)             tCCAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGGGCATGGTCAGCATGATCATgcg                                                                                                                                      <  3:317005/149‑1 (MQ=255)               cAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGATCATGCTGACCATgcac                                                                                                                                     <  5:350254/148‑3 (MQ=255)               cAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGATCATGCTGACCATgcac                                                                                                                                     <  5:344614/148‑3 (MQ=255)               cAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGATCATGCTGACCATgcac                                                                                                                                     >  6:344614/1‑146 (MQ=255)               cAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGATCATGCTGACCATgcac                                                                                                                                     >  6:350254/1‑146 (MQ=255)               cAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGGGCCGTCTTCGTTCTTGATCAGGTGCATGATCATGCTGACCATgcacc                                                                                                                                    <  4:234312/149‑4 (MQ=255)                aGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATgcgcgc                                                                                                                                   <  5:186231/149‑1 (MQ=255)                   gggTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAAc                                                                                                                                <  5:330829/149‑1 (MQ=255)                    ggTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGCGCATGATCATGCTGACCATgcacctgat                                                                                                                                >  5:231117/1‑141 (MQ=255)                      tCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAAcc                                                                                                                               >  3:173518/1‑147 (MQ=255)                         aGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCagaa                                                                                                                          <  1:53517/149‑1 (MQ=255)                         aGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCagaa                                                                                                                          <  3:55703/149‑1 (MQ=255)                              gTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCagaag                                                                                                                         >  2:221390/1‑145 (MQ=255)                               ttGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTTCACGTACTCGGTCTTGAACGGTACCTGCGGGGGGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAgat                                                                                                                    <  7:59060/149‑1 (MQ=255)                               ttGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAgat                                                                                                                    <  5:356956/149‑1 (MQ=255)                               ttGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAgat                                                                                                                    <  4:270883/149‑1 (MQ=255)                               ttGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAgat                                                                                                                    <  1:167005/149‑1 (MQ=255)                                   ccTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAgatgat                                                                                                                 <  7:98718/148‑1 (MQ=255)                                             ggacATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAg                                                                                                      >  6:144284/1‑149 (MQ=255)                                               acATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGt                                                                                                     >  3:170414/1‑148 (MQ=255)                                                 aTCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCAc                                                                                                  <  6:238833/149‑1 (MQ=255)                                                         gccctggccGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACgc                                                                                          >  5:3087/1‑148 (MQ=255)                                                           cctggccGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGc                                                                                                         >  6:264844/1‑132 (MQ=255)                                                           cctggccGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGc                                                                                                         <  5:264844/132‑1 (MQ=255)                                                             tggccGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACGCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGgt                                                                                      >  5:321891/1‑149 (MQ=255)                                                               gccGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGgtg                                                                                     <  2:83314/148‑1 (MQ=255)                                                                   tCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGgtggagtg                                                                                >  8:131406/1‑149 (MQ=255)                                                                    cACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGgtggagtg                                                                                >  7:49056/1‑148 (MQ=255)                                                                            aGGCCGTGCACGTACACGGTCTTGAACGGTAGCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAActtctt                                                                       <  4:170414/149‑1 (MQ=255)                                                                             ggCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCAc                                                                                                  >  2:338519/1‑121 (MQ=255)                                                                             ggCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCAc                                                                                                  <  1:338519/121‑1 (MQ=255)                                                                              gCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGa                                                                     <  2:75239/149‑1 (MQ=255)                                                                                      cGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCg                                                             <  7:131406/149‑1 (MQ=255)                                                                                          cacGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCt                                                         >  2:148330/1‑149 (MQ=255)                                                                                          cacGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGCCGGTGGAGTGGAACTTCTTGAGGAATTCGGTc                                                          <  6:3087/148‑1 (MQ=255)                                                                                                                ggggTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGga                                  <  3:185487/149‑1 (MQ=255)                                                                                                                    tGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAAt                              <  2:11303/149‑1 (MQ=255)                                                                                                                      ccGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGt                            >  1:200378/1‑149 (MQ=255)                                                                                                                      ccGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGt                            >  2:35162/1‑149 (MQ=255)                                                                                                                        gTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTc                           >  8:319894/1‑148 (MQ=255)                                                                                                                        gTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTc                           >  3:247172/1‑148 (MQ=255)                                                                                                                         tCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCca                         >  1:282284/1‑149 (MQ=255)                                                                                                                         tCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCca                         >  6:219412/1‑149 (MQ=255)                                                                                                                         tCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCca                         >  2:47259/1‑149 (MQ=255)                                                                                                                             cGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAATATGATGGCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCcaca                       <  7:274781/147‑1 (MQ=255)                                                                                                                               ttCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCGGGCC‑AGCCCAGGGTGGAGAATGTCCACAGGcc                   >  4:85967/1‑149 (MQ=255)                                                                                                                                 cTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCCACAGGccc                  <  1:47259/148‑1 (MQ=255)                                                                                                                                 cTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCCACAGGccc                  <  4:173518/148‑1 (MQ=255)                                                                                                                                  ttGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGGTGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCCACAGGccc                  >  7:8977/1‑147 (MQ=255)                                                                                                                                  ttGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCCACAGGccc                  >  8:140008/1‑147 (MQ=255)                                                                                                                                    gATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAATCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCCAAGCCCAGGGTGGAGAATGTCCACAGTccc                  <  1:221390/146‑1 (MQ=255)                                                                                                                                    gATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCCACAGGCCCGAg               <  1:148330/148‑1 (MQ=255)                                                                                                                                           tGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGgaga                                <  3:219024/124‑1 (MQ=255)                                                                                                                                           tGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGgaga                                >  4:219024/1‑124 (MQ=255)                                                                                                                                            gCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCCACAGGCCCGAGCTGAACCAg      >  7:45760/1‑149 (MQ=255)                                                                                                                                             caaggTCATCAAGATCATGCGCGAAACACAGAAGAAGATGATGTCGAAGAAAGACACAAGCACGACGGTCGACAGGAACATCAAGATTAATTCGGTCTCATACCTCC‑AGCCCAGCGTGCAGAAAGTCCACAGGCCCGAGCTGAACCAg      >  5:232795/4‑148 (MQ=255)                                                                                                                                             cATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCCACAGGCCCGAGCTGAACCAg      <  6:321891/148‑1 (MQ=255)                                                                                                                                                ggTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCCACAGGCCCGAGCTGAACCAGgtgt  >  8:183061/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CCGTCGACGATGTCCAGCGGGTCCAGCACGTTGCCCTTGGACTTGGACATCTTCTGGCCCTGGCCGTCACGTACCAGGCCGTGCACGTACACGGTCTTGAACGGTACCTGCGGGGTGCCGTCTTCGTTCTTGATCAGGTGCATGGTCAGCATGATCATGCGCGCAACCCAGAAGAAGATGATGTCGAAGCCAGTCACCAGCACGTCGGTGGAGTGGAACTTCTTGAGGAATTCGGTCTGCTCCGGCC‑AGCCCAGGGTGGAGAATGTCCACAGGCCCGAGCTGAACCAGGTGT  >  NC_002947/1116129‑1116420

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