Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	222,264	A→T	5.5%	Q40L (CAA→CTA) ‡	PP_0170 →	ABC transporter substrate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	222,264	0	A	T	5.5%	278.8 / 4.4	91	Q40L (CAA→CTA) ‡	PP_0170	ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (42/44);  new base T (2/3);  total (44/47)
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.02e-01

CCCGTTCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGG  >  NC_002947/222119‑222408                                                                                                                                                     |                                                                                                                                                 cccGTTCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAAtt                                                                                                                                               >  7:212739/1‑149 (MQ=255) cccGTTCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGCGCGAAAGCCTG‑GACGTCGGCACCCAAGCCGCCCCCGCCGTGGCCCAGCCCGAACGCAAGATC‑CGTATTGCCGA‑ACaaat                                                                                                                                               >  1:19478/1‑147 (MQ=255)  ccGTTCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTc                                                                                                                                              <  1:148224/149‑1 (MQ=255)  ccGTTCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTc                                                                                                                                              >  2:154820/1‑149 (MQ=255)   cGTTCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCg                                                                                                                                             <  5:41958/149‑1 (MQ=255)    gTTCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCgg                                                                                                                                            >  6:32997/1‑149 (MQ=255)        tGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATc                                                                                                                                        >  6:300284/1‑149 (MQ=255)          acgacTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ATTGTTCGGCaata                                                                                                                                       <  2:277729/148‑4 (MQ=255)          acgacTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCg                                                                                                                                       <  5:300284/148‑1 (MQ=255)            gacTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTc                                                                                                                                     >  8:37401/1‑148 (MQ=255)              cTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTAc                                                                                                                                  >  8:336732/1‑149 (MQ=255)                 aGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTGAAGCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACct                                                                                                                                >  1:295615/1‑149 (MQ=255)                 aGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACctg                                                                                                                               >  5:208758/1‑149 (MQ=255)                 aGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACct                                                                                                                                <  3:31764/148‑1 (MQ=255)                   gAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACctgct                                                                                                                             >  3:11371/1‑149 (MQ=255)                     aaCCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACctgctgc                                                                                                                           <  3:127170/149‑1 (MQ=255)                      aCCATCATTCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATCTCGTATTGCCGA‑ACAATTCGGCATCGTCTACctgctgc                                                                                                                           <  1:251700/149‑1 (MQ=255)                        catcatGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACctgctgctg                                                                                                                         >  6:72185/1‑148 (MQ=255)                        catcatGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACctgctgctg                                                                                                                         >  7:313865/1‑148 (MQ=255)                         atcatGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGaa                                                                                                                       <  6:208758/149‑1 (MQ=255)                             tGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACgtg                                                                                                                   <  5:32997/149‑1 (MQ=255)                               cgcAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACgtggt                                                                                                                 <  7:336732/149‑1 (MQ=255)                                gcAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACgtggtg                                                                                                                >  8:176245/1‑149 (MQ=255)                                 cAATTCCATCAGCCGCGTGGCGGCAAGCTTCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGGCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACATTTCGGCATCGTCTACCTGCTTCTGAACGTGGTGc                                                                                                               <  2:19478/149‑1 (MQ=255)                                 cAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGc                                                                                                               <  4:64641/149‑1 (MQ=255)                                       cATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGAAAGATA‑CGTATTGCCGA‑ACAATTAGGCATAGTCTACCTGCTGCTGAACGTGGTGCGTGAcc                                                                                                         >  6:358564/1‑149 (MQ=255)                                        aTCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACgtggt                                                                                                                 <  2:292677/140‑1 (MQ=255)                                        aTCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACgtggt                                                                                                                 >  1:292677/1‑140 (MQ=255)                                        aTCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCATGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGAc                                                                                                          >  7:237957/1‑147 (MQ=255)                                         tCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAg                                                                                                       <  2:188751/149‑1 (MQ=255)                                           aGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTT‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGc                                                                                                      <  7:37401/148‑1 (MQ=255)                                           aGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGc                                                                                                      <  3:352754/148‑1 (MQ=255)                                               gccTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCt                                                                                                  >  3:31100/1‑148 (MQ=255)                                                ccTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGa                                                                                                <  8:212739/149‑1 (MQ=255)                                                 cTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGAt                                                                                               <  4:11371/149‑1 (MQ=255)                                                  tGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGAtt                                                                                              >  2:120718/1‑149 (MQ=255)                                                      ggcAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTgaga                                                                                          <  1:120718/149‑1 (MQ=255)                                                        cAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGaaa                                                                                        >  3:253454/1‑149 (MQ=255)                                                         aaGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGaaa                                                                                        >  7:170620/1‑148 (MQ=255)                                                         aaGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAac                                                                                       >  3:169929/1‑149 (MQ=255)                                                          aGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCCGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAaca                                                                                      >  8:290319/1‑149 (MQ=255)                                                            cATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACg                                                                                    <  1:154820/149‑1 (MQ=255)                                                             aTCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCCGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACgg                                                                                   <  7:290319/149‑1 (MQ=255)                                                                 gCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAAc                                                                                                                      <  2:184967/110‑1 (MQ=255)                                                                 gCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAAc                                                                                                                      >  1:184967/1‑110 (MQ=255)                                                                 gCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAg                                                                               <  4:343080/149‑1 (MQ=255)                                                                 gCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAg                                                                               >  5:13097/1‑149 (MQ=255)                                                                  ccTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAacac                                                                                     >  4:201712/1‑142 (MQ=255)                                                                  ccTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAacac                                                                                     <  3:201712/142‑1 (MQ=255)                                                                  ccTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCATGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAgg                                                                              <  8:237957/149‑1 (MQ=255)                                                                   cTGGGCGCAAGCCTG‑GTCGTCGGCAGCATAGCCGCCCCCGCCGTTGACCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTTTTGAACGTTGTTCGTGACCAGCACCTGATTGAGAAACACGGCAATGa                                                                             <  5:358564/149‑1 (MQ=255)                                                                      ggCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCa                                                                          <  7:126226/149‑1 (MQ=255)                                                                        cgcAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAgg                                                                        <  4:350773/149‑1 (MQ=255)                                                                        cgcAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAgg                                                                        <  7:135227/149‑1 (MQ=255)                                                                        cgcAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAgg                                                                        >  3:350773/1‑149 (MQ=255)                                                                          cAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGc                                                                      <  7:39228/149‑1 (MQ=255)                                                                            aGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGcat                                                                    <  7:187470/149‑1 (MQ=255)                                                                             gCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGcatc                                                                   >  2:53527/1‑149 (MQ=255)                                                                               cTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGcatcga                                                                 >  5:243246/1‑149 (MQ=255)                                                                                    tcgtcgGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑AGAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGcatcgacatc                                                             <  8:313865/149‑1 (MQ=255)                                                                                         ggCAGCCTAGCCGCCCCCGCCGTGTCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGt                                                        <  5:166907/149‑1 (MQ=255)                                                                                         ggCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCTTCGACATCGAAGt                                                        <  5:72185/149‑1 (MQ=255)                                                                                          gCATCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTc                                                       <  6:178417/149‑1 (MQ=255)                                                                                            aGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGa                                                     <  8:170620/149‑1 (MQ=255)                                                                                                tAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTgg                                                 >  1:202118/1‑149 (MQ=255)                                                                                                  gccgccCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTgg                                                 <  4:253454/147‑1 (MQ=255)                                                                                                       ccccGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGAGACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAg                                                                            >  1:71624/1‑116 (MQ=255)                                                                                                       ccccGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAg                                                                            <  2:71624/115‑1 (MQ=255)                                                                                                        cccGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCt                                         >  6:260132/1‑149 (MQ=255)                                                                                                         ccgccgTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGaa                                                          >  6:54571/1‑131 (MQ=255)                                                                                                         ccgccgTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGaa                                                          <  5:54571/131‑1 (MQ=255)                                                                                                            ccgTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCg                                     <  4:31100/149‑1 (MQ=255)                                                                                                               tGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGt                                  <  7:176245/149‑1 (MQ=255)                                                                                                                   ccAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTga                                                                                            >  6:286053/1‑87 (MQ=255)                                                                                                                   ccAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTga                                                                                            <  5:286053/87‑1 (MQ=255)                                                                                                                   ccAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCt                              >  2:96071/1‑149 (MQ=255)                                                                                                                    cAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTc                             >  8:191037/1‑149 (MQ=255)                                                                                                                       gCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGc                          <  2:295615/149‑1 (MQ=255)                                                                                                                       gCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGAGTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGAAGCTGTCGGGTGGCTCGGc                          >  5:209024/1‑149 (MQ=255)                                                                                                                             ggCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCaa                    >  2:7204/1‑149 (MQ=255)                                                                                                                                        gTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCAttg          >  8:64748/1‑149 (MQ=255)                                                                                                                                          gatGCCGA‑ATTGTTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGAtt                                                   <  7:338772/104‑1 (MQ=255)                                                                                                                                          gatGCCGA‑ATTGTTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGAtt                                                   >  8:338772/3‑106 (MQ=255)                                                                                                                                           ttGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCAttgtt        <  6:13097/148‑1 (MQ=255)                                                                                                                                           atGCCGA‑ATTGTTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGc                               <  7:277330/124‑1 (MQ=255)                                                                                                                                           atGCCGA‑ATTGTTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGc                               >  8:277330/2‑125 (MQ=255)                                                                                                                                              ccGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGcc                         <  2:131293/128‑1 (MQ=255)                                                                                                                                              ccGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGcc                         >  1:131293/1‑128 (MQ=255)                                                                                                                                              ccGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTcc    >  8:119099/1‑149 (MQ=255)                                                                                                                                              ccGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTcc    >  2:213841/1‑149 (MQ=255)                                                                                                                                               cGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCg   >  4:155103/1‑149 (MQ=255)                                                                                                                                                gA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCgg  >  4:172318/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CCCGTTCTGACGACTGAGGAAACCATCATGCGCAAATCCATCAGCCGCCTGGCGGCAAGCATCGGCCTGGGCGCAAGCCTG‑GTCGTCGGCAGCCTAGCCGCCCCCGCCGTGGCCCAGGCCGAAGGCAAGATC‑CGTATTGCCGA‑ACAATTCGGCATCGTCTACCTGCTGCTGAACGTGGTGCGTGACCAGCACCTGATTGAGAAACACGGCAAGGAGCAGGGCATCGACATCGAAGTCGATTGGGCGCAGCTGTCGGGTGGCTCGGCCGTCAACGATGCATTGTTGTCCGG  >  NC_002947/222119‑222408

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