Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,032,844	Δ1 bp	6.2%	intergenic (+52/+84)	dppA‑IV → / ← PP_5284	dipeptide transport protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,032,844	0	G	.	6.2%	195.7 / 6.8	65	intergenic (+52/+84)	dppA‑IV/PP_5284	dipeptide transport protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/34);  new base . (2/2);  total (29/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.28e-01

GACGGCGGCGACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATC  >  NC_002947/6032699‑6032993                                                                                                                                                  |                                                                                                                                                        gACGGCGGCGACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTg                                                                                                                                                     <  2:381881/150‑1 (MQ=255)  aCGGCGGCGACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTgg                                                                                                                                                    <  2:1285796/150‑1 (MQ=255)    ggcggcGACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCa                                                                                                                                                  <  2:89172/150‑1 (MQ=255)     gcggcgACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAACAAACTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑g                                                                                                                                                <  1:1502049/150‑1 (MQ=39)          gACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTNGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGAc                                                                                                                                          <  1:1075566/150‑1 (MQ=255)             cTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTGCGGGCGCTGAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCt                                                                                                                                       >  2:1269593/1‑150 (MQ=255)              ttGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTa                                                                                                                                      <  2:2368104/150‑1 (MQ=255)              ttGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTa                                                                                                                                      <  1:1379943/150‑1 (MQ=255)                 ctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTAccc                                                                                                                                   <  2:1532476/150‑1 (MQ=255)                       cAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCCACGGGGTGGCA‑GT‑TGACGCTACCCCTGTAg                                                                                                                             <  1:1269593/150‑1 (MQ=39)                        aaGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTCA‑ACTGCCA‑GGCTGGCGCTACCCCTGTAgg                                                                                                                            >  1:499969/1‑150 (MQ=18)                         aGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTGAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGc                                                                                                                                                        >  2:2173764/1‑123 (MQ=255)                         aGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTGAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGc                                                                                                                                                        <  1:2173764/123‑1 (MQ=255)                              tCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCGAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCggg                                                                                                                      >  2:2296276/1‑150 (MQ=255)                                gAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGtt                                                                                                                    <  1:1361523/150‑1 (MQ=255)                                     gTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAACAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTAccc                                                                                                               <  1:364102/150‑1 (MQ=255)                                        ccAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCA‑ACTGCCAGGC‑TGGCGCTACCCCTGTAGGAGCggg                                                                                                                      >  1:2059517/1‑140 (MQ=21)                                        ccAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCA‑ACTGCCA‑GGCTGGCGCTACCCCTGTAGGAGCggg                                                                                                                      <  2:2059517/140‑1 (MQ=21)                                           ggTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGgaa                                                                                                         <  1:2402309/150‑1 (MQ=255)                                           ggTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGaaa                                                                                                         <  1:1529045/150‑1 (MQ=255)                                            gTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAg                                                                                                        >  2:1187304/1‑150 (MQ=255)                                                 cccGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTGA‑ACTGCCAGGC‑TGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGaagaggccg                                                                                                   <  2:499969/150‑8 (MQ=11)                                                   cGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGt                                                                                                 >  2:772855/1‑150 (MQ=255)                                                               tGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCtata                                                                                     <  1:1187304/150‑1 (MQ=255)                                                                  aTTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAgg                                                                                  <  1:1803332/150‑1 (MQ=255)                                                                                  ggtgaaCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTc                                                                  <  2:1922790/150‑1 (MQ=255)                                                                                          aGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACcaca                                                          >  1:1578522/1‑150 (MQ=255)                                                                                           gTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTa                                                                                                                              <  1:301073/81‑1 (MQ=255)                                                                                           gTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTa                                                                                                                              >  2:301073/1‑81 (MQ=255)                                                                                             aaGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAgg                                                                                  <  1:450309/123‑1 (MQ=255)                                                                                             aaGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAgg                                                                                  >  2:450309/1‑123 (MQ=255)                                                                                             aaGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGa                                                               <  1:1531101/142‑1 (MQ=255)                                                                                             aaGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGATTTACCCGCGAAGGGGCCCGTACTGAACCTATAAGGCTCTCAGTACTCCCTCAGa                                                               >  2:1531101/1‑142 (MQ=255)                                                                                                 cTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTCCAGGGGTAGCACCACCCTGGCC‑GC‑CGACGCAATCAAAGAAACAACGAGTTTACCCACGAAAGGGCCCGTACTGACCCTAGAAAGCCCTCATTGCTCCCTCAGCGCACAAaccacc                                                   >  2:887274/1‑150 (MQ=16)                                                                                                  tGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCAccc                                                  >  1:1633003/1‑150 (MQ=255)                                                                                                       cGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGc                                                                                                                         <  2:100690/74‑1 (MQ=255)                                                                                                       cGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGc                                                                                                                         >  1:100690/1‑74 (MQ=255)                                                                                                        ggCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGaaa                                                                                                         >  2:528295/1‑89 (MQ=255)                                                                                                        ggCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGaaa                                                                                                         <  1:528295/89‑1 (MQ=255)                                                                                                         gCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCg                                           <  2:1633003/150‑1 (MQ=255)                                                                                                              ttCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                                                     >  1:2213227/1‑87 (MQ=255)                                                                                                              ttCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                                                     <  2:2213227/87‑1 (MQ=255)                                                                                                              ttCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTAc                                                                                               >  2:1223590/1‑93 (MQ=255)                                                                                                              ttCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTAc                                                                                               <  1:1223590/93‑1 (MQ=255)                                                                                                                     gTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAg                                                                                   >  1:1526397/1‑98 (MQ=255)                                                                                                                     gTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAg                                                                                   <  2:1526397/98‑1 (MQ=255)                                                                                                                      tAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACg                              >  2:2176455/1‑150 (MQ=255)                                                                                                                        aaCCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGCCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGc                            >  1:1365260/1‑150 (MQ=255)                                                                                                                          cccGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACACCGGCTc                          >  2:2140169/1‑150 (MQ=255)                                                                                                                          cccGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTc                          >  2:2221965/1‑150 (MQ=255)                                                                                                                           ccGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTc                                            >  2:733230/1‑131 (MQ=255)                                                                                                                           ccGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTc                                            <  1:733230/131‑1 (MQ=255)                                                                                                                               tCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTCCCCCTGTAGGAGCGTTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCAt                     <  1:772855/150‑1 (MQ=255)                                                                                                                                  tACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                  <  2:1578522/150‑1 (MQ=255)                                                                                                                                          tAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTACGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCAACCCATGCTCGACCAGCGACAACGGCAACCCATCACCAatgacg          >  2:858023/1‑148 (MQ=255)                                                                                                                                            gcgcCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCggg                                                                                                                      >  2:1905325/1‑40 (MQ=255)                                                                                                                                            gcgcCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCggg                                                                                                                      <  1:1905325/40‑1 (MQ=255)                                                                                                                                             cgcCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                  <  1:671018/139‑1 (MQ=255)                                                                                                                                             cgcCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                  >  2:671018/1‑139 (MQ=255)                                                                                                                                              gcCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCa                                                                            <  2:1034426/80‑1 (MQ=255)                                                                                                                                              gcCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCa                                                                            >  1:1034426/1‑80 (MQ=255)                                                                                                                                              gcCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTc                                                                  <  2:2118238/90‑1 (MQ=255)                                                                                                                                              gcCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTc                                                                  >  1:2118238/1‑90 (MQ=255)                                                                                                                                                cAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGtt                                                                                                                    <  2:33405/38‑1 (MQ=255)                                                                                                                                                cAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGtt                                                                                                                    >  1:33405/1‑38 (MQ=255)                                                                                                                                                  gCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATc  >  1:2266099/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                        GACGGCGGCGACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCA‑GT‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATC  >  NC_002947/6032699‑6032993

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑