Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,503,607	A→G	5.4%	N186S (AAC→AGC)	oprC →	copper receptor OprC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,503,607	0	A	G	5.4%	218.9 / 2.0	74	N186S (AAC→AGC)	oprC	copper receptor OprC
Reads supporting (aligned to +/- strand):  ref base A (35/35);  new base G (2/2);  total (37/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.94e-01

CTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGTCGCGCTGGGAAAAGTGGAACGGCGATGT  >  NC_002947/5503461‑5503754                                                                                                                                                     |                                                                                                                                                     cttcGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCACGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑Cg                                                                                                                                                   <  4:124671/146‑1 (MQ=255)     gACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCg                                                                                                                                               >  4:83368/1‑149 (MQ=255)           cTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGa                                                                                                                                         >  2:1961/1‑149 (MQ=255)           cTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGAACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGa                                                                                                                                         <  1:267571/149‑1 (MQ=255)              aCCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGAc                                                                                                                                        <  6:162979/147‑1 (MQ=255)                cGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAgg                                                                                                                                    >  4:223980/1‑149 (MQ=255)                      cAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGGCAACGCAAGCCTGCTGGCCG‑GATCCA‑A‑CGGCCGCTTCGACAAGGTGCTg                                                                                                                               >  4:113528/1‑148 (MQ=255)                        aaGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGAt                                                                                                                            >  7:38750/1‑149 (MQ=255)                              cGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTg                                                                                                                               <  5:234493/140‑1 (MQ=255)                              cGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTg                                                                                                                               >  6:234493/1‑140 (MQ=255)                              cGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCcgcg                                                                                                                      >  3:332635/1‑149 (MQ=255)                                             gggggCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTg                                                                                                                               >  4:282797/1‑125 (MQ=255)                                             gggggCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTg                                                                                                                               <  3:282797/125‑1 (MQ=255)                                              ggggCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAgg                                                                                                      >  2:334517/1‑149 (MQ=255)                                              ggggCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAg                                                                                                       <  2:271033/148‑1 (MQ=255)                                                ggCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGc                                                                                                    >  7:284190/1‑149 (MQ=255)                                                ggCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGc                                                                                                    <  2:22515/149‑1 (MQ=255)                                                   ccGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTAt                                                                                                 >  7:317432/1‑149 (MQ=255)                                                   ccGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTAt                                                                                                 >  2:55714/1‑149 (MQ=255)                                                   ccGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTAt                                                                                                 <  3:223980/149‑1 (MQ=255)                                                      gcggcTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCATGGTTCCGGTGATCTAGAGCGCCCGCATCGACAATGAGCTGGATGCCGCCGCCCGCAACACACAGGGCTAt                                                                                                 >  1:112157/1‑149 (MQ=255)                                                               ccggcACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGc                                                                                           >  6:63482/1‑143 (MQ=255)                                                               ccggcACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGc                                                                                           <  5:63482/143‑1 (MQ=255)                                                               ccggcACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTtcgtc                                                                                     <  2:200724/149‑1 (MQ=255)                                                                cggcACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTtcgtcg                                                                                    <  4:332635/149‑1 (MQ=255)                                                                 ggcACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCgg                                                                                   <  8:317432/149‑1 (MQ=255)                                                                       aTCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCa                                                                             <  3:162879/149‑1 (MQ=255)                                                                        tCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAg                                                                            >  2:136554/1‑149 (MQ=255)                                                                         ccTCTTCGAGCTCGTACCTGAGAAATTCGGCACCCTCGGCAGGCGCGTCATCGCAATCCTGCTGGCCG‑GCTCTA‑A‑CGGCGGCTTCGACATGGTGCTGGATGCCGcggc                                                                                                                    <  4:148557/108‑1 (MQ=255)                                                                           tctTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCa                                                                         >  8:58014/1‑149 (MQ=255)                                                                               cGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCAc                                                                        <  3:83368/146‑1 (MQ=255)                                                                                gAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCt                                                                     <  1:55714/148‑1 (MQ=255)                                                                                      gAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGccc                                                                                             >  3:207590/1‑118 (MQ=255)                                                                                      gAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGccc                                                                                             <  4:207590/118‑1 (MQ=255)                                                                                          cTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTAc                                                          >  6:233871/1‑149 (MQ=255)                                                                                           tGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTa                                                           <  7:97500/147‑1 (MQ=255)                                                                                           tGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTa                                                           >  8:97500/1‑147 (MQ=255)                                                                                            gagaAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGga                                                                 <  5:332296/140‑1 (MQ=255)                                                                                            gagaAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGga                                                                 >  6:332296/1‑140 (MQ=255)                                                                                               aaaTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGgac                                                                <  5:260760/138‑1 (MQ=255)                                                                                               aaaTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGgac                                                                >  6:260760/1‑138 (MQ=255)                                                                                               aaaTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCacga                                                     >  1:48855/1‑149 (MQ=255)                                                                                               aaaTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCacga                                                     >  6:16383/1‑149 (MQ=255)                                                                                                 aTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑G‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGga                                                                 >  4:153841/1‑135 (MQ=255)                                                                                                 aTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑G‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGga                                                                 >  4:148621/1‑135 (MQ=255)                                                                                                 aTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑G‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGga                                                                 <  3:153841/135‑1 (MQ=255)                                                                                                 aTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑G‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGga                                                                 <  3:148621/135‑1 (MQ=255)                                                                                                 aTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCacgacg                                                   <  8:284190/149‑1 (MQ=255)                                                                                                 aTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCacgacg                                                   <  1:334517/149‑1 (MQ=255)                                                                                                 aTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCacgacg                                                   <  5:233871/149‑1 (MQ=255)                                                                                                 aTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCacgac                                                    >  5:124068/1‑148 (MQ=255)                                                                                                 aTTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCacgac                                                    <  6:124068/148‑1 (MQ=255)                                                                                                     ggCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGGCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGcaa                                               <  8:309900/149‑1 (MQ=255)                                                                                                        aCCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGcaacaa                                            <  2:112157/149‑1 (MQ=255)                                                                                                          ccTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGcaaca                                             >  6:166397/1‑146 (MQ=255)                                                                                                          ccTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGcaaca                                             >  2:63595/1‑146 (MQ=255)                                                                                                            tCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGa                                         >  5:158471/1‑148 (MQ=255)                                                                                                              ggCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGacac                                      <  2:48855/149‑1 (MQ=255)                                                                                                               gCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACAcc                                     <  1:63595/149‑1 (MQ=255)                                                                                                                    cgcgTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGcc                                <  1:1961/149‑1 (MQ=255)                                                                                                                    cgcgTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGcc                                <  5:166397/149‑1 (MQ=255)                                                                                                                      cgTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGt                              >  7:183232/1‑149 (MQ=255)                                                                                                                      cgTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGt                              >  8:37402/1‑149 (MQ=255)                                                                                                                       gTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGACCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGTc                             >  7:190000/1‑149 (MQ=255)                                                                                                                         cAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTAACACGACGGCAACAAAGac                                        >  8:127440/1‑136 (MQ=255)                                                                                                                         cAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTAACACGACGGCAACAAAGac                                        <  7:127440/136‑1 (MQ=255)                                                                                                                            cGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGTCGCGCt                        >  2:33071/1‑149 (MQ=255)                                                                                                                                      ctgGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGTCGCGCTGGGAAAAGt               <  8:183232/148‑1 (MQ=255)                                                                                                                                      ctgGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGTCGCGCTGGGAAAAGTg              <  5:16383/149‑1 (MQ=255)                                                                                                                                         gCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGTCGCGCTGGGAAAAGTGGaa           <  7:58014/149‑1 (MQ=255)                                                                                                                                          ccG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGTCGCGCTGGGAAAAGTGGAAc          <  8:38750/149‑1 (MQ=255)                                                                                                                                              gCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGTCGCGCTGGGAAAAGTGGAACGGc       >  7:72211/1‑149 (MQ=255)                                                                                                                                                tCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGTCGCGCTgg                      >  8:226470/1‑132 (MQ=255)                                                                                                                                                  cA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGTCGCGCTGGGAAAAGTGGAACGGCGATg   >  4:268944/1‑149 (MQ=255)                                                                                                                                                   a‑a‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGTCGCGCTGGGAAAAGTGGAACGGCGATGt  >  3:333283/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCG‑GCTCCA‑A‑CGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCAGGGCTATGCCCGCTTCGTCGGCAACCAGTCACGCTCGGACGACTACCACGACGGCAACAAAGACACCGTGCCGTCGCGCTGGGAAAAGTGGAACGGCGATGT  >  NC_002947/5503461‑5503754

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