Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,577,810	G→C	5.8%	G511A (GGT→GCT)	PP_4906 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,577,810	0	G	C	5.8%	227.7 / 3.4	69	G511A (GGT→GCT)	PP_4906	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/37);  new base C (2/2);  total (30/39)
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.16e-01

CTTCGCTCGCTGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCTGAGCAGCGC  >  NC_002947/5577662‑5577956                                                                                                                                                        |                                                                                                                                                   ctTCGCTCGCTGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAggg                                                                                                                                                    >  6:141877/1‑149 (MQ=255) ctTCGCTCGCTGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAggg                                                                                                                                                    >  4:280847/1‑149 (MQ=255) ctTCGCTCGCTGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAggg                                                                                                                                                    <  1:83612/149‑1 (MQ=255)  ttCGCTCGCTGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGt                                                                                                                                                   >  3:217385/1‑149 (MQ=255)     gctcgctGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGTCCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGa                                                                                                                                                <  5:174178/149‑1 (MQ=255)       tcgctGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGc                                                                                                                                              >  1:298769/1‑149 (MQ=255)          ctGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCtt                                                                                                                                           <  8:69343/149‑1 (MQ=255)             gCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTg                                                                                                                                          <  2:298769/147‑1 (MQ=255)              cAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGCTCAACCCTTgcgg                                                                                                                                       <  5:207103/149‑2 (MQ=255)                 ccAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTa                                                                                                                                     >  1:218820/1‑148 (MQ=255)                 ccAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTa                                                                                                                                     >  3:80666/1‑148 (MQ=255)                  cAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTAc                                                                                                                                    >  6:79847/1‑148 (MQ=255)                  cAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTAc                                                                                                                                    >  3:34929/1‑148 (MQ=255)                    aCGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAATCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAAt                                                                                                                                 <  2:304273/149‑1 (MQ=255)                         cAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGccc                                                                                                                            <  6:148465/149‑1 (MQ=255)                          aGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCCGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGccc                                                                                                                            <  1:88883/149‑1 (MQ=255)                            gTGCTGATGCTGCTGGTGGCCGACCATGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAACCCGCAAGGGTTGAGCCTTGCGCTACAATCGc                                                                                                                              <  1:196634/145‑1 (MQ=255)                            gTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGc                                                                                                                              >  2:196634/1‑144 (MQ=255)                               cTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCttt                                                                                                                      <  7:188139/149‑1 (MQ=255)                                  aTGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTc                                                                                                                     <  1:252480/147‑1 (MQ=255)                                  aTGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTc                                                                                                                     <  6:311843/147‑1 (MQ=255)                                           gtggCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTACCCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGga                                                                                                            <  1:246174/148‑1 (MQ=255)                                           gtggCCGACCAGGGCTGCCGGCCTTTGGACGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGTTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCTCCCCTCTTTCCAAACTGgaga                                                                                                          >  6:78947/1‑149 (MQ=255)                                            tggCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGAc                                                                                                         >  5:328002/1‑149 (MQ=255)                                              gCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGCTCAACCCTTgcggttggt                                                                                                                                  <  2:274503/122‑7 (MQ=255)                                              gCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGCTCAACCCTTgcggttggt                                                                                                                                  >  1:274503/1‑116 (MQ=255)                                                cGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTg                                                                                                     >  8:215857/1‑149 (MQ=255)                                                cGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTg                                                                                                     <  3:303875/149‑1 (MQ=255)                                                cGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTg                                                                                                     >  4:59440/1‑149 (MQ=255)                                                 gACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTgg                                                                                                    <  6:328002/149‑1 (MQ=255)                                                 gACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTgg                                                                                                    <  6:295722/149‑1 (MQ=255)                                                 gACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTgg                                                                                                    <  4:70982/149‑1 (MQ=255)                                                      gggCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGAc                                                                                               >  4:265164/1‑149 (MQ=255)                                                       ggCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGAcc                                                                                              <  8:73232/149‑1 (MQ=255)                                                       ggCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGAcc                                                                                              <  4:34929/149‑1 (MQ=255)                                                          tGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGAc                                                                                           <  4:217385/149‑1 (MQ=255)                                                            ccGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGAc                                                                                           <  2:218820/147‑1 (MQ=255)                                                              ggCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACtttt                                                                                       >  3:267989/1‑149 (MQ=255)                                                               gCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTc                                                                                      <  5:79847/149‑1 (MQ=255)                                                                ccTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTcc                                                                                     <  5:141877/149‑1 (MQ=255)                                                                ccTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTcc                                                                                     <  3:280847/149‑1 (MQ=255)                                                                       ccGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGAATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCAGACTCGAGACCGTGGATGACCGACTATTCCGGCCTGc                                                                              >  8:251234/1‑149 (MQ=255)                                                                             tCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTgg                                                                        <  7:105667/149‑1 (MQ=255)                                                                             tCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTgg                                                                        <  3:59440/149‑1 (MQ=255)                                                                               aGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTg                                                                      <  5:1752/149‑1 (MQ=255)                                                                                      aaGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAg                                                               >  1:172273/1‑149 (MQ=255)                                                                                         ccTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTg                                                            <  4:267989/149‑1 (MQ=255)                                                                                          cTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGAt                                                                    >  6:166559/1‑140 (MQ=255)                                                                                          cTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGAt                                                                    <  5:166559/140‑1 (MQ=255)                                                                                                  aGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACctgc                                                   <  2:172273/149‑1 (MQ=255)                                                                                                   gACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACctgct                                                  >  7:16005/1‑149 (MQ=255)                                                                                                     ccGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCtgctgc                                                <  7:256091/149‑1 (MQ=255)                                                                                                      cgcgCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGc                                                                                  >  7:266074/1‑114 (MQ=255)                                                                                                      cgcgCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGc                                                                                  <  8:266074/114‑1 (MQ=255)                                                                                                            aTGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGCCCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCt                                         >  7:206283/1‑149 (MQ=255)                                                                                                                aTCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGAt                                     >  5:164546/1‑149 (MQ=255)                                                                                                                    aG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTcgc                                 >  5:36451/1‑149 (MQ=255)                                                                                                                         ggCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTcaccgc                              >  8:244567/1‑148 (MQ=255)                                                                                                                                tCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCACCGCAACTGATc                      <  7:244567/149‑1 (MQ=255)                                                                                                                                  ggTGTTTA‑CCAA‑CCGCAATGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATcc                     <  4:80666/148‑1 (MQ=255)                                                                                                                                   gtgtTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTg                   >  5:192299/1‑149 (MQ=255)                                                                                                                                      tttA‑CCAA‑CCGGAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGAGCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTc                <  7:251234/149‑1 (MQ=255)                                                                                                                                      tttA‑CCAA‑CCGCAAGGGTTGAGCCTTTCGCTACAATCGCCCGTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGt                 <  5:78947/148‑1 (MQ=255)                                                                                                                                      tttA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGt                 <  6:192299/148‑1 (MQ=255)                                                                                                                                        tA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGa              >  7:174829/1‑149 (MQ=255)                                                                                                                                         a‑CCAA‑CCGCAAGGGTTGCGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGAc             >  8:236142/1‑149 (MQ=255)                                                                                                                                         a‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGAt                                                                    >  4:43668/1‑94 (MQ=255)                                                                                                                                         a‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGAt                                                                    <  3:43668/94‑1 (MQ=255)                                                                                                                                            cAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTCGAGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGACGACCt           >  3:108334/1‑149 (MQ=255)                                                                                                                                              tagcGCAAGGCTCAACCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCt           >  5:24475/4‑148 (MQ=255)                                                                                                                                                   cAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTcgc                                 <  6:10068/121‑1 (MQ=255)                                                                                                                                                   cAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTcgc                                 >  5:10068/1‑121 (MQ=255)                                                                                                                                                      gggTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCTGAGCAgcgc  <  7:215857/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                   CTTCGCTCGCTGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCTGAGCAGCGC  >  NC_002947/5577662‑5577956

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