Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,053,543	Δ1 bp	6.2%	intergenic (‑77/‑114)	PP_5305 ← / → exbB	epimerase/dehydratase‑like protein/biopolymer transport protein ExbB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,053,543	0	A	.	6.2%	198.3 / 7.2	64	intergenic (‑77/‑114)	PP_5305/exbB	epimerase/dehydratase‑like protein/biopolymer transport protein ExbB
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base . (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

GCCAACATCTGCTTGGCCAGACGCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCC  >  NC_002947/6053396‑6053685                                                                                                                                                    |                                                                                                                                                gCCAACATCTGCTTGGCCAGACGCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTa                                                                                                                                                 >  5:211056/1‑148 (MQ=255)  ccAACATCTGCTTGGCCAGACGCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTAc                                                                                                                                                >  7:216320/1‑148 (MQ=255)   cAACATCTGCTTGGCCAGACGCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGATACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACaa                                                                                                                                              >  1:152030/1‑149 (MQ=255)      cATCTGCTTGGCCAGACGCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAatat                                                                                                                                           <  3:56720/149‑1 (MQ=255)      cATCTGCTTGGCCAGACGCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAatat                                                                                                                                           <  8:176276/149‑1 (MQ=255)        tCTGCTTGGCCAGACGCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACCCGATGGTCGGCTAGCCCATCGGCGAAAAAAGCTACAAAAtt                                                                                                                                          >  8:114520/1‑148 (MQ=255)               ggCCAGACGAACACCCACATCAACACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTATACCAATTGATCAGCTTCGCCACTCTACACTATTGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACtttt                                                                                                                                    <  7:239988/147‑1 (MQ=255)               gcccAGACGCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGt                                                                                                                                  <  1:285611/147‑1 (MQ=255)                gCCAGACGCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGt                                                                                                                                  <  5:101852/148‑1 (MQ=255)                  cAGACGCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTaa                                                                                                                               <  4:177720/149‑1 (MQ=255)                   aGACGCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTaa                                                                                                                               >  8:253934/1‑148 (MQ=255)                       gCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCTAAAAAAGCTCTAATATTACTTGTGTTAACAAg                                                                                                                           >  2:303401/1‑148 (MQ=255)                         ccacccacATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAAt                                                                                                                        >  8:83098/1‑149 (MQ=255)                           acccacATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTa                                                                                                                      <  7:204264/149‑1 (MQ=255)                               acaTCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTg                                                                                                                  <  3:15913/149‑1 (MQ=255)                                      cacaACCCTCCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCaataat                                                                                                           <  7:264254/149‑1 (MQ=255)                                            ccaccaccaTCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGc                                                                                                     <  6:211056/149‑1 (MQ=255)                                                  ccctCAAGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAAtt                                                                                               <  7:114520/146‑1 (MQ=255)                                                      cACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTtc                                                                                           >  2:64918/1‑149 (MQ=255)                                                       aCGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTtc                                                                                           <  7:253934/148‑1 (MQ=255)                                                         ggAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGAcc                                                                                                       >  1:77005/1‑134 (MQ=255)                                                         ggAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGAcc                                                                                                       <  2:77005/134‑1 (MQ=255)                                                           aaaGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGAcc                                                                                                       <  3:92478/132‑1 (MQ=255)                                                           aaaGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGAcc                                                                                                       >  4:92478/1‑132 (MQ=255)                                                               gTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGATACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGc                                                                                      >  1:295577/1‑145 (MQ=255)                                                               gTCTGACATCTCTAAACTCCCTAGACCAATTGATCACCTTCGCCACGCAACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCAACAATAATACTTTTGTAAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGc                                                                                      <  2:295577/145‑1 (MQ=255)                                                                  tGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGAcc                                                                                                       <  2:74946/125‑1 (MQ=255)                                                                  tGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGAcc                                                                                                       >  1:74946/1‑125 (MQ=255)                                                                         tctAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCaa                                                                                                               <  6:294205/110‑1 (MQ=255)                                                                         tctAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCaa                                                                                                               >  5:294205/1‑110 (MQ=255)                                                                         tctAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCg                                                                                        <  2:55831/133‑1 (MQ=255)                                                                         tctAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCg                                                                                        >  1:55831/1‑133 (MQ=255)                                                                           tAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTg                                                                      <  7:83098/149‑1 (MQ=255)                                                                             aaCTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTg‑g                                                                    <  8:216320/148‑1 (MQ=255)                                                                                tCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGaa                                                                >  8:83007/1‑149 (MQ=255)                                                                                           aTTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCt                                                            >  3:48777/1‑142 (MQ=255)                                                                                           aTTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCt                                                            <  4:48777/142‑1 (MQ=255)                                                                                                    cTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTtcatc                                            <  7:25531/149‑1 (MQ=255)                                                                                                         ccACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGt                                        <  1:131753/148‑1 (MQ=255)                                                                                                         ccACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGt                                        >  2:131753/1‑148 (MQ=255)                                                                                                         ccACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTc                                       >  8:247613/1‑149 (MQ=255)                                                                                                          cACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCT‑‑AATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑ga                                                                   >  4:162632/1‑118 (MQ=255)                                                                                                          cACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCT‑‑AATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑ga                                                                   >  4:158264/1‑118 (MQ=255)                                                                                                          cACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCT‑‑AATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑ga                                                                   <  3:158264/118‑1 (MQ=255)                                                                                                          cACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCT‑‑AATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑ga                                                                   <  3:162632/118‑1 (MQ=255)                                                                                                                    gATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCAt                             >  6:240615/1‑148 (MQ=255)                                                                                                                        ggcggcTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTGAGAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGaca                         <  1:64918/149‑1 (MQ=255)                                                                                                                         gcggcTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGacac                        <  7:247613/148‑1 (MQ=255)                                                                                                                              tAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACt                   >  2:199273/1‑148 (MQ=255)                                                                                                                              tAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCATTG‑GAGAACGCTCACTTTTTTCTTCAAAAGTTCCGTCCAGCCTGACACGAACt                   >  3:126710/1‑148 (MQ=255)                                                                                                                                gACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAg                                <  3:238400/133‑1 (MQ=255)                                                                                                                                gACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAg                                >  4:238400/1‑133 (MQ=255)                                                                                                                                gACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCaa                >  8:193020/1‑149 (MQ=255)                                                                                                                                   cATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTcc                                      >  3:210237/1‑124 (MQ=255)                                                                                                                                   cATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTcc                                      <  4:210237/124‑1 (MQ=255)                                                                                                                                   cATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCt             >  4:289457/1‑149 (MQ=255)                                                                                                                                    aTAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTCCTCAACCtt            >  2:145581/1‑149 (MQ=255)                                                                                                                                     tAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTc           >  6:132452/1‑149 (MQ=255)                                                                                                                                      aGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTcc          >  7:58053/1‑149 (MQ=255)                                                                                                                                        gCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGc        >  1:189032/1‑149 (MQ=255)                                                                                                                                           aaaaaaGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTc     >  1:102811/1‑149 (MQ=255)                                                                                                                                           aaaaaaGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTc     <  1:199273/149‑1 (MQ=255)                                                                                                                                           aaaaaaGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTc     >  2:198294/1‑149 (MQ=255)                                                                                                                                           aaaaaaGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTc     <  7:83007/149‑1 (MQ=255)                                                                                                                                               aaGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGcc  >  2:8293/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                GCCAACATCTGCTTGGCCAGACGCCCACCCACATCACCACAACCCACCACCATCACGGAAAGGTCTGACATCTCTAAACTCCCTAGACCAATTGATCAGCTTGGCCACGCTACACGATGGGCGGCTAGACCATAGGCGAAAAAAGCTACAATATTACTTTTGTTAACAAGAATTACTTGCAATAATGACCGCCAAATTGTTCTCGGCCTACATCGAGGCCTTG‑GAGAACGCTCACTTTTTTCTTCATCAGGTCCGGCCAGCATGACACGTACTCAACCTTCCGCTTCGCC  >  NC_002947/6053396‑6053685

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