Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,423,807	A→T	6.0%	T14S (ACA→TCA) ‡	PP_3915 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,423,807	0	A	T	6.0%	195.5 / 2.2	67	T14S (ACA→TCA) ‡	PP_3915	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base T (2/2);  total (31/36)
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.26e-01

GCGGCCCAGACCAGGCCGGCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATGGGCA  >  NC_002947/4423659‑4423951                                                                                                                                                     |                                                                                                                                                 gcgGCCCAGACCAGGCCGGCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGca                                                                                                                                                  >  8:383584/1‑149 (MQ=255)      ccAGACCAGGCCGGCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACagag                                                                                                                                             <  3:416568/149‑1 (MQ=255)           gctggCCGGCAAGGAGCGCCAGCTCTCACAGCACGAAACGCATACCTAACCTACACCCGCTGCCTCGTCGGCGGAATTGATTGATGAGAGGTATCAGCTGTGTGTCGGAATAAAAGAAACTTGACCATTGTTTCAAGCACTGAGGTgga                                                                                                                                        <  4:156816/146‑1 (MQ=255)           ccAGGCCGGCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTgga                                                                                                                                        <  2:409424/149‑1 (MQ=255)            cAGGCCGGCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTggag                                                                                                                                       <  4:16847/149‑1 (MQ=255)             aGGCCGGCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTggag                                                                                                                                       >  7:353658/1‑148 (MQ=255)                 cGGCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCTCTGTGCTggaaacaat                                                                                                                                  <  1:67271/149‑7 (MQ=37)                 cGGCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAAtt                                                                                                                                  >  4:119759/1‑149 (MQ=255)                   gCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTgc                                                                                                                                <  4:27643/149‑1 (MQ=255)                   gCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCAGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTgc                                                                                                                                <  6:10060/149‑1 (MQ=255)                    cAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCTCTGTGCTggaa                                                                                                                                       >  7:205836/1‑140 (MQ=255)                    cAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCTCTGTGCTggaa                                                                                                                                       <  8:205836/141‑2 (MQ=255)                     aaCGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTgcgc                                                                                                                              <  3:259323/149‑1 (MQ=255)                        gAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGcgccgc                                                                                                                           <  8:215715/149‑1 (MQ=255)                             ccAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCg                                                                                                                      >  1:351139/1‑149 (MQ=255)                              cAGCTCTCACCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCGAGTTCAACATTGTTTCCAGCACAGAGCTGGCGTAATAGCGCCGCGCCCGt                                                                                                                     >  4:51540/1‑149 (MQ=37)                                gCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCAAGGATTGATGAGAGCTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTTTTTCCAGCACAGAGCTCTAGGTATTGAGCCGCGCCCGt                                                                                                                     >  5:92799/1‑147 (MQ=38)                                 ctctCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAg                                                                                                                  >  2:204520/1‑149 (MQ=255)                                  gctcCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGc                                                                                                                 <  6:265326/148‑1 (MQ=255)                                   caccccGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGc                                                                                                                 >  8:403011/3‑148 (MQ=255)                                     ccccGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGAc                                                                                                              >  4:90635/1‑149 (MQ=255)                                      cccGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGAc                                                                                                              <  5:143956/148‑1 (MQ=255)                                       ccGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACt                                                                                                             >  1:74218/1‑148 (MQ=255)                                          cACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAAc                                                                                                         <  5:391241/149‑1 (MQ=255)                                             gAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGcc                                                                                                      >  1:438201/1‑149 (MQ=255)                                             gAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGcc                                                                                                      <  1:204520/149‑1 (MQ=255)                                                  cgaacccGATACTACAAGCGCTGCCTAGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACGATTGTTTCCAGCAGAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGa                                                                                                 <  6:92799/147‑1 (MQ=255)                                                    aacccGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGt                                                                                               >  3:284285/1‑149 (MQ=255)                                                       ccGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCAt                                                                                            >  2:427893/1‑149 (MQ=255)                                                          aaCCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCAtgct                                                                                         >  6:123343/1‑149 (MQ=255)                                                          aaCCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCAtgct                                                                                         <  4:46441/149‑1 (MQ=255)                                                           agcTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGTTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCAtgctg                                                                                        <  2:416811/147‑1 (MQ=255)                                                             cTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGCCCGTGCGCGACTGAACGCCCTGGCGTCATgcggct                                                                                      >  7:203707/1‑149 (MQ=255)                                                                cacCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGctgctg                                                                                     <  1:427893/147‑1 (MQ=255)                                                                cacCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTggg                                                                                   <  8:353658/149‑1 (MQ=255)                                                                  cccGCTGCCTCGTCGGCGGCATGGAGTGATGAGAGGTTTCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGAt                                                                                 <  1:366828/149‑1 (MQ=255)                                                                   ccGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATg                                                                                >  4:242595/1‑149 (MQ=255)                                                                    cGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGt                                                                               <  3:119759/149‑1 (MQ=255)                                                                      cTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTc                                                                             <  3:61698/149‑1 (MQ=255)                                                                      cTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTc                                                                             >  3:387440/1‑149 (MQ=255)                                                                          cTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCCCCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCGGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGCTGTTCGgtt                                                                         >  6:290939/1‑149 (MQ=255)                                                                             gCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTc                                                                        <  3:325182/147‑1 (MQ=255)                                                                             gCCGGCGGAAGGGATTTATGTGTGGTATCAGCTGTGAGTCGTATTAAAAGCCACTTCACCTTTGTTTCCAGCACAGAGCTGGTGGATTTGCGCCGCGACCGTGAGCTACTTAACGCCCTGGTGTCATGCTGCTGGGTTGTTCGTTTCga                                                                      <  8:176245/149‑1 (MQ=255)                                                                                     cATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCAtt                                                              >  8:356702/1‑149 (MQ=255)                                                                                      aTGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTc                                                             >  8:7505/1‑149 (MQ=255)                                                                                              atgaGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGgga                                                     <  1:7451/149‑1 (MQ=255)                                                                                               tgaGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGggat                                                    <  2:74218/149‑1 (MQ=255)                                                                                               tgaGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGgga                                                     <  7:403011/148‑1 (MQ=255)                                                                                                 agagGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGggatg                                                   <  7:383584/148‑1 (MQ=255)                                                                                                     gTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGatggatgg                                              <  4:387440/149‑1 (MQ=255)                                                                                                           gCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTa                                        <  5:123343/149‑1 (MQ=255)                                                                                                                 aGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGGACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGa                                   <  2:438201/148‑1 (MQ=255)                                                                                                                  gTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATc                                 <  2:351139/149‑1 (MQ=255)                                                                                                                   tCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAAt                                                         <  7:259459/124‑1 (MQ=255)                                                                                                                   tCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAAt                                                         >  8:259459/1‑124 (MQ=255)                                                                                                                   tCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCg                                >  6:396396/1‑149 (MQ=255)                                                                                                                   tCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGACCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTAGAATGGGATGGATGGCGACTACACGATCg                                >  6:176838/1‑149 (MQ=255)                                                                                                                          aaaaGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCAtt                           >  8:335431/1‑147 (MQ=255)                                                                                                                           aaaGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCAtt                           >  8:398564/1‑146 (MQ=255)                                                                                                                           aaaGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCAtt                           >  8:158440/1‑146 (MQ=255)                                                                                                                           aaaGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTgag                        <  7:335431/149‑1 (MQ=255)                                                                                                                                         aTTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGATTGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTa          <  1:319072/149‑1 (MQ=255)                                                                                                                                             tttCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATg      >  1:411272/1‑149 (MQ=255)                                                                                                                                             tttCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATg      >  1:394913/1‑149 (MQ=255)                                                                                                                                              ctCCAGCTCTGTGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATgg     >  7:228883/2‑149 (MQ=255)                                                                                                                                               tCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATgg     <  5:338276/148‑1 (MQ=255)                                                                                                                                                 cAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATGGGc   >  7:239121/1‑148 (MQ=255)                                                                                                                                                  aGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATGGGCa  >  7:46964/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GCGGCCCAGACCAGGCCGGCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATGGGCA  >  NC_002947/4423659‑4423951

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