Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,970,525	A→G	6.2%	W257R (TGG→CGG)	PP_5236 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,970,525	0	A	G	6.2%	236.2 / 4.4	82	W257R (TGG→CGG)	PP_5236	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/36);  new base G (3/2);  total (44/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

CGACGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACTGTGG  >  NC_002947/5970381‑5970669                                                                                                                                                  |                                                                                                                                                 cGACGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCgg                                                                                                                                              <  6:32953/149‑1 (MQ=255) cGACGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCgg                                                                                                                                              <  2:257176/149‑1 (MQ=255)  gACGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCggg                                                                                                                                             >  3:234842/1‑149 (MQ=255)   aCGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCgggg                                                                                                                                            <  8:195103/149‑1 (MQ=255)   aCGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATACAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCgggg                                                                                                                                            >  5:70291/1‑149 (MQ=255)    cGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCgggg                                                                                                                                            >  3:287621/1‑148 (MQ=255)     ggCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCgggg                                                                                                                                            >  5:266074/1‑147 (MQ=255)     ggCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGgca                                                                                                                                          <  5:294210/149‑1 (MQ=255)           cACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCTGCCCCGCTGGGGcggcag                                                                                                                                      <  8:303126/147‑1 (MQ=255)             ccTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCAATGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCCGCGGGGCAGCAGGTac                                                                                                                                  >  1:277020/1‑149 (MQ=255)                 gTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCTGCCCCGCTGGGGcggcagccttgat                                                                                                                               <  2:201320/148‑8 (MQ=255)                 gTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTacac                                                                                                                                <  1:319927/147‑1 (MQ=255)                 gTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTAcaca                                                                                                                               <  7:76276/148‑1 (MQ=255)                 gTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTAcaca                                                                                                                               <  4:287621/148‑1 (MQ=255)                    aGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGc                                                                                                                             >  5:278390/1‑147 (MQ=255)                      ccTGTGGCAGCTTCCGCCGGCCTGCGCCCGCCATCCAGCGCAGAGCATTACGCAACCGCTCGTCCCCGGTTGCGCACCGCGGGCCATTGCCCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAgg                                                                                                                         <  6:343519/149‑1 (MQ=255)                      ccTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAgg                                                                                                                         <  5:161821/149‑1 (MQ=255)                       cTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAggg                                                                                                                        >  6:209854/1‑149 (MQ=255)                                 ttCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGc                                                                                                              >  2:264762/1‑149 (MQ=255)                                   cAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGc                                                                                                              <  6:278390/147‑1 (MQ=255)                                    aGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCtt                                                                                                            <  7:103186/148‑1 (MQ=255)                                       cGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCAGGCCATAGACCAACAGCTCGAACACAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCCGGTACACAGCAAGCGCCAGGCTCGCATCGaa                                                                                                        >  1:272518/1‑149 (MQ=255)                                               gACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTgg                                                                                                                   >  2:47936/1‑130 (MQ=255)                                               gACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTgg                                                                                                                   <  1:47936/130‑1 (MQ=255)                                                aCAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCggg                                                                                               >  7:306456/1‑149 (MQ=255)                                                  aGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCgggg                                                                                              >  1:169/1‑148 (MQ=255)                                                  aGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCgggg                                                                                              <  2:169/148‑1 (MQ=255)                                                    ccATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCTAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGgggcgg                                                                                           >  6:265788/1‑149 (MQ=255)                                                    ccATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGgggcgg                                                                                           <  6:292020/149‑1 (MQ=255)                                                     cATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGgggcggg                                                                                          <  4:234842/149‑1 (MQ=255)                                                       tCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAagc                                                                                                      <  8:292489/135‑1 (MQ=255)                                                       tCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAagc                                                                                                      >  7:292489/1‑135 (MQ=255)                                                       tCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGgggcgggc                                                                                         >  8:30001/1‑148 (MQ=255)                                                           gctcagcgcaTTACGCAAACGCTCGTCCGAGGTTGCGCCCAGCGGGCGATTGACAAACAGCTCGAAAGCAT‑CAAGGCTGCCGCCCCAGCGGTGCAGGAGGTTCACAGCAAGGGCCTGTATAGCTTCGAATCTGCGGGGGCGGGCGATCa                                                                                    <  2:157160/146‑1 (MQ=255)                                                                 cgcaTTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCAt                                                                              <  3:101835/149‑1 (MQ=255)                                                                     ttACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATgg                                                                            <  2:1649/147‑1 (MQ=255)                                                                     ttACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATgg                                                                            >  1:1649/1‑147 (MQ=255)                                                                        cGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCa                                                                        <  6:182255/148‑1 (MQ=255)                                                                        cGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAg                                                                       >  2:3789/1‑149 (MQ=255)                                                                        cGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAg                                                                       >  6:42967/1‑149 (MQ=255)                                                                               gCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTaa                                                                >  6:262948/1‑149 (MQ=255)                                                                                    tCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGc                                                            >  5:213231/1‑148 (MQ=255)                                                                                       cAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAgg                                                        >  8:64003/1‑149 (MQ=255)                                                                                       cAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAgg                                                        <  1:264762/149‑1 (MQ=255)                                                                                                cacaGCGGGCCATTGACCAACAGCTCGAACGCATGCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGt                                                >  1:218389/1‑149 (MQ=255)                                                                                                       ggCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATgg                                        >  1:175173/1‑149 (MQ=255)                                                                                                        gCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGt                                       <  5:209854/149‑1 (MQ=255)                                                                                                            ttGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGAc                                   >  1:262033/1‑149 (MQ=255)                                                                                                              gACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGAccc                                 >  1:300786/1‑149 (MQ=255)                                                                                                              gACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGAcc                                  <  8:306456/148‑1 (MQ=255)                                                                                                              gACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGAcc                                  <  2:262033/148‑1 (MQ=255)                                                                                                               aCCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCg                                >  1:3524/1‑149 (MQ=255)                                                                                                                  aaCAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCtt                             <  5:311536/149‑1 (MQ=255)                                                                                                                   aCAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTg                            >  8:293862/1‑149 (MQ=255)                                                                                                                   aCAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTg                            >  3:129596/1‑149 (MQ=255)                                                                                                                       cTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCa                        <  1:3789/149‑1 (MQ=255)                                                                                                                          gAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACATCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGTATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGAt                     <  6:70291/149‑1 (MQ=255)                                                                                                                          gAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATg                                         <  3:41687/129‑1 (MQ=255)                                                                                                                          gAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATg                                         >  4:41687/1‑129 (MQ=255)                                                                                                                          gAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGa                      >  4:285047/1‑148 (MQ=255)                                                                                                                               cAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAgcgc                >  4:135355/1‑149 (MQ=255)                                                                                                                                aT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGcgcg               >  4:113902/1‑149 (MQ=255)                                                                                                                                aT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGcgcg               <  7:30001/149‑1 (MQ=255)                                                                                                                                aT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTGTTCATGGTCGACCCGCTTGCCCCGATAGcgcg               >  6:235073/1‑149 (MQ=255)                                                                                                                                    aaGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCa            <  5:262948/149‑1 (MQ=255)                                                                                                                                      ggCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCgg              <  7:134147/145‑1 (MQ=255)                                                                                                                                      ggCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCgg              >  8:134147/1‑145 (MQ=255)                                                                                                                                      ggCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGc          <  8:238152/149‑1 (MQ=255)                                                                                                                                       gCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGCCCCGCTTGCCCAGATAGCGCGGCAGCg         >  5:26731/1‑149 (MQ=255)                                                                                                                                        ctgctgCCCCGCTGGGGCGGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGa        >  3:310289/1‑149 (MQ=255)                                                                                                                                        ctgctgCCCCGCTGGGGCGGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCACATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGa        >  8:20028/1‑149 (MQ=255)                                                                                                                                        ctgctgCCCCGCGGGGGCGGCAGGTACACAGCAAGGGCCTGGATCGCTTCGCAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGGTGTCCTGGTCGCCCCGCTGGCCCAGATAGCGCGGCAGCgc        >  5:72/1‑148 (MQ=255)                                                                                                                                         tGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGAc       >  3:279860/1‑149 (MQ=255)                                                                                                                                          gccgccCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACt      >  4:195764/1‑149 (MQ=255)                                                                                                                                           ccgccCCAGCGGGGCAGCAGGGACACCGCAAGGGCCTGGATCGCGGCGACGCCGCCGGGGCGGGCGCTCAGCGCATGGCGCCGCCAGGAAGAGCGAGGGTCCGCGGGGTCCTGGGCGCCCCGCTGGCCCAGAGAGCGCGGCAGCgccgg     >  5:200492/1‑145 (MQ=255)                                                                                                                                             gccCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACtgtg   <  6:266074/149‑1 (MQ=255)                                                                                                                                              ccccAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACTGTgg  <  3:195764/149‑1 (MQ=255)                                                                                                                                                ccAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCa                        >  6:185698/1‑125 (MQ=255)                                                                                                                                                ccAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCa                        <  5:185698/125‑1 (MQ=255)                                                                                                                                                ccAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACTGTgg  >  5:155376/1‑147 (MQ=255)                                                                                                                                                 cAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACtgt    <  5:97223/144‑1 (MQ=255)                                                                                                                                                 cAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACtgt    >  6:97223/1‑144 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CGACGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAT‑CAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACTGTGG  >  NC_002947/5970381‑5970669

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑