Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,909,632	T→C	6.4%	L527L (CTT→CTC)	PP_3450 →	TPR repeat‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,909,632	0	T	C	6.4%	178.4 / 2.3	62	L527L (CTT→CTC)	PP_3450	TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base C (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

AACTGGCCGGCGACAGTGCGCAGAGCACCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCTGGTCGACGGTGATGA  >  NC_002947/3909484‑3909778                                                                                                                                                     |                                                                                                                                                   aaCTGGCCGGCGACAGTGCGCAGAGCACCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCtt                                                                                                                                                    <  4:255092/149‑1 (MQ=255) aaCTGGCCGGCGACAGTGCGCAGAGCACCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCtt                                                                                                                                                    <  4:250697/149‑1 (MQ=255) aaCTGGCCGGCGACAGTGCGCAGAGCACCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCtt                                                                                                                                                    >  3:89150/1‑149 (MQ=255) aaCTGGCCGGCGACAGTGCGCAGAGCACCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCtt                                                                                                                                                    >  2:229610/1‑149 (MQ=255)    tGGCCGGCGACAGTGCGCAGAGCACCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAACGCTCAGCGCATGCAATCGATCCCCTACCTTGaa                                                                                                                                                 >  6:57504/1‑149 (MQ=255)        cGGCGACAGTGCGCAGAGCACCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTCGAACACg                                                                                                                                             >  3:152193/1‑149 (MQ=255)        cGGCGACAGTGCGCAGAGCACCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTCGAACACg                                                                                                                                             >  3:147333/1‑149 (MQ=255)                gTGCGCAGAGCACCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCAc                                                                                                                                     <  1:169892/149‑1 (MQ=255)                     cAGAGCACCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACtt                                                                                                                                >  4:211579/1‑149 (MQ=255)                      agagCACCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTc                                                                                                                               >  1:320862/1‑149 (MQ=255)                           aCCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTcccc                                                                                                                            >  1:162904/1‑147 (MQ=255)                            cccccTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGCCACGACAGACAAGCCTGAGCGCATGCAATCGATCCCCTACCTTGAACACCCCACTCACGACTTCCCCCAg                                                                                                                         >  1:266107/4‑149 (MQ=255)                            ccGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAg                                                                                                                         <  8:126950/149‑1 (MQ=255)                                        ccGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCtt                                                                                                             <  7:210888/149‑1 (MQ=255)                                         cGAAGCCGTGCGCCTGGCCCCCGACAAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCAACTCACGACTTCCCCGAGGACTACCGCCTTg                                                                                                            >  7:125142/1‑149 (MQ=255)                                            aGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCg                                                                                                         <  8:88847/149‑1 (MQ=255)                                               cGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACGAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGaaa                                                                                                      <  8:191793/149‑1 (MQ=255)                                                 tGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCCTGGCGAAaca                                                                                                    >  2:52785/1‑149 (MQ=255)                                                    gcCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACAGGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTc                                                                                                 <  1:229610/149‑1 (MQ=255)                                                         gCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCtt                                                                                            <  8:242701/149‑1 (MQ=255)                                                          cccccGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTg                                                                                           >  2:138761/1‑149 (MQ=255)                                                           ccccGACCAGCCGCGCTATCGCGCCGAGTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGAGTTCCCCGAGGACTACCGCCTTGGCGAAACACTc                                                                                                 <  5:41267/142‑1 (MQ=255)                                                           ccccGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTgc                                                                                          <  2:132294/149‑1 (MQ=255)                                                           ccccGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTgc                                                                                          >  5:58009/1‑149 (MQ=255)                                                           ccccGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTCGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTgc                                                                                          <  4:152193/149‑1 (MQ=255)                                                           ccccGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTCGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTgc                                                                                          <  4:147333/149‑1 (MQ=255)                                                           ccccGACCAGCCGCGCTATCGCGCCCACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCACGAGCACTACCGCCATGGCGAAACACTCgagatcg                                                                                          >  6:41267/1‑143 (MQ=255)                                                            cccGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTgcg                                                                                         >  2:184270/1‑149 (MQ=255)                                                              cGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCt                                                                                       <  1:138761/149‑1 (MQ=255)                                                                    gCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCTACACACTCGACTTTCGAAacgag                                                                                  >  6:223105/1‑147 (MQ=255)                                                                    gCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTacgac                                                                                  >  4:51373/1‑148 (MQ=255)                                                                     ccGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTacgacg                                                                                 >  2:227845/1‑148 (MQ=255)                                                                       gcgcTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAAc                                                                              <  3:211579/149‑1 (MQ=255)                                                                       gcgcTATCGCGCCGAATACGGCATGCGCCTTGCTGGCTCTACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAAGACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAAc                                                                              <  2:266107/149‑1 (MQ=255)                                                                                cgcCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCCTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACa                                                                     <  1:52785/149‑1 (MQ=255)                                                                                    gACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATAGATCCCCTACCTTGAACAAGCCACTCACGACTTCCCCGAGTACTACCGCCTTTGCGAAACACTCGCCTTGCGCTACGACGAACAAGCCGACAgcgc                                                                 >  3:285407/1‑149 (MQ=255)                                                                                            cATGCGCCTGGCTGGCTCGACAGACAAGGCTCATCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCAc                                                         <  6:203157/149‑1 (MQ=255)                                                                                             aTGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCAcg                                                        >  7:114666/1‑149 (MQ=255)                                                                                                      gctggctCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGaa                                                 <  3:51373/147‑1 (MQ=255)                                                                                                        tggctCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTgc                                             >  7:41004/1‑149 (MQ=255)                                                                                                           ctCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGcgcc                                          <  2:162904/149‑1 (MQ=255)                                                                                                            tCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGcgccg                                         >  3:247179/1‑149 (MQ=255)                                                                                                            tCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGcgccg                                         >  6:128041/1‑149 (MQ=255)                                                                                                               acagacaAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCAt                                      >  2:293402/1‑149 (MQ=255)                                                                                                                cagacaAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATc                                     <  2:320862/149‑1 (MQ=255)                                                                                                                   acaAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTc                                  >  4:319084/1‑149 (MQ=255)                                                                                                                   acaAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTc                                  >  6:2487/1‑149 (MQ=255)                                                                                                                   acaAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTc                                  >  8:226720/1‑149 (MQ=255)                                                                                                                    caAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCg                                 <  5:128041/149‑1 (MQ=255)                                                                                                                    caAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCAGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCg                                 <  6:70245/149‑1 (MQ=255)                                                                                                                     aaGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGa                                >  6:21142/1‑149 (MQ=255)                                                                                                                     aaGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGa                                >  6:25715/1‑149 (MQ=255)                                                                                                                         cTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACgca                            <  6:71702/149‑3 (MQ=255)                                                                                                                            aGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTagag                         <  5:223105/149‑1 (MQ=255)                                                                                                                            aGCGCATGCAATCGATCCCCTACATTGAACACTCCACTCACGACTTCCCCGACGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTaga                          <  8:125142/148‑1 (MQ=255)                                                                                                                             gcgcATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGc                        <  3:92009/149‑1 (MQ=255)                                                                                                                                    cAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCt                 <  1:227845/149‑1 (MQ=255)                                                                                                                                     aaTCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCTg                <  3:319084/149‑1 (MQ=255)                                                                                                                                      aTCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCTgg               >  5:86832/1‑149 (MQ=255)                                                                                                                                       tCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGAGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCTGGt              >  2:223724/1‑149 (MQ=255)                                                                                                                                         gATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCTGGTc             <  6:58009/148‑1 (MQ=255)                                                                                                                                              ccTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCTGGTCGACGGt       >  2:181454/1‑149 (MQ=255)                                                                                                                                                  ccTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCTGGTCGACGGtga     <  4:247179/147‑1 (MQ=255)                                                                                                                                                   cTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCt                 >  7:80091/1‑134 (MQ=255)                                                                                                                                                   cTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCt                 <  8:80091/134‑1 (MQ=255)                                                                                                                                                   cTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATACCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCTGGTCGACGGtgatga  >  6:38302/1‑149 (MQ=255)                                                                                                                                                    ttGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCg                                 >  4:289695/1‑117 (MQ=255)                                                                                                                                                    ttGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCTGGt              <  4:25393/136‑1 (MQ=255)                                                                                                                                                    ttGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCTGGt              >  3:25393/1‑136 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AACTGGCCGGCGACAGTGCGCAGAGCACCGCCTGGCTGGCCGAAGCCGTGCGCCTGGCCCCCGACCAGCCGCGCTATCGCGCCGACTACGGCATGCGCCTGGCTGGCTCGACAGACAAGGCTCAGCGCATGCAATCGATCCCCTACCTTGAACACGCCACTCACGACTTCCCCGAGGACTACCGCCTTGGCGAAACACTCGCCTTGCGCTACGACGAACATGCCGACAGCGCTGCGGCACGCCGCGAACTGCGCCGCATCCTCGACGTAGAGCAGAACCTGGTCGACGGTGATGA  >  NC_002947/3909484‑3909778

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