Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,457,146	G→A	6.2%	Y111Y (TAC→TAT)	PP_1275 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,457,146	0	G	A	6.2%	189.2 / 2.4	64	Y111Y (TAC→TAT)	PP_1275	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/28);  new base A (2/2);  total (34/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 = 8.12e-01

ATGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGT  >  NC_002947/1457000‑1457288                                                                                                                                                    |                                                                                                                                               aTGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTa                                                                                                                                              >  8:31031/1‑149 (MQ=255)  tGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTAt                                                                                                                                             >  3:36641/1‑149 (MQ=255)  tGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTAt                                                                                                                                             >  3:110378/1‑149 (MQ=255)  tGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTAt                                                                                                                                             >  2:166018/1‑149 (MQ=255)    gaggaAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATcg                                                                                                                                           >  1:287453/1‑149 (MQ=255)     aggaAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATcgc                                                                                                                                          <  1:166018/149‑1 (MQ=255)      ggaAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCGATTCACCTGGCCCG‑CGGTATcgcg                                                                                                                                         >  7:283217/1‑149 (MQ=255)        taaCACATCAATTGGCAAACGGACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCgg                                                                                                                                        <  2:248179/147‑1 (MQ=255)            acaTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATcgcg                                                                                                                                         >  8:277693/1‑143 (MQ=255)            acaTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATcgcg                                                                                                                                         <  7:277693/143‑1 (MQ=255)            acaTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGt                                                                                                                                   >  5:170799/1‑149 (MQ=255)            acaTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGt                                                                                                                                   >  5:364990/1‑149 (MQ=255)             caTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCGATTCACCTGGCCCG‑CGGTATCGCGGGGCGtt                                                                                                                                  <  8:283217/149‑1 (MQ=255)             caTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGtt                                                                                                                                  >  8:202667/1‑149 (MQ=255)                 aaTTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAg                                                                                                                              <  4:168945/149‑1 (MQ=255)                  aTTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGATGTAgg                                                                                                                             >  4:215457/1‑149 (MQ=255)                    tGGCAAACGGACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCGGCAATTCACCTGGCCCG‑CGATACCGCgggccag                                                                                                                                   <  3:282053/141‑5 (MQ=255)                      gCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGtt                                                                                                                         <  6:364990/149‑1 (MQ=255)                      gCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGtt                                                                                                                         >  4:50371/1‑149 (MQ=255)                      gCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCATTGTAGGTGtt                                                                                                                         <  7:31031/149‑1 (MQ=255)                          aCGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCg                                                                                                                     <  3:6801/149‑1 (MQ=255)                            gTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATc                                                                                                                                            >  7:340568/1‑124 (MQ=255)                            gTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATc                                                                                                                                            <  8:340568/124‑1 (MQ=255)                               ccTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCa                                                                                                                 <  4:36641/148‑1 (MQ=255)                               ccTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAg                                                                                                                >  3:340438/1‑149 (MQ=255)                                cTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAgg                                                                                                               >  4:206361/1‑149 (MQ=255)                                 ttGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGt                                                                                                              >  8:134988/1‑149 (MQ=255)                                    cAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCa                                                                                                           >  3:123796/1‑149 (MQ=255)                                         gctcggctcgGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTa                                                                                                      >  2:93999/1‑149 (MQ=255)                                             ggctcggAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGc                                                                                                  >  1:67210/1‑149 (MQ=255)                                                tcggAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGTCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CCGTATCGCGGGGCGTTGTAGGAGTTCTCGGCCAGGTTCACGGTATCGCCCg                                                                                               >  2:127642/1‑149 (MQ=255)                                                tcggAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGTAGTTCTAGGCCAGGTTCACGGTATCGCCCg                                                                                               >  8:345115/1‑149 (MQ=255)                                                   gAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGt                                                                                             <  1:93999/149‑1 (MQ=255)                                                    aaTTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGt                                                                                             <  3:27435/147‑1 (MQ=255)                                                    aaTTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGt                                                                                             >  4:27435/1‑147 (MQ=255)                                                     aTTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCg                                                                                          <  4:110378/149‑1 (MQ=255)                                                         cAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGc                                                                                                  <  8:111674/137‑1 (MQ=255)                                                         cAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGc                                                                                                  >  7:111674/1‑137 (MQ=255)                                                          aGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTAcgc                                                                                     <  4:340438/149‑1 (MQ=255)                                                                   ggCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGTAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCAcgcg                                                                            >  2:118881/1‑149 (MQ=255)                                                                      cggcccccGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACgcgc                                                                           <  8:169925/141‑1 (MQ=255)                                                                            ccGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAg                                                                   <  6:170799/149‑1 (MQ=255)                                                                                aaCGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCg                                                               <  3:206361/149‑1 (MQ=255)                                                                                 aCGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCgg                                                              >  4:54377/1‑149 (MQ=255)                                                                                  cGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGt                                                             <  7:134988/149‑1 (MQ=255)                                                                                     aGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAg                                                          >  5:220305/1‑149 (MQ=255)                                                                                      gcggcgGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGc                                                         <  3:54377/149‑1 (MQ=255)                                                                                              aaTGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGt                                                                                             >  4:279714/1‑105 (MQ=255)                                                                                              aaTGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGt                                                                                             <  3:279714/105‑1 (MQ=255)                                                                                                         gTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGttct                                      <  3:50371/149‑1 (MQ=255)                                                                                                            ttGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCg                                          >  3:355946/1‑142 (MQ=255)                                                                                                            ttGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCg                                          <  4:355946/142‑1 (MQ=255)                                                                                                              gCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCAt                                  <  7:345115/148‑1 (MQ=255)                                                                                                              gCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCAt                                  <  2:287453/148‑1 (MQ=255)                                                                                                                cACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCg                               >  3:159276/1‑149 (MQ=255)                                                                                                                 agtCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCa                                                                    <  8:229396/109‑1 (MQ=255)                                                                                                                          ccGGCTTCAATTCACCTGGCCCG‑CGGTAGCGCGGGGCGGTGTCGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGGCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAg                     <  3:215457/149‑1 (MQ=255)                                                                                                                               tGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGAGTACCTGCAAGGCACg                >  8:16264/1‑149 (MQ=255)                                                                                                                                gCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACgg               >  5:107797/1‑149 (MQ=255)                                                                                                                                  aaTTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTc             <  4:123796/149‑1 (MQ=255)                                                                                                                                         acgccccG‑CGATACCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTcggcgg        <  5:41705/143‑1 (MQ=255)                                                                                                                                         acgccccG‑CGATACCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGATGCCGTTCTTCATGCGTACCTGCAAGGCACGGTcggcgg        >  6:41705/5‑147 (MQ=255)                                                                                                                                           gccccG‑CGATACCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGa      >  3:309880/3‑147 (MQ=255)                                                                                                                                             cccG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAgt  >  5:169172/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                               ATGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGT  >  NC_002947/1457000‑1457288

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