Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,423,809	A→T	7.0%	T14T (ACA→ACT) ‡	PP_3915 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,423,809	0	A	T	7.0%	185.6 / 2.5	66	T14T (ACA→ACT) ‡	PP_3915	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base T (2/3);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCAGACCAGGCCGGCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATGGGCA  >  NC_002947/4423664‑4423951                                                                                                                                                  |                                                                                                                                               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:387440/1‑149 (MQ=255)                                                                 cTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTc                                                                             <  3:61698/149‑1 (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                                                         >  8:259459/1‑124 (MQ=255)                                                                                                              tCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAAt                                                         <  7:259459/124‑1 (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:158440/1‑146 (MQ=255)                                                                                                                      aaaGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCAtt                           >  8:398564/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)                                                                                                                                                  |                                                                                                                                               CCAGACCAGGCCGGCAACGAGCGCCAGCTCTCCCCGCACGAAACCCAACCCGAACCTACACCCGCTGCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATGGGCA  >  NC_002947/4423664‑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: ‑