Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,503,190	C→G	11.6%	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,190	0	C	G	11.6%	174.9 / 13.9	61	A31G (GCC→GGC)	oprC	copper receptor OprC
Reads supporting (aligned to +/- strand):  ref base C (24/29);  new base G (5/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

TTGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGA  >  NC_002947/5503042‑5503334                                                                                                                                                     |                                                                                                                                                 ttGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGgccgcc                                                                                                                                                 <  1:1044579/150‑1 (MQ=255)          aaGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGc                                                                                                                                        <  1:796975/150‑1 (MQ=255)          aaGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGc                                                                                                                                        <  2:1191204/150‑1 (MQ=255)             ccAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCAt                                                                                                                                     >  1:744557/1‑150 (MQ=255)               aggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATgg                                                                                                                                   <  1:703953/150‑2 (MQ=255)               aggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATgg                                                                                                                                   <  1:7859/150‑2 (MQ=255)                  aggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGa                                                                                                                                <  2:194313/150‑1 (MQ=255)                  aggCTTCCCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGa                                                                                                                                <  2:194469/150‑1 (MQ=255)                    gCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGAcc                                                                                                                              >  2:793412/1‑150 (MQ=255)                      ttCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCGTGAGGAccgt                                                                                                                            <  2:744557/150‑3 (MQ=255)                      ttCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAt                                                                                                                            >  1:546394/1‑150 (MQ=255)                                       cTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGa                                                                                                           <  1:987863/150‑1 (MQ=255)                                        tGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAg                                                                                                          >  2:563008/1‑150 (MQ=255)                                            cGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGa                                                                                                      >  1:1005956/1‑150 (MQ=255)                                               cTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGcc                                                                                                   <  2:497743/150‑1 (MQ=255)                                                tGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGccc                                                                                                  <  2:646091/150‑1 (MQ=255)                                                  gAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTTGCCGCCGGCCCTGGCCATGAGGACCATGTGCACGACGCACCCGAGTTGAGCCCGa                                                                                                >  1:917950/1‑150 (MQ=255)                                                   aaTCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGAc                                                                                               >  1:180557/1‑150 (MQ=255)                                                       aTTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTg                                                                                           <  2:1106208/150‑1 (MQ=255)                                                        ttATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGa                                                                                          >  1:457913/1‑150 (MQ=255)                                                            gTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATc                                                                                        >  1:78737/1‑148 (MQ=255)                                                            gTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATc                                                                                        <  2:78737/148‑1 (MQ=255)                                                             tCCGGCTGCACCCCTGTTTTTGCGTTGTCTATGCGTGGGACCATCGCCGCCATGTGCGGCTCGCTGTTCGCGCGCATGGCTCTGGCAGACGGCCCTGGCCATGAATACCATGTGAAGGACGGAGACGAGTTGAGGGCGTAGGTGATCaaa                                                                                     >  2:467070/1‑148 (MQ=255)                                                                     cACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGc                                                                             <  1:563008/150‑1 (MQ=255)                                                                        cccTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGcccc                                                                          >  1:690252/1‑150 (MQ=255)                                                                          cTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCGGCTGGAGATCACACTGCTGAACTCCAGCACGGCTCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGa                                                                        >  1:1047543/1‑150 (MQ=255)                                                                          cTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGa                                                                        <  2:690252/150‑1 (MQ=255)                                                                                 tGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTcgcc                                                                 >  2:230225/1‑150 (MQ=255)                                                                                   ccTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTcgccgc                                                               >  2:923574/1‑150 (MQ=255)                                                                                     ttGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGCGTTCGCCGCTg                                                             >  2:219369/1‑150 (MQ=255)                                                                                          ccTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACtgt                                                        >  2:981682/1‑150 (MQ=255)                                                                                              cGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTc                                                    <  1:793412/150‑1 (MQ=255)                                                                                               gTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCa                                                   <  2:180557/150‑1 (MQ=255)                                                                                                tGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCAc                                                  <  1:260171/150‑1 (MQ=255)                                                                                                tGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCAc                                                  <  2:709808/150‑1 (MQ=255)                                                                                                       aTCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTgg                                                      <  2:1064937/139‑1 (MQ=255)                                                                                                       aTCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTgg                                                      >  1:1064937/1‑139 (MQ=255)                                                                                                        tCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACt                                                          <  2:82913/134‑1 (MQ=255)                                                                                                        tCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACt                                                          >  1:82913/1‑134 (MQ=255)                                                                                                        tCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGa                                          >  1:842130/1‑150 (MQ=255)                                                                                                        tCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGAAGGCGTGTTTATCCCGACGGTGATCACCGCGGTGGGCCCGAGTTCGCCGCTGACTGTGGGCACCAACCCGa                                          >  1:811605/1‑150 (MQ=255)                                                                                                               ccTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGAccc                                   <  2:1005956/150‑1 (MQ=255)                                                                                                                 tgtgCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACtgtg                                                       >  2:507747/1‑128 (MQ=255)                                                                                                                 tgtgCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACtgtg                                                       <  1:507747/128‑1 (MQ=255)                                                                                                                           gctgctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCa                                                                                                                >  1:740047/1‑61 (MQ=255)                                                                                                                           gctgctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCa                                                                                                                <  2:740047/61‑1 (MQ=255)                                                                                                                           gctgctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGt                                                     <  2:273865/120‑1 (MQ=255)                                                                                                                           gctgctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGt                                                     >  1:273865/1‑120 (MQ=255)                                                                                                                                 cgcgCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGc                 <  2:943016/150‑1 (MQ=255)                                                                                                                                   tcttcATGGCCAGGGCCGGCGGCCAGGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCa               >  1:89760/5‑150 (MQ=255)                                                                                                                                   cgcCCATGGCCCTGGACGCCCGACCTGGCCAAGAAGCTCATGCGCACTTCGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAGCCCGAAAGACCCGCGCCAGCCAGTGCCGGCCa               <  2:917950/150‑1 (MQ=255)                                                                                                                                      ccATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGCTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCg            >  2:780812/1‑150 (MQ=255)                                                                                                                                       cATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGaaa                                        <  2:753714/121‑1 (MQ=255)                                                                                                                                       cATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGaaa                                        >  1:753714/1‑121 (MQ=255)                                                                                                                                       cATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGa           <  2:759318/150‑1 (MQ=255)                                                                                                                                         tggccagggCCGGCGGCCAGGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACg         >  2:306434/8‑150 (MQ=255)                                                                                                                                         tGGCCAGGGCCGGCGGCCAGGGCCATGGAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACg         >  2:195263/1‑150 (MQ=255)                                                                                                                                         tGGCCAGGGCCGGCGGCCAGGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTgg                                                                              >  2:1143090/1‑81 (MQ=37)                                                                                                                                           gCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGgc       <  2:1047543/150‑1 (MQ=255)                                                                                                                                            ccagggCCGGCGGCCAGGGCCATGGAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGgcg      >  2:344947/5‑150 (MQ=255)                                                                                                                                             ccTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGgcgc     <  1:923574/150‑1 (MQ=255)                                                                                                                                                ggCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGa  <  1:230225/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TTGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGA  >  NC_002947/5503042‑5503334

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