Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,428,065	A→T	6.3%	S672C (AGC→TGC)	PP_3042 →	terminase large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,428,065	0	A	T	6.3%	192.5 / 2.8	63	S672C (AGC→TGC)	PP_3042	terminase large subunit
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base T (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

CCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACTGATCAAGGCCCG  >  NC_002947/3427918‑3428211                                                                                                                                                    |                                                                                                                                                   ccGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCcagcag                                                                                                                                                   >  6:3136/1‑149 (MQ=255)  cGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCagcagc                                                                                                                                                  <  8:145040/149‑1 (MQ=255)   ggTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCg                                                                                                                                                 >  4:271405/1‑149 (MQ=255)       gTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCCGCACCTCCCCCCATGCAACGCCGCAGCTCCAGCAGCGGCTa                                                                                                                                             >  5:4069/1‑149 (MQ=255)        ttGAGCACCAACACGATGACGACGATGACCCGGTCGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCTGCTGGAGCTGCGGcgtt                                                                                                                                            <  5:344110/149‑4 (MQ=255)            gCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTATCACCGCAAACACCACCGGCTCGGCCCGCTAAACGCTCCACTCCACCCCCGCCTCGCCCGGCACCTCCCCCCAAGCAACCCCGCAGCTCCAGCAGCGGCTATcag                                                                                                                                         >  4:38332/1‑146 (MQ=255)                cAACACGATGACGACGATGACCCGGACGACGATGCAGTTTCACCGCCAACAGCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGTCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAagag                                                                                                                                    <  3:38332/149‑1 (MQ=255)                    acGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCagcagc                                                                                                                                                  >  4:182680/1‑131 (MQ=255)                    acGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCagcagc                                                                                                                                                  <  3:182680/131‑1 (MQ=255)                      gATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCt                                                                                                                               >  5:259499/1‑148 (MQ=255)                           cgacgaTGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTg                                                                                                                              <  4:202410/144‑1 (MQ=255)                           cgacgaTGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTg                                                                                                                              >  3:202410/1‑144 (MQ=255)                                aTGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCata                                                                                                                    >  2:7826/1‑149 (MQ=255)                                     ccGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATAccc                                                                                                               <  1:7826/149‑1 (MQ=255)                                       ggACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCt                                                                                                                               <  7:328650/131‑1 (MQ=255)                                       ggACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCt                                                                                                                               >  8:328650/1‑131 (MQ=255)                                             cgaTCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCacac                                                                                                       <  3:58059/149‑1 (MQ=255)                                             agaTCCAGTTTCACCGGCAACACCCCCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCATCGCCGCAGCTCCAGCAGCGGCTATCTGTAGAGACGGTGACATGGCATATACCCATGCacac                                                                                                       <  6:4069/148‑1 (MQ=255)                                                 ccAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTcg                                                                                                   >  5:147093/1‑149 (MQ=255)                                                 ccAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTcg                                                                                                   >  8:86887/1‑149 (MQ=255)                                                  cAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATCCCCTGGCACACCTcgc                                                                                                  >  1:123720/1‑149 (MQ=255)                                                   aggtgcACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTcgcg                                                                                                 <  8:230442/144‑1 (MQ=255)                                                     tttCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGc                                                                                                           <  6:10090/137‑1 (MQ=255)                                                     tttCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGc                                                                                                           >  5:10090/1‑137 (MQ=255)                                                     tttCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACCCCTCGCGGc                                                                                               >  2:334143/1‑149 (MQ=255)                                                     tttCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGc                                                                                               >  4:209825/1‑149 (MQ=255)                                                       tCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGcc                                                                                              <  5:3136/148‑1 (MQ=255)                                                          ccGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTc                                                                                           >  2:206252/1‑148 (MQ=255)                                                                    caccGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGa                                                                                 >  3:176843/1‑148 (MQ=255)                                                                    caccGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGa                                                                                 <  4:176843/148‑1 (MQ=255)                                                                     accGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGAtt                                                                               <  8:105534/149‑1 (MQ=255)                                                                      ccGGCTCGGCCCGCTCCACGCTCCACTCCTCCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTg                                                                              <  3:271405/149‑1 (MQ=255)                                                                      ccGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTg                                                                              >  2:283497/1‑149 (MQ=255)                                                                                ccGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGaa                                                                                                                                       >  6:128251/1‑82 (MQ=255)                                                                                ccGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGaa                                                                                                                                       <  5:128251/82‑1 (MQ=255)                                                                                      cacgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGc                                                                      <  5:102894/141‑1 (MQ=255)                                                                                      cacgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGc                                                                      >  6:102894/1‑141 (MQ=255)                                                                                      cacgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGgt                                                              >  5:149543/1‑149 (MQ=255)                                                                                       acgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGgtt                                                             >  1:156166/1‑149 (MQ=255)                                                                                         gctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGCCATAAACCCAGGCACACCACCCGGCCGTCGAGCGTGCGAATGCGCGTGGCTAACGGgttgt                                                           >  8:315301/1‑149 (MQ=255)                                                                                          ctcctctccaCCGCGGCCTCGCCCGGAACCTCCACCCATTCAACGCCGCAGCTCCAGCAGAGGCTATCTGAAGAGACGCTGACATGGCATATAGCCAGGCAGACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGgttgtt                                                          <  8:149971/149‑1 (MQ=255)                                                                                          ctccactccaCCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGgttgtt                                                          <  3:209825/149‑1 (MQ=255)                                                                                           tccactccacCGCCGCCTCGCCCGGCACCTCCGCCCATGCAACGCCGCAGCTCCAGCAGCGGCCATCTGAAGAGACGCTGACATCGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGATGTTc                                                         >  3:199901/1‑149 (MQ=255)                                                                                             cactccacCGCCGCCTCGCCCGGCACCTCCACCCCTGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGa                                                       <  2:326359/149‑1 (MQ=255)                                                                                                tccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCACGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTACGAtc                                                     >  8:121466/1‑147 (MQ=255)                                                                                                  cacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGGCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATAc                                                    <  6:259499/147‑1 (MQ=255)                                                                                                   acCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGc                                                 <  1:283497/149‑1 (MQ=255)                                                                                                      gccgccTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGAc                                              >  7:12/1‑149 (MQ=255)                                                                                                             cGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGgttgtt                                                          >  5:334825/1‑130 (MQ=255)                                                                                                             cGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGgttgtt                                                          <  6:334825/130‑1 (MQ=255)                                                                                                              gCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCg                                       <  7:315301/148‑1 (MQ=255)                                                                                                                 cGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGa                                   <  8:12/149‑1 (MQ=255)                                                                                                                     aCCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTAt                               <  2:156166/149‑1 (MQ=255)                                                                                                                          cACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGAc                          <  6:147093/149‑1 (MQ=255)                                                                                                                           aCCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGTCGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGCCg                         >  3:208844/1‑149 (MQ=255)                                                                                                                           aCCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGa                                   <  1:4908/139‑1 (MQ=255)                                                                                                                           aCCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGa                                   >  2:4908/1‑139 (MQ=255)                                                                                                                             ccATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACgg                        >  3:190326/1‑148 (MQ=255)                                                                                                                             ccATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGt                       <  7:86887/149‑1 (MQ=255)                                                                                                                                tGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGa                    <  8:132791/149‑1 (MQ=255)                                                                                                                                   aaCGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACg                  >  8:286148/1‑148 (MQ=255)                                                                                                                                      gccgcTGCTGGAGCTGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGa                                                       <  5:158464/108‑1 (MQ=255)                                                                                                                                      gccgcTGCTGGAGCTGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGa                                                       >  6:158464/1‑108 (MQ=255)                                                                                                                                       ccgcTGCTGGAGCTGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACt              >  2:181850/1‑148 (MQ=255)                                                                                                                                        cgcAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAATATCGGa                           <  1:227688/134‑1 (MQ=255)                                                                                                                                        cgcAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAATATCGGa                           >  2:227688/1‑134 (MQ=255)                                                                                                                                           aGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACTGATc          >  5:333573/1‑148 (MQ=255)                                                                                                                                                  gcagcGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACTGATCAAGGCCcg  >  1:15846/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACTGATCAAGGCCCG  >  NC_002947/3427918‑3428211

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