Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,985,854	Δ1 bp	9.7%	coding (96/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,854	0	G	.	9.7%	189.0 / 14.7	62	coding (96/753 nt)	flgA	flagella basal body P‑ring formation protein
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base . (2/4);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.15e-01

GTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCG  >  NC_002947/4985710‑4985994                                                                                                                                                  |                                                                                                                                             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                                                                                                                                   <  4:331379/146‑4 (MQ=255)           aGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCCACCGT‑GCGAAcgccgg                                                                                                                                   >  3:331379/1‑143 (MQ=255)           aGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCAACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCAc                                                                                                                                >  6:342781/1‑149 (MQ=255)             gTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAAc                                                                                                                                        <  6:378849/139‑1 (MQ=255)             gTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAAc                                                                                                                                        >  5:378849/1‑139 (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                                                                                                         <  7:245439/146‑1 (MQ=255)                                     ttttttGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGgtg                                                                                                         >  8:245439/1‑146 (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                                                                                     >  8:225494/1‑149 (MQ=255)                                                      acacGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGaa                                                                                     <  7:244980/149‑1 (MQ=255)                                                      acacGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGaa                                                                                     >  1:251921/1‑149 (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)                                                                                                                                                  |                                                                                                                                             GTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCG‑CACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCG  >  NC_002947/4985710‑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: ‑