Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,503,142	C→G	9.2%	A31G (GCC→GGC)	oprC →	copper receptor OprC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,503,142	0	C	G	9.2%	188.9 / 10.0	65	A31G (GCC→GGC)	oprC	copper receptor OprC
Reads supporting (aligned to +/- strand):  ref base C (25/33);  new base G (2/4);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

CACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTAC  >  NC_002947/5502999‑5503290                                                                                                                                                 |                                                                                                                                                      cACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGcc                                                                                                                                                 <  3:158191/149‑1 (MQ=255)  aCGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCCCCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGccc                                                                                                                                                >  8:153526/1‑149 (MQ=255)      aGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGc                                                                                                                                            <  3:225896/149‑1 (MQ=255)       gCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGcc                                                                                                                                           <  8:327853/149‑1 (MQ=255)       gCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCAGGGCCGG‑CGGccagggc                                                                                                                                            >  7:57997/1‑143 (MQ=255)          aggaggCTTCGCGCTAACTCATGCCTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCAGGGCCGG‑CGGCCAGGGCCATg                                                                                                                                        <  4:207927/149‑1 (MQ=255)          aggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTTTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATg                                                                                                                                        <  1:263030/149‑1 (MQ=255)          aggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATg                                                                                                                                        <  8:175472/149‑1 (MQ=255)           ggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGa                                                                                                                                       <  1:125995/149‑1 (MQ=255)           ggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCAGGGCCGG‑CGGCCAGGGCCATgg                                                                                                                                       <  8:57997/149‑2 (MQ=255)            gaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGaa                                                                                                                                      <  5:34826/149‑1 (MQ=255)                 ttCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCa                                                                                                                                 >  1:218107/1‑149 (MQ=255)                   cgcgCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCTCAT‑GGCCAGGGCCGG‑CGGCCAGGGCCATgggcgcga                                                                                                                                 <  6:169438/147‑8 (MQ=38)                     cgcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCAtgt                                                                                                                              <  1:53001/149‑1 (MQ=255)                     cgcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCAtgtg                                                                                                                             >  8:91360/1‑149 (MQ=255)                      gcTAACTCATGACTGATCGACTGGAATCATTATGTACGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCAtgtg                                                                                                                             >  6:369153/1‑148 (MQ=255)                      gcTAACTCATGACTGATCGACT‑GAATAATTATGTACCGGGGCAACCCGGTTTTTGCCATTTGCCATCGTGGGAACATACCCGCCCATTTCCGCAAGCATCACGCCCCCAA‑GGCCCTTGCCGC‑CGGCCCAGGCCAATAAAACCAAGTGCa                                                                                                                           >  6:3021/1‑149 (MQ=255)                       cTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCa                                                                                                                           <  7:153526/149‑1 (MQ=255)                           cTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCacgac                                                                                                                       <  7:295526/149‑1 (MQ=255)                            tCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCAGGGCCGG‑CGGCCAGGGCCATg                                                                                                                                        <  6:320308/131‑1 (MQ=255)                            tCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCAGGGCCGG‑CGGCCAGGGCCATg                                                                                                                                        >  5:320308/1‑131 (MQ=255)                                 aCTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCAcc                                                                                                                  <  6:71793/148‑1 (MQ=255)                                        gACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGa                                                                                                          <  5:369153/149‑1 (MQ=255)                                        gACTGGAAACATTATGTCCCCCTGAACCCCTGTTTTTCACTAGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGA‑CGGCCCAGGCCAAGAACACCATGTACACGAAGCACCCCAGATGa                                                                                                          >  2:34493/1‑149 (MQ=255)                                         aCTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAg                                                                                                         <  7:207879/149‑1 (MQ=255)                                          cTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGc                                                                                                        <  7:91360/149‑1 (MQ=255)                                          cTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGc                                                                                                        >  1:233507/1‑149 (MQ=255)                                             gAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGccc                                                                                                      >  6:205154/1‑148 (MQ=255)                                                tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGccgc‑cg                                                                                                                                                    <  6:104720/99‑1 (MQ=255)                                                tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCACTGGccgc‑cg                                                                                                                                                    >  5:104720/1‑99 (MQ=255)                                                   ttATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCacga                                                                                                                        >  5:257590/1‑124 (MQ=255)                                                   ttATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCacga                                                                                                                        <  6:257590/124‑1 (MQ=255)                                                     aTGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCTTGCGTGGGACCATCGCCGCACTGTGCGGCTCGCTGCTCGCGCCCTT‑GGCCCTGGCCGC‑CGTCCCTGTCCATGAAGTGCATGTTCACGAGTGACGCGTGTTGAGCCCGACGGTGa                                                                                              <  1:34493/148‑1 (MQ=255)                                                       gTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCa                                                                                           >  5:102817/1‑149 (MQ=255)                                                                cACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTgg                                                                                  >  8:241047/1‑149 (MQ=255)                                                                           ttGCCTTGTCCCTGCGGGGGACCATCGCGGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCGGGCCATGAAGCCCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTc                                                                        <  8:284938/148‑1 (MQ=255)                                                                             gCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTcgc                                                                      <  6:102817/148‑1 (MQ=255)                                                                               cTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGCCGCACCCGAGTTGAGCCCGTCGGTGATCACCGCCGTGGCCCCGAGTTcgccgc                                                                   >  4:273580/1‑149 (MQ=255)                                                                               cTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGt                                                                                                             >  8:233353/1‑107 (MQ=255)                                                                               cTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGt                                                                                                             <  7:233353/107‑1 (MQ=255)                                                                                    cccTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACt                                                              >  2:284234/1‑149 (MQ=255)                                                                                     ccTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGtt                                                                                                            >  1:279898/1‑102 (MQ=255)                                                                                     ccTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGtt                                                                                                            <  2:279898/102‑1 (MQ=255)                                                                                     ccTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACtg                                                             >  8:166321/1‑149 (MQ=255)                                                                                     ccTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACtg                                                             >  1:235820/1‑149 (MQ=255)                                                                                         cGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACtgtg                                                           <  2:233507/147‑1 (MQ=255)                                                                                          gTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTc                                                        <  5:3021/149‑1 (MQ=255)                                                                                             ggACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCAc                                                      <  7:241047/148‑1 (MQ=255)                                                                                               aCCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCaa                                                   >  8:275647/1‑149 (MQ=255)                                                                                                 cATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAAcc                                                 <  7:166321/149‑1 (MQ=255)                                                                                                  aTCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGCCGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAAccc                                                >  3:188767/1‑149 (MQ=255)                                                                                                      ccgccCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGaaa                                            <  5:69678/149‑1 (MQ=255)                                                                                                      ccgccCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGaaa                                            <  2:235820/149‑1 (MQ=255)                                                                                                         cccTGTGCGGCTGGCTGCTCGCGGCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGACGACCTTGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCTCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGAc                                         <  4:191176/149‑1 (MQ=255)                                                                                                           cTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGAccc                                       <  7:275647/149‑1 (MQ=255)                                                                                                            tgtgCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCcg                                      >  5:359251/1‑149 (MQ=255)                                                                                                                cGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGCCGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCgcc                                  >  5:196396/1‑149 (MQ=255)                                                                                                                   cTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCgccagc                               >  5:55534/1‑149 (MQ=255)                                                                                                                    tCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCgccagc                               >  7:209051/1‑148 (MQ=255)                                                                                                                       ctgctCGCGCCCAT‑GGCCCTGGCCGCGCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGccagccag                            >  1:185625/1‑149 (MQ=255)                                                                                                                           tCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGt                           >  1:328469/1‑145 (MQ=255)                                                                                                                           tCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGt                           <  2:328469/145‑1 (MQ=255)                                                                                                                             gcgcCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGc                     <  7:338154/149‑1 (MQ=255)                                                                                                                             gcgcCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGc                     <  6:196396/149‑1 (MQ=255)                                                                                                                             gcgcCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGc                     >  3:94732/1‑149 (MQ=255)                                                                                                                                        cccTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGgc           <  6:55534/149‑1 (MQ=255)                                                                                                                                         ccTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGgcg          <  4:94732/149‑1 (MQ=255)                                                                                                                                               cgc‑cgGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGAc     >  3:47780/1‑148 (MQ=255)                                                                                                                                                 c‑cgGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTAc  >  4:177522/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      CACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAT‑GGCCCTGGCCGC‑CGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTAC  >  NC_002947/5502999‑5503290

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