Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,152,779	A→G	5.2%	D168G (GAT→GGT)	PP_1907 →	haloacid dehalogenase‑like family hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,152,779	0	A	G	5.2%	230.5 / 2.2	76	D168G (GAT→GGT)	PP_1907	haloacid dehalogenase‑like family hydrolase
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base G (2/2);  total (39/37)
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.64e-01

TCGGGTGCTCAAGGCCAATGGCTGGGAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTATGCATCGATCACTTTTCCCAGTTGCGCGAGTGGCTGGGT  >  NC_002947/2152633‑2152924                                                                                                                                                     |                                                                                                                                                  tCGGGTGCTCAAGGCCAATGGCTGGGAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGAtt                                                                                                                                                 >  6:78312/1‑149 (MQ=255) tCGGGTGCTCAAGGCCAATGGCTGGGAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGAtt                                                                                                                                                 <  3:281681/149‑1 (MQ=255)   gggTGCTCAAGGCCAATGGCTGGGAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTc                                                                                                                                                <  5:78312/148‑1 (MQ=255)         tCAAGGCCAATGGCTGGGAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTc                                                                                                                                         >  3:272678/1‑149 (MQ=255)         tCAAGGCCAATGGCTGGGAGCGGTTCTTCAACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGGTTCGGCATTc                                                                                                                                         >  2:171076/1‑149 (MQ=255)         tCAAGGCCAATGGCTGGGAGCGGTTCTTCAACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGGTTCGGCATTc                                                                                                                                         >  2:165576/1‑149 (MQ=255)            aGGCCAATGGCTGGGAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGAc                                                                                                                                      <  3:102732/149‑1 (MQ=255)             ggCTAATGGCTGGGAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACt                                                                                                                                     <  4:156727/149‑1 (MQ=255)                     gCTGGGAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGAt                                                                                                                              <  2:67795/148‑1 (MQ=255)                      cTGGGAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATgg                                                                                                                            >  1:19617/1‑149 (MQ=255)                      cTGGGAGAGGTTCTTCGACATCACCCCCGCCGCTGCCGAGACACGTGGCAACCAACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTGGGCATTCTACATGGAGATgg                                                                                                                            >  4:301789/1‑149 (MQ=255)                         ggAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACC‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGc                                                                                                                           >  2:320443/1‑147 (MQ=255)                         ggAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCa                                                                                                                         <  4:272678/149‑1 (MQ=255)                         ggAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCAGTCGACTTGCAGATGGccc                                                                                                                         >  6:43666/1‑148 (MQ=255)                               gTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATg                                                                                                                   >  6:236826/1‑149 (MQ=255)                               gTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATg                                                                                                                   >  6:242254/1‑149 (MQ=255)                               gTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATg                                                                                                                   >  4:31496/1‑149 (MQ=255)                                     ctaCATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGGTTCGGCATTCGACTTGCAGATGGCCAGCAATGCCgg                                                                                                               <  1:171076/145‑1 (MQ=255)                                     cGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGGTTCGGCATTCGACTTGCAGATGGCCAGCAATGCCgg                                                                                                               <  1:165576/147‑1 (MQ=255)                                     cGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGAtt                                                                                                                                                 >  8:220136/1‑113 (MQ=255)                                       aCATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTc                                                                                                                                         <  8:327142/119‑1 (MQ=255)                                       aCATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTc                                                                                                                                         >  7:327142/1‑119 (MQ=255)                                         aTCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGgcatgca                                                                                                         >  1:354932/1‑149 (MQ=255)                                             cccGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGGTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTcc                                                                                                     >  5:48564/1‑148 (MQ=255)                                             cccGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCg                                                                                                     >  1:147766/1‑149 (MQ=255)                                                gcgcCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTg                                                                                                  <  2:354932/149‑1 (MQ=255)                                                  gccgcTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGc                                                                                                <  4:5739/149‑1 (MQ=255)                                                     gcTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGGACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGt                                                                                             >  3:220321/1‑149 (MQ=255)                                                         aCGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGACGAAATCCTCTGGCATTTTGGCGTCGAACCGTCGCGTGCGTTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTgg                                                                                          <  8:351784/148‑1 (MQ=255)                                                            agaCCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTAt                                                                                      <  5:242254/149‑1 (MQ=255)                                                            agaCCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTAt                                                                                      <  5:236826/149‑1 (MQ=255)                                                             gaCCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATg                                                                                     <  1:271929/149‑1 (MQ=255)                                                              gggtgTGGCTTGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATgg                                                                                    >  6:175912/5‑149 (MQ=255)                                                               cccGTGTCAATCCCCATCCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑TAGATTCGGCNTTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATg                                                                                     <  1:316566/148‑1 (MQ=255)                                                                ccGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTg                                                                                  >  5:236082/1‑149 (MQ=255)                                                                ccGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTg                                                                                  >  2:253104/1‑149 (MQ=255)                                                                  gTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGGACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGcc                                                                                <  4:220321/149‑1 (MQ=255)                                                                  gTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGc                                                                                 >  4:223863/1‑148 (MQ=255)                                                                    ggCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGgcat                                                                                                            <  5:213930/119‑1 (MQ=255)                                                                    ggCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGgcat                                                                                                            >  6:213930/1‑119 (MQ=255)                                                                     gCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATg                                                                             >  1:159719/1‑149 (MQ=255)                                                                          ccACA‑GCCGCTGATGCTTGAGGAATTCCTCGGGCATTGTGGGGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGcc                                                                                <  8:142985/141‑1 (MQ=255)                                                                          ccACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGcc                                                                                >  7:142985/1‑141 (MQ=255)                                                                           caca‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTg                                                                       >  7:316965/1‑149 (MQ=255)                                                                            aca‑cCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTg                                                                       >  3:270679/1‑148 (MQ=255)                                                                                   gCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGc                                                                         >  3:172679/1‑140 (MQ=255)                                                                                   gCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGc                                                                         <  4:172679/140‑1 (MQ=255)                                                                                      gATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTgg                                                             >  1:342568/1‑149 (MQ=255)                                                                                        tGCTTGAGGAAATCCTCGGGCATTGTGGCGTCTAACCGGCGCGTGCGCGGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGcc                                                           <  6:48564/149‑1 (MQ=255)                                                                                         gCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTgg                                                             >  4:301635/1‑146 (MQ=255)                                                                                         gCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTgg                                                             <  3:301635/146‑1 (MQ=255)                                                                                            tGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGt                                                       >  3:128452/1‑149 (MQ=255)                                                                                               ggAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTg                                                    <  6:236082/149‑1 (MQ=255)                                                                                               ggAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTg                                                    >  4:16390/1‑149 (MQ=255)                                                                                                gAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCg                                                                                                                                               >  4:277316/1‑57 (MQ=255)                                                                                                gAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCg                                                                                                                                               <  3:277316/57‑1 (MQ=255)                                                                                                gAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGttt                                                     <  1:320443/147‑1 (MQ=255)                                                                                                    tCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCg                                               <  4:128452/149‑1 (MQ=255)                                                                                                       tCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAg                                            >  1:229450/1‑149 (MQ=255)                                                                                                       tCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAg                                            >  2:232062/1‑149 (MQ=255)                                                                                                       tCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGAAGGCGCTGGCCGAGTTTGGCCCGCAg                                            >  4:340932/1‑149 (MQ=255)                                                                                                         gggCATTGTGGCGTCGAACCGGCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAg                                            <  3:301789/147‑1 (MQ=255)                                                                                                          ggCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCGCCCGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGt                                          <  1:125766/149‑1 (MQ=255)                                                                                                           gCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTAt                                        <  4:270679/149‑1 (MQ=255)                                                                                                                gtgGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTATGCAt                                    <  2:229450/148‑1 (MQ=255)                                                                                                                   gCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTATGCATCGAt                                <  3:16390/149‑1 (MQ=255)                                                                                                                     gTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTATGCATCGATCa                              <  1:253104/149‑1 (MQ=255)                                                                                                                            cGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCAtt                                                                                                       >  7:271941/1‑69 (MQ=255)                                                                                                                            cGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCAtt                                                                                                       <  8:271941/69‑1 (MQ=255)                                                                                                                                gcgTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTATGCATCGATCACTTTTCCCAg                    <  3:223863/148‑1 (MQ=255)                                                                                                                                   tGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGAGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTATGCATCGATCACTTTTCCCAGTTgc                <  3:31496/149‑1 (MQ=255)                                                                                                                                    gcgcTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTATGCATCGATCACTTTTCCCAGTTgcg               <  8:316965/149‑1 (MQ=255)                                                                                                                                        tGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGACGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTATGCATCGATCACTTTTCCCAGTTGCGCGAg           >  2:296072/1‑149 (MQ=255)                                                                                                                                            gttCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGACGCTGGCCGAGTTTGGCCCGCAGGTATGCATCGATCACTTTTCCCAGTTGCGCGAGTGGc       >  1:113954/3‑149 (MQ=255)                                                                                                                                            ggTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTATGCATCGATCACTTTTCCCAGTTGCGCGAGTGGc       >  8:145008/1‑149 (MQ=255)                                                                                                                                               cG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTATGCATCGATCACTTTTCCCAGTTGCGCGAGTGGCTgg    <  2:342568/149‑1 (MQ=255)                                                                                                                                                  gagaTTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTATGCATCGATCACTTTTCCCAGTTGCGCGAGTGGCTGggt  <  2:159719/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TCGGGTGCTCAAGGCCAATGGCTGGGAGCGGTTCTTCGACATCACCCGCGCCGCTGACGAGACCCGTGGCAAGCCACA‑CCCGCTGATGCTTGAGGAAATCCTCGGGCATTGTGGCGTCGAACCGTCGCGTGCGCTGATGGTCG‑GAGATTCGGCATTCGACTTGCAGATGGCCAGCAATGCCGGCATGCATTCGGTGGCGGTTGGCTATGGTGCCATGTCGCTGCAGGCGCTGGCCGAGTTTGGCCCGCAGGTATGCATCGATCACTTTTCCCAGTTGCGCGAGTGGCTGGGT  >  NC_002947/2152633‑2152924

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