Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,245,444	G→A	6.5%	Y277Y (TAC→TAT)	alr ←	alanine racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,245,444	0	G	A	6.5%	180.6 / 2.6	61	Y277Y (TAC→TAT)	alr	alanine racemase
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CCCGGCGGTAGCCATCGGAGTACCCGACCGTAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCGGTCCAGCCTGGCGTGCTTGATC  >  NC_002947/4245296‑4245591                                                                                                                                                     |                                                                                                                                                    cccGGCGGTAGCCATCGGAGTACCCGACCGTAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTg                                                                                                                                                     >  3:97427/1‑149 (MQ=255)  ccGGCGGTAGCCATCGGAGTACCCGACCGTAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTg                                                                                                                                                     >  1:99198/1‑148 (MQ=255)                  gAGTACCCGACCGTAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCgg                                                                                                                                    >  5:66324/1‑149 (MQ=255)                       cccGACCGTAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACgg                                                                                                                               >  8:114491/1‑149 (MQ=255)                       cccGACCGTAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACgg                                                                                                                               <  6:324802/149‑1 (MQ=255)                       cccGACCGTAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCGGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACgg                                                                                                                               <  5:307633/149‑1 (MQ=255)                        ccGACCGTAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGgt                                                                                                                              >  1:89384/1‑149 (MQ=255)                        ccGACCGTAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGCGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGgt                                                                                                                              <  3:145824/149‑1 (MQ=255)                        ccGACAGTAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCTGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCTCCTCCACGTGCGATTTGAAATGCATCGCACGTTTGTACTCGGTGCGAGCCGGAACGgt                                                                                                                              >  4:54656/1‑149 (MQ=255)                           aCCGTAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTc                                                                                                                           <  4:97427/149‑1 (MQ=255)                               tAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGTGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGcc                                                                                                                        >  5:166272/1‑148 (MQ=255)                               tAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGcc                                                                                                                        >  7:46302/1‑148 (MQ=255)                                  tGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAACGTGCGGTCATAGACCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTAGGTGCGCGCCGGCACGGTGTCGCCGaa                                                                                                                     >  8:116910/1‑148 (MQ=255)                                   gTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGaa                                                                                                                     >  7:131820/1‑147 (MQ=255)                                       gCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAgcgc                                                                                                               <  8:308700/149‑1 (MQ=255)                                            ccGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCAc                                                                                                          <  2:99198/149‑1 (MQ=255)                                            ccGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCAc                                                                                                          >  3:261990/1‑149 (MQ=255)                                               cgAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCg                                                                                                       >  6:352279/1‑149 (MQ=255)                                                 aaTCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGtt                                                                                                     <  4:261990/149‑1 (MQ=255)                                                    cACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCg                                                                                                   >  8:266420/1‑148 (MQ=255)                                                     aCGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTa                                                                                                 >  3:316146/1‑149 (MQ=255)                                                      cGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTAc                                                                                                <  8:218981/149‑1 (MQ=255)                                                       gggCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTAcc                                                                                               <  7:114491/149‑1 (MQ=255)                                                                  aaGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGgc                                                                                     <  6:329210/148‑1 (MQ=255)                                                                    ggTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTgcgcgc                                                                                                                                       <  7:253211/96‑1 (MQ=255)                                                                    ggTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTgcgcgc                                                                                                                                       >  8:253211/1‑96 (MQ=255)                                                                    ggTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTATACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGgcgcg                                                                                  >  4:195811/1‑149 (MQ=255)                                                                    ggTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTATACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGgcgcg                                                                                  >  4:319299/1‑149 (MQ=255)                                                                     gTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGgcgcgc                                                                                 >  6:96984/1‑149 (MQ=255)                                                                           tCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCgg                                                                           <  6:66324/149‑1 (MQ=255)                                                                           tCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCgg                                                                           >  6:139197/1‑149 (MQ=255)                                                                           tCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCgg                                                                           >  1:264717/1‑149 (MQ=255)                                                                             aTAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGc                                                                          <  8:46302/148‑1 (MQ=255)                                                                             aTAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCa                                                                         >  2:282863/1‑149 (MQ=255)                                                                                   cACGGTGTTGCCGGCCGGATAGCTGTGCACCGCAGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCTCCGGCACGGTTTCGCCGAACAGAGCGCCACCCGTTCGTACCATGTCCAGGAGCGCTTCCGGCACTTc                                                                     >  4:20425/1‑147 (MQ=255)                                                                                         gtTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGGACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAg                                                             >  7:15932/1‑149 (MQ=255)                                                                                            gccggccgGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCg                                                           >  5:8807/1‑148 (MQ=255)                                                                                              cggccggATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAAc                                                        <  5:96984/149‑1 (MQ=255)                                                                                                gccggATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGa                                                      <  6:166272/149‑1 (MQ=255)                                                                                                gccggATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGa                                                      <  4:316146/149‑1 (MQ=255)                                                                                                 ccggATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAg                                                     >  2:52417/1‑149 (MQ=255)                                                                                                       aGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTgg                                                 <  6:136336/147‑1 (MQ=255)                                                                                                       aGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTgg                                                 >  5:136336/1‑147 (MQ=255)                                                                                                            tgCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTg                                          >  2:34471/1‑149 (MQ=255)                                                                                                               aCCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAg                                       >  5:90244/1‑149 (MQ=255)                                                                                                               aCCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAg                                       <  7:266420/149‑1 (MQ=255)                                                                                                                  gccgccACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGt                                    <  7:116910/149‑1 (MQ=255)                                                                                                                   ccgccACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGt                                    >  1:232581/1‑148 (MQ=255)                                                                                                                   ccgccACGTGCGATTTGAACTGCATCGCACGTTTATACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTg                                   <  3:195811/149‑1 (MQ=255)                                                                                                                   ccgccACGTGCGATTTGAACTGCATCGCACGTTTATACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTg                                   <  3:319299/149‑1 (MQ=255)                                                                                                                    cgccACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGa                                  <  5:139197/149‑1 (MQ=255)                                                                                                                              gATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCg                        <  5:352279/149‑1 (MQ=255)                                                                                                                              gATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGGTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCg                        <  4:290334/149‑1 (MQ=255)                                                                                                                               aTTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCgg                       <  2:89384/149‑1 (MQ=255)                                                                                                                               aTTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCgg                       <  1:282863/149‑1 (MQ=255)                                                                                                                               aTTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCgg                       >  8:324505/1‑149 (MQ=255)                                                                                                                                tttGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGGACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCGGt                      <  8:15932/149‑1 (MQ=255)                                                                                                                                     aCTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGGGCCACCCGTTCGTACCATGTCCAGGCGCGGTTCCGGCTCTTGCAGCGTAGAGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCGGTCCAGc                 <  3:54656/149‑1 (MQ=255)                                                                                                                                          aTCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCGGTCCAGCCTGGc            >  7:288062/1‑149 (MQ=255)                                                                                                                                              cACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCGGTCCAGCCTGGCGTGc        <  6:8807/149‑1 (MQ=255)                                                                                                                                              cACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCGGTCCAGCCTGGCGTGc        <  8:131820/149‑1 (MQ=255)                                                                                                                                                 gTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCGGTCCAGCCTGGCGTGCtt      >  2:22682/1‑148 (MQ=255)                                                                                                                                                  tttGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCGGTCCAGCCTGGCGTGCTTGa    <  1:52417/149‑1 (MQ=255)                                                                                                                                                     gTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGGTTGCTGCGGTCCAGCCTGGCGTGCTTGATc  <  1:34471/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCCGGCGGTAGCCATCGGAGTACCCGACCGTAATGTTGGCCAGCCGCGAATCACGGGCCAGGGTGAAGGTGCGGTCATAGCCCACGGTGTTGCCGGCCGGATAGCTGTGCACCGCCGCCACGTGCGATTTGAACTGCATCGCACGTTTGTACTCGGTGCGCGCCGGCACGGTGTCGCCGAACAGCGCGCCACCCGTTCGTACCATGTCCAGGCGCGCTTCCGGCACTTCCAGCGTAGCGAACGAGTTGGCGGCGTGCAGGGTGAGCTTGCTGCGGTCCAGCCTGGCGTGCTTGATC  >  NC_002947/4245296‑4245591

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