Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,428,056	A→T	6.0%	S669C (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,056	0	A	T	6.0%	198.8 / 2.4	67	S669C (AGC→TGC)	PP_3042	terminase large subunit
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base T (2/2);  total (34/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 = 7.61e-01

CGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACTGATC  >  NC_002947/3427910‑3428203                                                                                                                                                   |                                                                                                                                                    cGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGc                                                                                                                                                   <  5:349505/149‑1 (MQ=255)         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                                                                                                                                          <  3:182680/131‑1 (MQ=255)                            acGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCagcagc                                                                                                                                          >  4:182680/1‑131 (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                                                                                           >  8:86887/1‑149 (MQ=255)                                                         ccAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTcg                                                                                           >  5:147093/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                                                                         <  4:176843/148‑1 (MQ=255)                                                                            caccGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGa                                                                         >  3:176843/1‑148 (MQ=255)                                                                             accGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGAtt                                                                       <  8:105534/149‑1 (MQ=255)                                                                              ccGGCTCGGCCCGCTCCACGCTCCACTCCTCCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTg                                                                      <  3:271405/149‑1 (MQ=255)                                                                              ccGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTg                                                                      >  2:283497/1‑149 (MQ=255)                                                                                        ccGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGaa                                                                                                                               <  5:128251/82‑1 (MQ=255)                                                                                        ccGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGaa                                                                                                                               >  6:128251/1‑82 (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                                                  <  6:334825/130‑1 (MQ=255)                                                                                                                     cGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGgttgtt                                                  >  5:334825/1‑130 (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                   >  2:227688/1‑134 (MQ=255)                                                                                                                                                cgcAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAATATCGGa                   <  1:227688/134‑1 (MQ=255)                                                                                                                                                   aGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACTGATc  >  5:333573/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CGAAGCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACTGATC  >  NC_002947/3427910‑3428203

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