Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,266,114	A→T	10.5%	intergenic (+31/+28)	PP_2863 → / ← PP_2864	membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,266,114	0	A	T	10.5%	178.5 / 11.0	67	intergenic (+31/+28)	PP_2863/PP_2864	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base T (4/3);  total (34/33)
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.06e-01

CACAAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGACC  >  NC_002947/3265970‑3266262                                                                                                                                                 |                                                                                                                                                     cacaAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGcc                                                                                                                                                  >  1:184920/1‑149 (MQ=255)  acaAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCcg                                                                                                                                                 <  3:155113/149‑1 (MQ=255)  acaAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCcg                                                                                                                                                 >  8:408015/1‑149 (MQ=255)  acaAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCcg                                                                                                                                                 >  7:201187/1‑149 (MQ=255)    aaGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCcgcg                                                                                                                                               <  8:365903/149‑1 (MQ=255)      gACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGa                                                                                                                                              >  8:19523/1‑148 (MQ=255)      gACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGCTTTACCCGCGaa                                                                                                                                             <  8:69670/149‑1 (MQ=255)         ccGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGCTTTACCCGCGaagag                                                                                                                                          <  8:13344/149‑4 (MQ=255)         ccGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGCTTTACCCGCGaagag                                                                                                                                          >  7:13344/1‑146 (MQ=255)          cGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAggg                                                                                                                                         <  7:408015/149‑1 (MQ=255)             aaGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCg                                                                                                                                      <  8:86116/149‑1 (MQ=38)              atgTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCgg                                                                                                                                     <  2:332842/147‑1 (MQ=255)               ggTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCgg                                                                                                                                     >  2:245210/1‑148 (MQ=255)                     gTTCCAGGTGGTGTGCTGGTGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAgg                                                                                                                               <  8:201187/148‑1 (MQ=255)                         cAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGAt                                                                                                                          >  5:326572/1‑149 (MQ=255)                           ggtggcGGGCTGGGTGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATca                                                                                                                        <  7:272232/149‑1 (MQ=255)                              ggtgTGCTGGGGGATTGTGCTGCTGCACGAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACAt                                                                                                                     <  8:264811/149‑1 (MQ=255)                               gtgtGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACAt                                                                                                                     >  1:229571/1‑148 (MQ=255)                                       ggggATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAAc                                                                                                              <  1:245210/147‑1 (MQ=255)                                          gATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCAc                                                                                                         <  1:72961/149‑1 (MQ=255)                                           aTTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCAc                                                                                                         <  6:303270/148‑1 (MQ=255)                                             tgtgCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCt                                                                                                       <  3:271993/148‑1 (MQ=255)                                              gtgCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTg                                                                                                      <  4:195571/148‑1 (MQ=255)                                               tgctgctgcACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGc                                                                                                     >  3:351761/1‑148 (MQ=255)                                                  tgctgcACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAAcc                                                                                                 >  1:232903/1‑149 (MQ=255)                                                  tgctacACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGCCCAACCCCACCTGCAAcc                                                                                                 >  2:88260/1‑149 (MQ=255)                                                     tgcACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCt                                                                                              >  5:184747/1‑149 (MQ=255)                                                       cACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTg                                                                                            >  1:311073/1‑149 (MQ=255)                                                       cACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTg                                                                                            <  6:326572/149‑1 (MQ=255)                                                          cAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATcaca                                                                                                                      >  1:373395/1‑120 (MQ=255)                                                          cAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATcaca                                                                                                                      <  2:373395/120‑1 (MQ=255)                                                          cAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTggg                                                                                          >  4:178529/1‑148 (MQ=255)                                                                ttCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGa                                                                                   >  1:50973/1‑149 (MQ=255)                                                                           cTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAg                                                                        >  6:344991/1‑149 (MQ=255)                                                                                tGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGt                                                                   >  1:281270/1‑149 (MQ=255)                                                                                     cTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACg                                                                                 <  3:334651/130‑1 (MQ=255)                                                                                     cTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACg                                                                                 >  4:334651/1‑130 (MQ=255)                                                                                     cTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGAc                                                              >  3:266527/1‑149 (MQ=255)                                                                                       tGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCg                                                            >  4:217018/1‑149 (MQ=255)                                                                                         aCGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGc                                                           <  5:336709/148‑1 (MQ=255)                                                                                           gggCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGc                                                        <  4:351761/149‑1 (MQ=255)                                                                                             gCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCaa                                                      >  1:237938/1‑149 (MQ=255)                                                                                             gCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCaa                                                      <  2:50973/149‑1 (MQ=255)                                                                                                       ccccTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACCCCTcaca                                             >  1:274589/1‑148 (MQ=255)                                                                                                       ccccTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTc                                                <  2:229571/145‑1 (MQ=255)                                                                                                         ccTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAg                                          >  3:103284/1‑149 (MQ=255)                                                                                                             ggTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGtt                                        <  2:311073/147‑1 (MQ=255)                                                                                                                          ttGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                         <  2:237938/149‑1 (MQ=255)                                                                                                                          ttGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                         >  2:11621/1‑149 (MQ=255)                                                                                                                           tGTTCCGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCGGTACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACCCAGTTCATCACGCCCCATggg                        >  7:48456/1‑149 (MQ=37)                                                                                                                              aCCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGcc                               >  5:69112/1‑139 (MQ=255)                                                                                                                              aCCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGcc                               <  6:69112/139‑1 (MQ=255)                                                                                                                              aCCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgggg                       >  2:127451/1‑147 (MQ=255)                                                                                                                               ccGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAg                    <  3:178529/149‑1 (MQ=255)                                                                                                                                 ggccctttCGCGGCTTTACCCGCGAAGAGGCCGGTACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAggg                  >  5:390533/7‑149 (MQ=38)                                                                                                                                  gccctttCGCGGCTTTACCCGCGAAGAGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGggc                 >  7:43653/6‑149 (MQ=255)                                                                                                                                   ccTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGggcg                >  6:297093/1‑149 (MQ=255)                                                                                                                                     tctTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAggg                  >  5:37044/1‑145 (MQ=255)                                                                                                                                     tctTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAggg                  <  6:37044/145‑1 (MQ=255)                                                                                                                                       ttCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGggcggc              <  7:19523/147‑1 (MQ=255)                                                                                                                                         cgcgGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTcaca                                             <  4:231540/114‑1 (MQ=255)                                                                                                                                         cgcgGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTcaca                                             >  3:231540/1‑114 (MQ=255)                                                                                                                                         cgcgGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGgcggcg             <  2:281270/146‑1 (MQ=255)                                                                                                                                            gggTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGc       <  2:184920/149‑1 (MQ=255)                                                                                                                                               tAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGa    >  8:245657/1‑149 (MQ=255)                                                                                                                                               tAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGa    <  5:344991/149‑1 (MQ=255)                                                                                                                                               tAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGa    <  2:232903/149‑1 (MQ=255)                                                                                                                                                aaaGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGAc   >  6:319487/1‑149 (MQ=255)                                                                                                                                                 aaGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGAcc  <  4:103284/149‑1 (MQ=255)                                                                                                                                                 aaGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGAc   <  3:217018/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CACAAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGACC  >  NC_002947/3265970‑3266262

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