Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,873,748	T→G	7.9%	intergenic (‑129/+29)	aqpZ ← / ← PP_4283	aquaporin Z/GntR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,873,748	0	T	G	7.9%	192.0 / 7.4	63	intergenic (‑129/+29)	aqpZ/PP_4283	aquaporin Z/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (25/33);  new base G (2/3);  total (27/36)
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.34e-01

CCAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCATGTGATGCA  >  NC_002947/4873610‑4873895                                                                                                                                           |                                                                                                                                                      ccAGGGCCATGCTCATGAATCCATTCCTCAAGAAGATGAACGCAACGAATGCGCAAGCCCAAATGATTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGTGCTGCCTTGCAGCCCCTc                                                                                                                                             >  4:182196/1‑149 (MQ=255)  cAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCtt                                                                                                                                                       >  3:170029/1‑138 (MQ=255)  cAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCtt                                                                                                                                                       <  4:170029/138‑1 (MQ=255)    gggACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGCCCTTccag                                                                                                                                          <  3:343384/149‑4 (MQ=255)      gACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGccct                                                                                                                                               <  6:120865/142‑2 (MQ=255)      gACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGccct                                                                                                                                               >  5:120865/1‑141 (MQ=255)            cTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAg                                                                                                                                 <  8:322928/148‑1 (MQ=255)            cTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAg                                                                                                                                 <  5:85943/148‑1 (MQ=255)                      cATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATc                                                                                                                       >  5:78830/1‑148 (MQ=255)                      cATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATc                                                                                                                       >  2:135553/1‑148 (MQ=255)                      cATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATc                                                                                                                       <  6:78830/148‑1 (MQ=255)                       aTTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑a                                                                                                                                          >  4:36653/1‑129 (MQ=255)                       aTTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑a                                                                                                                                          <  3:36653/129‑1 (MQ=255)                       aTTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCaa                                                                                                                     <  7:173663/149‑1 (MQ=255)                       aTTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCaa                                                                                                                     <  1:186639/149‑1 (MQ=255)                           cTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAAccc                                                                                                                  <  1:91292/148‑1 (MQ=255)                            tCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGCCCTTc                                                                                                                                             >  4:62705/1‑122 (MQ=255)                            tCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGCCCTTc                                                                                                                                             <  3:62705/122‑1 (MQ=255)                              aagaagTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCtt                                                                                                                <  5:65388/147‑1 (MQ=255)                              aagaagTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGa                                                                                                              >  3:183143/1‑149 (MQ=255)                                   gTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCa                                                                                                                                      >  1:153462/1‑120 (MQ=255)                                   gTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCa                                                                                                                                      <  2:153462/120‑1 (MQ=255)                                        aCGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCtt                                                                                                    <  5:33006/149‑1 (MQ=255)                                         cGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTg                                                                                                   <  6:233310/149‑1 (MQ=255)                                          gCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCaa                                                                                                                     <  2:222458/130‑1 (MQ=255)                                          gCAACGAATGCGCAAGCCCAACTGCTTGATGAAAGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCaa                                                                                                                     >  1:222458/1‑130 (MQ=255)                                            aaCGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGt                                                                                                >  2:289705/1‑149 (MQ=255)                                            aaCGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGt                                                                                                >  2:213832/1‑149 (MQ=255)                                                 aTGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATc                                                                                                                       >  2:89010/1‑121 (MQ=255)                                                 aTGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATc                                                                                                                       <  1:89010/121‑1 (MQ=255)                                                    cgcAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCaa                                                                                        <  4:183143/149‑1 (MQ=255)                                                             aaCTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCGAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCT‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTgg                                                                                                  <  6:286220/130‑1 (MQ=255)                                                             aaCTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTgg                                                                                                  >  5:286220/1‑130 (MQ=255)                                                             aaCTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCga                                                                               <  7:104354/149‑1 (MQ=255)                                                             aaCTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCga                                                                               <  5:326139/149‑1 (MQ=255)                                                             aaCTGCTTGATGAAGGTTTGCTCAAAAGAATTTAGCAGACCTTTGGGTTATTGCAATTAGTGGCTGGAATGGCTTCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCga                                                                               <  3:182196/149‑1 (MQ=255)                                                             aaCGGCTTGATGAAGGTTTGCTCACAAGATTTGAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCGTGCAGCCCCTCG‑A‑TCCCGCAGGCATACGATCCACCCTTGACGGAAACCTTGGGTTTTTTCAAGAGATGCga                                                                               <  4:138770/149‑1 (MQ=255)                                                                tGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGaa                                                                            >  7:291440/1‑149 (MQ=255)                                                                   ttGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAacac                                                                         >  6:313691/1‑149 (MQ=255)                                                                        gaAGGTTTGCTCACAAGAATTAAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGCAAGGGCTGCCTTGCAGACACTCG‑A‑TCAAGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTATTTTCAACAGATGCGAGAACACGGCAt                                                                    >  1:263015/1‑149 (MQ=255)                                                                                  tCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAAc                                                                                       <  2:124816/120‑1 (MQ=255)                                                                                    acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGaa                                                                                                         >  8:245180/1‑100 (MQ=255)                                                                                    acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGaa                                                                                                         <  7:245180/100‑1 (MQ=255)                                                                                    acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                         >  5:26372/1‑148 (MQ=255)                                                                                    acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                         >  6:244688/1‑148 (MQ=255)                                                                                      aaGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAg                                                      <  1:213832/149‑1 (MQ=255)                                                                                               gCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTctcc                                             <  8:291440/149‑1 (MQ=255)                                                                                                cAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCga                                                                               <  7:318269/114‑1 (MQ=255)                                                                                                cAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCga                                                                               >  8:318269/1‑114 (MQ=255)                                                                                                cAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTctcct                                            >  3:262689/1‑149 (MQ=255)                                                                                                 aGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTctcctc                                           >  8:310479/1‑149 (MQ=255)                                                                                                         gggTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGa                                    <  5:244688/148‑1 (MQ=255)                                                                                                                 ggcAGTTGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACTCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTg                           <  2:263015/147‑1 (MQ=255)                                                                                                                 gCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTg                           <  1:135553/149‑1 (MQ=255)                                                                                                                  ccAGTATTGGCTGGAATGGCTGCCTTTGAGCCCCTCGCACTCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCtt                            <  1:289705/149‑1 (MQ=255)                                                                                                                     gtagtGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGt                       >  7:324276/1‑149 (MQ=255)                                                                                                                          ggCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAtg                  >  1:316016/1‑149 (MQ=255)                                                                                                                           gCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAtg                  <  2:316016/148‑1 (MQ=255)                                                                                                                            cTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCAGGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAtgtg                >  1:213012/1‑149 (MQ=255)                                                                                                                            cTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAtgtg                <  7:310479/149‑1 (MQ=255)                                                                                                                             tGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTgg               >  6:120358/1‑149 (MQ=255)                                                                                                                                 aGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCa           <  6:26372/149‑1 (MQ=255)                                                                                                                                        ccTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCATGTGATg    <  5:313691/149‑1 (MQ=255)                                                                                                                                         cTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAgg                                                     >  7:97787/1‑99 (MQ=255)                                                                                                                                         cTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAgg                                                     <  8:97787/99‑1 (MQ=255)                                                                                                                                          ttGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAg                                                      <  3:223236/97‑1 (MQ=255)                                                                                                                                          ttGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAg                                                      >  4:223236/1‑97 (MQ=255)                                                                                                                                          ttGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGCTTATTTCAACAGCTGCGAGAACACCGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGAAGTGGGCCATGTGATGca  >  7:257376/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                      CCAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCATGTGATGCA  >  NC_002947/4873610‑4873895

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