Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,916,544	C→G	5.9%	A94A (GCC→GCG)	mexB →	multidrug resistance protein MexB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,916,544	0	C	G	5.9%	217.4 / 4.5	68	A94A (GCC→GCG)	mexB	multidrug resistance protein MexB
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base G (2/2);  total (34/34)
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.32e-01

TACGCCGTGTACCCGGGTGCCTCGGCGGCGACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCTTCCTGCTGTTGGCGACCCTCACCTCC  >  NC_002947/3916398‑3916691                                                                                                                                                   |                                                                                                                                                    tACGCCGTGTACCCGGGTGCCTCGGCGGCGACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCg                                                                                                                                                  >  1:159426/1‑150 (MQ=255)  aCGCCGTGTACCCGGGTGCCTCGGCGGCGACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCgg                                                                                                                                                 <  1:1269764/150‑1 (MQ=255)        tgtACCCGGGTGCCTCGGCGGCGACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTaccc                                                                                                                                           <  1:1032218/150‑1 (MQ=255)          tACCCGGGTGCCTCGGCGGCGACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTacccac                                                                                                                                         <  2:693791/150‑1 (MQ=255)             ccGGGTGCCTCGGCGGCGACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCg                                                                                                                                      >  1:111003/1‑150 (MQ=255)               gggTGCCTCGGCGGCGACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCATGGGCAGCGCCACCATCACCGACAACTTCGCCCCGGTTACACACCCggc                                                                                                                                    >  2:620054/1‑149 (MQ=255)                    ccTCGGCGGCGACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTgg                                                                                                                               <  2:969250/150‑1 (MQ=255)                    ccTCGGCGGCGACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTgg                                                                                                                               <  2:1185741/150‑1 (MQ=255)                     cTCGGCGGCGACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGc                                                                                                                              >  2:859933/1‑150 (MQ=255)                         gcggcgACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAg                                                                                                                          <  1:414157/150‑1 (MQ=255)                            gcgACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTc                                                                                                                       >  1:1357602/1‑150 (MQ=255)                              gACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGa                                                                                                                     >  1:971816/1‑150 (MQ=255)                                   tGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGc                                                                                                                >  1:357692/1‑150 (MQ=255)                                        gAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAAc                                                                                                           >  2:136966/1‑150 (MQ=255)                                                 gTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGaaa                                                                                                  >  2:474019/1‑150 (MQ=255)                                                 gTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGaaa                                                                                                  <  1:529379/150‑1 (MQ=255)                                                 gTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGaaa                                                                                                  <  2:1239277/150‑1 (MQ=255)                                                    aGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTg                                                                                               >  2:776812/1‑150 (MQ=255)                                                      ccTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGt                                                                                             <  2:1357602/150‑1 (MQ=255)                                                                 aGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGaa                                                                                                            <  1:8182/124‑1 (MQ=255)                                                                 aGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGaa                                                                                                            >  2:8182/1‑124 (MQ=255)                                                                    aggtCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGc                                                                               <  2:1109536/147‑1 (MQ=255)                                                                           cGGCGCCGACAACCTTCTGCACTTCCAGTCGCACATCAGCCTGGCCAGCGCCACCATCACCGCCACCTTCGCCCCGGTGACCCACCCGGACCGGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGccc                                                                        <  1:620054/150‑1 (MQ=255)                                                                                ccGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCa                                                                   <  1:474019/150‑1 (MQ=255)                                                                                  gACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCgg                                                                                                                                                 >  1:1111709/1‑70 (MQ=255)                                                                                  gACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCgg                                                                                                                                                 <  2:1111709/70‑1 (MQ=255)                                                                                  gACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTAGCGCCTGCCGCGCCCGGTCAcc                                                                 >  1:636911/1‑150 (MQ=255)                                                                                         tgctgTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAgg                                                          <  2:357692/150‑1 (MQ=255)                                                                                         tgctgTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAgg                                                          >  1:1160734/1‑150 (MQ=255)                                                                                           ctgTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGc                                                        >  2:852862/1‑150 (MQ=255)                                                                                                cTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGt                                                                                                                  <  2:1387943/87‑1 (MQ=255)                                                                                                cTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGt                                                                                                                  >  1:1387943/1‑87 (MQ=255)                                                                                                 ttCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCTCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAg                                                  <  2:952688/150‑1 (MQ=255)                                                                                                 ttCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGcccc                                                                                                                                                   <  2:175227/53‑1 (MQ=255)                                                                                                 ttCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGcccc                                                                                                                                                   >  1:175227/1‑53 (MQ=255)                                                                                                 ttCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTg                                                                                               >  2:335015/1‑105 (MQ=255)                                                                                                 ttCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTg                                                                                               <  1:335015/105‑1 (MQ=255)                                                                                                 ttCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAgg                                                          >  1:722899/1‑142 (MQ=255)                                                                                                 ttCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAgg                                                          <  2:722899/142‑1 (MQ=255)                                                                                                 ttCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCGGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAg                                                  >  1:810723/1‑150 (MQ=255)                                                                                                   cGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGt                                                >  1:911545/1‑150 (MQ=255)                                                                                                   cGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGt                                                <  2:1160734/150‑1 (MQ=255)                                                                                                    gAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTg                                               <  1:776812/150‑1 (MQ=255)                                                                                                    gAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTg                                               <  1:852862/150‑1 (MQ=255)                                                                                                          cAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGaaa                                         <  1:136966/150‑1 (MQ=255)                                                                                                           agagCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAg                                        >  1:283927/1‑150 (MQ=255)                                                                                                            gagCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCGCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGc                                                    <  2:346796/137‑1 (MQ=255)                                                                                                            gagCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCGCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGc                                                    >  1:346787/1‑137 (MQ=255)                                                                                                            gagCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCGCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGc                                                    >  1:346796/1‑137 (MQ=255)                                                                                                            gagCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCGCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGc                                                    <  2:346787/137‑1 (MQ=255)                                                                                                                 gcCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGc                                                        >  1:963622/1‑128 (MQ=255)                                                                                                                 gcCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGc                                                        <  2:963622/128‑1 (MQ=255)                                                                                                                 gcCTGGGCAGCGCCACCATCACCGCAACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGGGGAGAAAGTGTCCa                                  >  2:636228/1‑150 (MQ=255)                                                                                                                    tGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCg                               >  1:1142364/1‑150 (MQ=255)                                                                                                                            gcCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCcagc                                                             <  1:670553/112‑1 (MQ=255)                                                                                                                            gcCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCcagc                                                             >  2:670553/1‑112 (MQ=255)                                                                                                                             ccaccaTAACCGCCACCTTCGACCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCTTCctgc                      <  1:506252/150‑1 (MQ=255)                                                                                                                                ccaTCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCGGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGCTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCTTCCTGCTGt                   >  1:87298/1‑150 (MQ=255)                                                                                                                                 caTCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGa                                                                                          >  1:85863/1‑78 (MQ=255)                                                                                                                                 caTCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGa                                                                                          <  2:85863/78‑1 (MQ=255)                                                                                                                                 caTCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACGCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCTTCCTGCTGtt                  >  1:1224615/1‑150 (MQ=255)                                                                                                                                     aCCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCTTCCTGCTGTTGGCg              >  1:109198/1‑150 (MQ=255)                                                                                                                                      ccGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCtt                           >  2:618519/1‑136 (MQ=255)                                                                                                                                      ccGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCtt                           <  1:618519/136‑1 (MQ=255)                                                                                                                                        gCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCTTCCTGCTGTTGGCGAcc           >  2:1007735/1‑150 (MQ=255)                                                                                                                                         ccACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCTTCCTGCTGTTGGCGAccc          <  2:1120768/150‑1 (MQ=255)                                                                                                                                             cTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCGGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCTTCCTGCTGTTGGCGACCCTCAc      <  2:87298/150‑1 (MQ=255)                                                                                                                                             cTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCTTCCTGCTGTTGGCGACCCTCAc      >  1:1332663/1‑150 (MQ=255)                                                                                                                                             cTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCTTCCTGCTGTTGGCGACCCTCAc      <  2:1142364/150‑1 (MQ=255)                                                                                                                                                 gCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCTTCCTGCTGTTGGCGACCCTCACCTcc  <  1:859933/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TACGCCGTGTACCCGGGTGCCTCGGCGGCGACCATGGACGAAAGCGTGGTCAGCCTGATCGAACAGGAGCTCAACGGCGCCGACAACCTGCTGTACTTCGAGTCGCAGAGCAGCCTGGGCAGCGCCACCATCACCGCCACCTTCGCCCCGGGTACCCACCCGGACCTGGCCCAGGTCGACGTGCAGAACCGCCTGAAAGTGGTCGAGTCGCGCCTGCCGCGCCCGGTCACCCAGCAGGGCCTGCAGGTGGAGAAAGTGTCCACCGGCTTCCTGCTGTTGGCGACCCTCACCTCC  >  NC_002947/3916398‑3916691

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