Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,577,812	T→C	5.8%	*512R (TGA→CGA)	PP_4906 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,577,812	0	T	C	5.8%	207.2 / 2.6	69	*512R (TGA→CGA)	PP_4906	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/38);  new base C (2/2);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

GCTCGCTGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCTGAGCAGCGC  >  NC_002947/5577666‑5577956                                                                                                                                                      |                                                                                                                                                 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                                                                                                     <  3:303875/149‑1 (MQ=255)                                            cGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTg                                                                                                     >  8:215857/1‑149 (MQ=255)                                            cGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTg                                                                                                     >  4:59440/1‑149 (MQ=255)                                             gACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTgg                                                                                                    <  6:295722/149‑1 (MQ=255)                                             gACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTgg                                                                                                    <  6:328002/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                                                                    <  5:166559/140‑1 (MQ=255)                                                                                      cTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGAt                                                                    >  6:166559/1‑140 (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)                                                                                                                                                      |                                                                                                                                                 GCTCGCTGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAG‑CGGGCCTTGTCGGTGTTTA‑CCAA‑CCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCTGAGCAGCGC  >  NC_002947/5577666‑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: ‑