Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,914,857	C→A	8.1%	intergenic (‑71/‑175)	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,857	0	C	A	8.1%	196.6 / 6.2	62	intergenic (‑71/‑175)	PP_3454/PP_3455	transcriptional regulator/multidrug RND transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base A (2/3);  total (31/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 = 9.03e-01

CAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCA  >  NC_002947/3914724‑3915003                                                                                                                                       |                                                                                                                                                   cAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTGAACCCGCTCCTACGGGGACCTGCAGCAGCggg                                                                                                                                      <  7:163846/148‑1 (MQ=37) cAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTGAACCCGCTCCTACGGGGACCTGCAGCAGCggg                                                                                                                                      <  1:33343/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                                                                                                                        <  8:146378/132‑1 (MQ=255)                               gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                        >  1:342848/1‑132 (MQ=255)                               gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                        >  7:146378/1‑132 (MQ=255)                               gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                        <  2:342848/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                                                                                                                        >  8:27909/1‑124 (MQ=255)                                       cTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                        <  7:27909/124‑1 (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                                                                                       >  5:312541/1‑124 (MQ=255)                                                                        gggCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAccc                                                                                       <  6:312541/124‑1 (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                                                      >  5:68830/4‑115 (MQ=21)                                                                                                                   gtgaaCCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                                      <  6:68830/112‑1 (MQ=21)                                                                                                                   gtgaaCCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAGTCTTCACACTTTTACACACTCGCCCTACAGGTATTAAAAGCCGCCCAg                    >  8:243953/4‑149 (MQ=25)                                                                                                                       aCCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAg                >  1:172522/1‑149 (MQ=255)                                                                                                                                  cAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCt       >  4:195366/1‑147 (MQ=255)                                                                                                                                  cAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCt       >  2:39097/1‑147 (MQ=255)                                                                                                                                  cAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTa     <  5:161059/149‑1 (MQ=255)                                                                                                                                     gTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCa  >  1:260594/1‑149 (MQ=255)                                                                                                                                       |                                                                                                                                                   CAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCG‑CGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCA  >  NC_002947/3914724‑3915003

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