Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,882,254	A→G	8.6%	T305A (ACC→GCC)	oprN →	multidrug RND transporter outer membrane protein OprN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,882,254	0	A	G	8.6%	153.4 / 4.6	59	T305A (ACC→GCC)	oprN	multidrug RND transporter outer membrane protein OprN
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base G (3/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

CTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCTGGGCAGCGTGCGG  >  NC_002947/3882107‑3882400                                                                                                                                                       |                                                                                                                                                   cTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACCTTGGC‑GTGGCCAc                                                                                                                                                   >  5:53001/1‑149 (MQ=255)       gTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCAGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCCAGCGCCAGCTGGCGCCGGCCACCGCCAACCTTGGC‑GTGGCCACCGCCGa                                                                                                                                             >  1:259580/1‑149 (MQ=255)           aCCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTg                                                                                                                                         >  4:210098/1‑149 (MQ=255)             cTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGAGGTGCGCCGCCGCCCGGACGTCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTg                                                                                                                                         >  1:107464/1‑149 (MQ=255)                  gCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCcgc                                                                                                                                  <  6:276746/149‑1 (MQ=255)                   cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GAGGCCACCGCCGACCTGTTCCcgcg                                                                                                                                 >  5:56414/1‑149 (MQ=255)                    ccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCcgcgc                                                                                                                                <  5:222621/149‑1 (MQ=255)                        ggCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCa                                                                                                                            <  5:282522/149‑1 (MQ=255)                                cGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGt                                                                                                                    >  4:30196/1‑149 (MQ=255)                                cGGCAATCGCCAATGCGCTGCCGGTCGGTGACCCGGTTTAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGACAACGTTTGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGt                                                                                                                    <  4:24079/149‑1 (MQ=255)                                  gCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTgg                                                                                                                  <  3:30196/149‑1 (MQ=255)                                        gCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCcgcgc                                                                                                                                <  1:244033/129‑1 (MQ=255)                                        gCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCcgcgc                                                                                                                                >  2:244033/1‑129 (MQ=255)                                         ccAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTc                                                                                                           >  2:288152/1‑149 (MQ=255)                                          cAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCAccgc                                                                                                                                                >  4:20970/1‑111 (MQ=255)                                          cAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCAccgc                                                                                                                                                <  3:20970/111‑1 (MQ=255)                                                gcTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTa                                                                                                    >  2:249509/1‑149 (MQ=255)                                                  tGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTAcc                                                                                                  >  8:21655/1‑149 (MQ=255)                                                    ccGGTCGGTGACGCGGGTGAGCT‑‑GCTGCGCCGCGGCCGGGACATCCGCAGCGGCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGAGCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCGTCGGCTTTAccg                                                                                                 <  2:259580/148‑1 (MQ=255)                                                    ccGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTAccgc                                                                                                >  2:113268/1‑149 (MQ=255)                                                      ggtcggtGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTAccgccg                                                                                              >  3:255713/1‑149 (MQ=255)                                                           gtGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGccgcccgc                                                                                         <  2:110167/149‑1 (MQ=255)                                                            tGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCcgcg                                                                                        <  2:165361/149‑1 (MQ=255)                                                               cccGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCt                                                                                     >  8:273926/1‑149 (MQ=255)                                                                   ggTGAGCT‑‑GCTGCGCCGCCGCCCGGACGTCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCa                                                                                 <  2:107464/149‑1 (MQ=255)                                                                      gAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGTGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGCCCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCg                                                                                   >  5:44183/1‑144 (MQ=255)                                                                        gct‑‑gctGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCg                                                                            <  7:27960/149‑1 (MQ=255)                                                                          t‑‑gctgcGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGc                                                                          >  5:337614/1‑149 (MQ=255)                                                                                 cgccgccgcCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGctct                                                                       <  6:337614/147‑1 (MQ=255)                                                                                 cgccgccgcCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTc                                                                     >  7:241172/1‑149 (MQ=255)                                                                                   ccgccgccCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATc                                                                             <  5:71080/139‑1 (MQ=255)                                                                                   ccgccgccCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATc                                                                             >  6:71080/1‑139 (MQ=255)                                                                                         cccGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGcc                                                              >  1:155065/1‑148 (MQ=255)                                                                                            ggACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACg                                                          >  6:217265/1‑149 (MQ=255)                                                                                                     cAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATc                                                                             >  4:123916/1‑121 (MQ=255)                                                                                                     cAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATc                                                                             <  3:123916/121‑1 (MQ=255)                                                                                                          ccGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGc                                            >  3:13255/1‑149 (MQ=255)                                                                                                           cGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGcc                                           <  3:210098/149‑1 (MQ=255)                                                                                                              gcgcCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGa                                        <  1:129782/149‑1 (MQ=255)                                                                                                              gcgcCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGa                                        <  4:13255/149‑1 (MQ=255)                                                                                                                   aGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCAtt                                   <  7:21655/149‑1 (MQ=255)                                                                                                                    gCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTa                                  <  1:288152/149‑1 (MQ=255)                                                                                                                       ggcggtggcCA‑CGCCAACGTTGGCGGTGGCCGCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCt                               >  4:39557/1‑149 (MQ=39)                                                                                                                        gcggtggcCA‑CGCCAACGTTGGCGGTGGCCGCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTg                              >  7:299189/1‑149 (MQ=39)                                                                                                                          ggcggcCACCGCCAACGTTGGC‑CTGGCCACCGCCGACCTGTTCCCCCGCGTAAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGCCTCGCAGATATGCTCTTCGGCCGCCAACTCCTCGGCACTGGGCCAGAGCATTACCTggt                            >  2:132258/1‑148 (MQ=255)                                                                                                                              gcCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGggcggc                        <  6:53001/149‑1 (MQ=255)                                                                                                                               ccACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGcc                       >  1:324267/1‑149 (MQ=255)                                                                                                                                cACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCATCGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCt                      >  4:287248/1‑149 (MQ=255)                                                                                                                                 aCCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGTCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCtt                     <  1:249509/149‑1 (MQ=255)                                                                                                                                         cGTTGGCGGTGGCCGCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTc                    <  8:232504/143‑1 (MQ=255)                                                                                                                                         agTTGGCGGTGGCCGCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTc                    >  7:232504/2‑143 (MQ=255)                                                                                                                                             ggcggTGGCCGCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCTGGGc           <  3:39557/144‑1 (MQ=255)                                                                                                                                             ggC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCTGGGCa          <  6:56414/148‑1 (MQ=255)                                                                                                                                             ggC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCTGGGCa          <  7:273926/148‑1 (MQ=255)                                                                                                                                                   tgCCTCCGCCGACGTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACGGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGTCCCGAGCATTACGTGTGCGGCCTTCGACCTGGGCAGCGTGc    <  3:182644/148‑1 (MQ=255)                                                                                                                                                   ggCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCTGGGCAGCGTGc    >  6:166024/1‑149 (MQ=255)                                                                                                                                                    gCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCTGGGCAGCGTGCg   <  2:324267/149‑1 (MQ=255)                                                                                                                                                     ccACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCTGGGCAGCGTGCgg  >  4:187282/1‑149 (MQ=255)                                                                                                                                                     ccACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCTGGGCAGCGTGCgg  <  5:217265/149‑1 (MQ=255)                                                                                                                                                       aCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGc                                            >  2:321914/1‑105 (MQ=255)                                                                                                                                                       aCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGc                                            <  1:321914/105‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                   CTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCT‑‑GCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCTGGGCAGCGTGCGG  >  NC_002947/3882107‑3882400

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