Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,914,848	G→C	7.9%	intergenic (‑62/‑184)	PP_3454 ← / → PP_3455	transcriptional regulator/multidrug RND transporter membrane fusion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,914,848	0	G	C	7.9%	200.2 / 6.5	63	intergenic (‑62/‑184)	PP_3454/PP_3455	transcriptional regulator/multidrug RND transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base C (2/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAG  >  NC_002947/3914700‑3914989                                                                                                                                                      |                                                                                                                                              cTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCgctg                                                                                                                                               >  4:328672/1‑149 (MQ=255)       aTGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAg                                                                                                                                         >  5:247713/1‑149 (MQ=255)       aTGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAg                                                                                                                                         <  5:322943/149‑1 (MQ=255)       aTGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAg                                                                                                                                         >  2:73262/1‑149 (MQ=255)         ggAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGt                                                                                                                                       >  5:31523/1‑149 (MQ=255)         ggAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGt                                                                                                                                       >  1:71530/1‑149 (MQ=255)                         cAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTGAACCCGCTCCTACGGGGACCTGCAGCAGCggg                                                                                                                        <  1:33343/148‑1 (MQ=37)                         cAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTGAACCCGCTCCTACGGGGACCTGCAGCAGCggg                                                                                                                        <  7:163846/148‑1 (MQ=37)                         cAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCggg                                                                                                                        <  2:71530/148‑1 (MQ=255)                          aGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGtt                                                                                                                      >  8:298378/1‑149 (MQ=255)                          aGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGt                                                                                                                       >  3:42065/1‑148 (MQ=255)                           ggTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCgg                                                                                                                         <  1:181096/145‑1 (MQ=255)                           ggTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCgg                                                                                                                         >  2:181096/1‑145 (MQ=255)                              aGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCAcc                                                                                                                  >  7:132979/1‑149 (MQ=255)                                cTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCAccc                                                                                                                 <  1:73262/149‑1 (MQ=255)                                           cTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCa                                                                                                     <  6:247713/149‑1 (MQ=255)                                                     gcgcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGccc                                                                                           <  7:298378/149‑1 (MQ=255)                                                       gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                          <  2:342848/132‑1 (MQ=255)                                                       gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                          >  7:146378/1‑132 (MQ=255)                                                       gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                          >  1:342848/1‑132 (MQ=255)                                                       gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                          <  8:146378/132‑1 (MQ=255)                                                       gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCAt                                                                                         >  8:83359/1‑149 (MQ=255)                                                             gtcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCAcc                                                                                                                  <  2:201304/118‑1 (MQ=255)                                                             gtcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCAcc                                                                                                                  >  1:201304/1‑118 (MQ=255)                                                              tcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAat                                                                                  <  4:42065/149‑1 (MQ=255)                                                               cTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                          <  7:27909/124‑1 (MQ=255)                                                               cTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                          >  8:27909/1‑124 (MQ=255)                                                                ttCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAata                                                                                 >  7:150386/1‑148 (MQ=255)                                                                  cGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAatat                                                                                <  3:91773/147‑1 (MQ=255)                                                                  cGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAata                                                                                 <  8:150386/146‑1 (MQ=255)                                                                    aCCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGaa                                                                            >  5:60432/1‑149 (MQ=255)                                                                      cAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAcc                                                                          >  4:143026/1‑149 (MQ=255)                                                                      cAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAcc                                                                          >  1:217690/1‑149 (MQ=255)                                                                       aGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAcc                                                                          >  5:96247/1‑148 (MQ=255)                                                                         gAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCaa                                                                       >  5:328371/1‑149 (MQ=255)                                                                          aaaaTGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAAc                                                                      <  4:48871/149‑1 (MQ=255)                                                                                   ggCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACa                                                             <  3:328672/149‑1 (MQ=255)                                                                                        aGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAAtata                                                                               >  5:118578/1‑126 (MQ=255)                                                                                        aGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAAtata                                                                               <  6:118578/126‑1 (MQ=255)                                                                                          gTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAt                                                      >  2:309340/1‑149 (MQ=255)                                                                                                gggCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAccc                                                                         <  6:312541/124‑1 (MQ=255)                                                                                                gggCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAccc                                                                         >  5:312541/1‑124 (MQ=255)                                                                                                gggCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTca                                                <  8:132979/149‑1 (MQ=255)                                                                                                  gCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTc                                                                 <  4:250908/130‑1 (MQ=255)                                                                                                  gCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTc                                                                 >  3:250908/1‑130 (MQ=255)                                                                                                       tGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACtttt                                         >  7:164249/1‑149 (MQ=255)                                                                                                         ggTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                        >  5:53579/1‑148 (MQ=255)                                                                                                         ggTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACtttt                                         <  6:328371/147‑1 (MQ=255)                                                                                                         ggTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTc                                       <  8:164249/149‑1 (MQ=255)                                                                                                             gAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcaca                                    <  6:53579/148‑1 (MQ=255)                                                                                                             gAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcaca                                    <  1:309340/148‑1 (MQ=255)                                                                                                              aaGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACtt                                           <  8:47670/140‑1 (MQ=255)                                                                                                              aaGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACtt                                           >  7:47670/1‑140 (MQ=255)                                                                                                                  aaGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGcc                              >  7:16967/1‑149 (MQ=255)                                                                                                                   aGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGcc                              <  2:217690/148‑1 (MQ=255)                                                                                                                           gCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGt                      <  3:143026/148‑1 (MQ=255)                                                                                                                             gAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTAt                    >  6:161059/1‑148 (MQ=255)                                                                                                                                cgcTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAAc                <  6:60432/149‑1 (MQ=255)                                                                                                                                 gcTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACa               <  6:96247/149‑1 (MQ=255)                                                                                                                                     cg‑cgGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAgccg           >  6:17968/1‑149 (MQ=255)                                                                                                                                           gtgaaCCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                        <  6:68830/112‑1 (MQ=21)                                                                                                                                           gtgaaCCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                        >  5:68830/4‑115 (MQ=21)                                                                                                                                           gtgaaCCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAGTCTTCACACTTTTACACACTCGCCCTACAGGTATTAAAAGCCGCCCAg      >  8:243953/4‑149 (MQ=25)                                                                                                                                               aCCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAg  >  1:172522/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                              CTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAG  >  NC_002947/3914700‑3914989

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