Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,985,847:1	+C	9.2%	coding (89/753 nt)	flgA →	flagella basal body P‑ring formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,985,847	1	.	C	9.2%	237.3 / 21.0	65	coding (89/753 nt)	flgA	flagella basal body P‑ring formation protein
Reads supporting (aligned to +/- strand):  ref base . (31/28);  new base C (2/4);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

CCTTTCGTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCG  >  NC_002947/4985704‑4985994                                                                                                                                                 |                                                                                                                                                    ccTTTCGTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACgg                                                                                                                                                >  7:99646/1‑149 (MQ=255) ccTTTCGTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACgg                                                                                                                                                >  6:299916/1‑149 (MQ=255) ccTTTCGTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACgg                                                                                                                                                <  2:295057/149‑1 (MQ=255)  cTTTCGTTGCAACTCCGGTAGCGAGGCGGCTCGGGCTTTGCGTGTTTGAGGGTCATGCACACGAAAACGACATTTTCCCGCCGATTGACCCGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGt                                                                                                                                               <  5:64752/149‑1 (MQ=255)       gTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTGTTTTGAGTGTCCTGCACACGAAAACGACATTTTCCCGACGATTGACCCGCCTGCAGACTGGCCCCCTGGCCCCGCGGTGGCTGCTGATGACCGGCGTTCG‑CCAGGTGGCg                                                                                                                                           >  4:244101/1‑148 (MQ=255)          gCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAAcg                                                                                                                                       >  2:336122/1‑149 (MQ=255)           cAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCCACCGT‑GCGAAcgc                                                                                                                                      <  8:36084/149‑1 (MQ=255)            aaCTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGGTGATGCCCGGCGTTCGCCACCGT‑GCGAAcgcc                                                                                                                                     <  6:120697/149‑2 (MQ=255)            aaCTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAAcgcg                                                                                                                                     <  6:270225/149‑1 (MQ=255)             aCTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGt                                                                                                                                    >  4:239844/1‑149 (MQ=255)              cTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCCACCGT‑GCGAAcgccgg                                                                                                                                   <  6:331702/149‑4 (MQ=255)                 aGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCCACCGT‑GCGAAcgccgg                                                                                                                                   >  3:331379/1‑143 (MQ=255)                 aGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCCACCGT‑GCGAAcgccgg                                                                                                                                   <  4:331379/146‑4 (MQ=255)                 aGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCAACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCAc                                                                                                                                >  6:342781/1‑149 (MQ=255)                   gTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAAc                                                                                                                                        >  5:378849/1‑139 (MQ=255)                   gTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAAc                                                                                                                                        <  6:378849/139‑1 (MQ=255)                       cGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGcc                                                                                                                          <  7:40777/149‑1 (MQ=255)                             ggCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACa                                                                                                                    <  5:224812/149‑1 (MQ=255)                              gCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGt                                                                                                                                    >  6:215315/1‑132 (MQ=255)                                aCGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCt                                                                                                                 >  8:244980/1‑149 (MQ=255)                                      tttgctttTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCgg                                                                                                           >  7:115345/1‑149 (MQ=255)                                           ttttttGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGgtg                                                                                                         >  8:245439/1‑146 (MQ=255)                                           ttttttGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGgtg                                                                                                         <  7:245439/146‑1 (MQ=255)                                                      tCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGttt                                                                                           <  1:336122/149‑1 (MQ=255)                                                        aTGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCt                                                                                         <  8:216467/149‑1 (MQ=255)                                                         tGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCtt                                                                                        <  8:99646/149‑1 (MQ=255)                                                         tGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCtt                                                                                        <  3:239844/149‑1 (MQ=255)                                                            acacGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGaa                                                                                     >  1:251921/1‑149 (MQ=255)                                                            acacGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGaa                                                                                     >  8:225494/1‑149 (MQ=255)                                                            acacGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGaa                                                                                     <  7:244980/149‑1 (MQ=255)                                                              acGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAAtt                                                                                   >  5:248786/1‑149 (MQ=255)                                                                 aaaaCGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTctc                                                                                <  6:354423/149‑3 (MQ=255)                                                                    aCGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGt                                                                             >  8:62963/1‑149 (MQ=255)                                                                     cGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGtt                                                                            <  7:62963/149‑1 (MQ=255)                                                                       aCATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGa                                                                          <  5:299916/149‑1 (MQ=255)                                                                        cATTTTCCCGACGATTGTCACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGGCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGaa                                                                         <  3:244101/149‑1 (MQ=255)                                                                           tttCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGACCAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGAt                                                                      >  8:240992/1‑149 (MQ=255)                                                                              cccGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTa                                                                    <  6:32057/148‑1 (MQ=255)                                                                               ccGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATc                                                                  >  7:30743/1‑149 (MQ=255)                                                                               ccGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTGCTTGAATTCACCGTTGAAGATTATc                                                                  >  2:232233/1‑149 (MQ=255)                                                                                  acgaTTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGt                                                               <  5:342781/149‑1 (MQ=255)                                                                                           acGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCaccacc                                                       >  6:25523/1‑148 (MQ=255)                                                                                              ccTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGa                                                   <  7:142023/149‑1 (MQ=255)                                                                                                 gctgACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTg                                                >  3:197932/1‑149 (MQ=255)                                                                                                   tgACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTAt                                                                   >  1:223401/1‑128 (MQ=255)                                                                                                   tgACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTAt                                                                   <  2:223401/128‑1 (MQ=255)                                                                                                     aCTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTgccgg                                            >  2:7270/1‑149 (MQ=255)                                                                                                      cTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTgccggc                                           >  7:281762/1‑149 (MQ=255)                                                                                                       tGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAGCAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCAGCGTTGAAGATTATCTGTCCACCACCCAGACTgccggcc                                          >  6:81169/1‑149 (MQ=255)                                                                                                            cGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTAt                                     <  4:121467/149‑1 (MQ=255)                                                                                                                  ccACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATc                               >  6:171646/1‑149 (MQ=255)                                                                                                                   cACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGATCGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATc                               >  8:177071/1‑148 (MQ=255)                                                                                                                     cGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAg                            <  6:248786/149‑1 (MQ=255)                                                                                                                     cGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAg                            <  8:115345/149‑1 (MQ=255)                                                                                                                               gctgATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGc                  >  2:14036/1‑149 (MQ=255)                                                                                                                                ctgATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGc                  <  8:281762/148‑1 (MQ=255)                                                                                                                                 tgATGCCCGGCGTTCA‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTAt                                                                   <  1:325040/98‑1 (MQ=255)                                                                                                                                 tgATGCCCGGCGTTCA‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTAt                                                                   >  2:325040/1‑98 (MQ=255)                                                                                                                                    tGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGa              <  7:225494/148‑1 (MQ=255)                                                                                                                                    tGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGAc             >  8:165414/1‑149 (MQ=255)                                                                                                                                       ccGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCa          >  8:226308/1‑149 (MQ=255)                                                                                                                                        cGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGTCCCAc         <  1:232233/149‑1 (MQ=255)                                                                                                                                        cGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCAc         <  5:25523/149‑1 (MQ=255)                                                                                                                                         ggCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACg        >  5:188831/1‑149 (MQ=255)                                                                                                                                          gCGTTCGCCACCGT‑GCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCCCCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGc       >  7:341909/1‑149 (MQ=255)                                                                                                                                                g‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTgcg  >  5:71588/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CCTTTCGTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCG  >  NC_002947/4985704‑4985994

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