Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,873,746	C→A	7.4%	intergenic (‑127/+31)	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,746	0	C	A	7.4%	196.5 / 4.5	63	intergenic (‑127/+31)	aqpZ/PP_4283	aquaporin Z/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (26/32);  new base A (2/3);  total (28/35)
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.57e-01

ATGCGTTTACCCAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCATGTGATG  >  NC_002947/4873600‑4873893                                                                                                                                                   |                                                                                                                                                      aTGCGTTTACCCAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCtg                                                                                                                                                     >  8:5265/1‑148 (MQ=255)           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                                                                                                                                             >  5:120865/1‑141 (MQ=255)                gACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCAAGGCAGccct                                                                                                                                             <  6:120865/142‑2 (MQ=255)                      cTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAg                                                                                                                               <  5:85943/148‑1 (MQ=255)                      cTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAg                                                                                                                               <  8:322928/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                                                                                                       <  7:245180/100‑1 (MQ=255)                                                                                              acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGaa                                                                                                       >  8:245180/1‑100 (MQ=255)                                                                                              acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                       >  6:244688/1‑148 (MQ=255)                                                                                              acaAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTcc                                                       >  5:26372/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)                                                                                                                                                   |                                                                                                                                                      ATGCGTTTACCCAGGGACATGCTCATGAATCCATTCCTCAAGAAGTTGAACGCAACGAATGCGCAAGCCCAACTGCTTGATGAAGGTTTGCTCACAAGAATTTAGCAGACCTTTGGGTTATGGCCAGTAGTGGCTGGAAGGGCTGCCTTGCAGCCCCTCG‑A‑TCACGCAGGCAAACGATCAACCCTTGACGGAAACCTTGGGTTTTTTCAACAGATGCGAGAACACGGCATGCAAATCGTCCGAGGCGCTCTCCTCGTCCAGATTCAGCTTGCTGTCGATGTGGTCCATGTGATG  >  NC_002947/4873600‑4873893

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