Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,421,952	G→C	5.9%	P80P (CCG→CCC) ‡	moeA →	molybdenum::molybdopterin ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,421,952	0	G	C	5.9%	218.9 / 3.4	68	P80P (CCG→CCC) ‡	moeA	molybdenum::molybdopterin ligase
Reads supporting (aligned to +/- strand):  ref base G (37/27);  new base C (2/2);  total (39/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 = 2.13e-01

ATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGATGGCC  >  NC_002947/2421806‑2422100                                                                                                                                                   |                                                                                                                                                       aTCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCt                                                                                                                                                      >  3:74786/1‑149 (MQ=255)    cgcgAGACCGAAACCCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                                    <  5:301630/148‑1 (MQ=255)     gcgAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCg                                                                                                                                                  >  7:16666/1‑149 (MQ=255)         gaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGT‑ca                                                                                                                                             <  4:22168/149‑1 (MQ=255)            cGAAAACCTGCCGCTGGCCGAAGCCGATGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑cagc                                                                                                                                           >  7:62276/1‑148 (MQ=255)                       cGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTAt                                                                                                                                   <  4:175095/145‑1 (MQ=255)                       cGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTAt                                                                                                                                   >  3:175095/1‑145 (MQ=255)                             ccGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTc                                                                                                                               >  5:302478/1‑143 (MQ=255)                             ccGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTc                                                                                                                               <  6:302478/143‑1 (MQ=255)                             ccGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGc                                                                                                                         >  8:239602/1‑149 (MQ=255)                             ccGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGc                                                                                                                         >  6:134050/1‑149 (MQ=255)                              cGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGc                                                                                                                         >  2:206928/1‑148 (MQ=255)                              cGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGc                                                                                                                         >  4:100675/1‑148 (MQ=255)                               gAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGT‑ca                                                                                                                                             <  1:250570/127‑1 (MQ=38)                               gAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGT‑ca                                                                                                                                             >  2:250570/1‑127 (MQ=38)                                       gggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                                    >  8:106845/1‑113 (MQ=255)                                       gggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                                    <  7:106845/113‑1 (MQ=255)                                            gcgTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCg                                                                                                          >  1:205246/1‑149 (MQ=255)                                              gTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                                        <  4:68539/149‑1 (MQ=255)                                              gTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                                        >  7:252231/1‑149 (MQ=255)                                                 cTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGc                                                                                                      >  6:305732/1‑148 (MQ=255)                                                    gCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑c                                                                                                                                              >  6:201883/1‑105 (MQ=255)                                                    gCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑c                                                                                                                                              <  5:201883/105‑1 (MQ=255)                                                      cACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCgg                                                                                                <  3:216273/149‑1 (MQ=255)                                                      cACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCgg                                                                                                <  8:16666/149‑1 (MQ=255)                                                      cACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCgg                                                                                                <  8:62276/149‑1 (MQ=255)                                                         cGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                                        >  6:30505/1‑138 (MQ=255)                                                         cGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                                        <  5:30505/138‑1 (MQ=255)                                                         cGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCT‑GCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTAc                                                                                            >  5:46513/1‑149 (MQ=255)                                                                gtggCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTgcgc                                                                                      >  3:79691/1‑149 (MQ=255)                                                                   gCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTgcgcgc                                                                                    <  7:239602/148‑1 (MQ=255)                                                                   gCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTgcgcgc                                                                                    <  5:134050/148‑1 (MQ=255)                                                                   gCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGcgcgcg                                                                                   <  4:74786/149‑1 (MQ=255)                                                                    ccTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCTCAGCCGCTGCAACCGGGTACCTGCGCGCgg                                                                                  <  1:206928/149‑1 (MQ=255)                                                                     ctcgctcgACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGa                                                                                 >  4:240417/1‑149 (MQ=255)                                                                     ctcgctcgACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGa                                                                                 >  3:254797/1‑149 (MQ=255)                                                                        gctcgACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCt                                                                             >  6:50916/1‑149 (MQ=255)                                                                            gACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTc                                                                           <  5:302907/147‑1 (MQ=255)                                                                                    gccCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGgcgcgc                                                                 >  8:60004/1‑149 (MQ=255)                                                                                       cTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCtg                                                              >  3:67993/1‑149 (MQ=255)                                                                                            cGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCg                                                         <  8:252231/149‑1 (MQ=255)                                                                                               aCAGCGCAATGGACGGTTATGCCTTGCGCCAGGCTGACTGGCAGGGCCAGCCGCTGACGGT‑CAGCCAGCGTATTTTCGCCGGCCATGAGCCCCAGCCGCTGAAACCGGGTACCTGCGCGCGCA‑TCTTCACCGGCGCGCCTGTGCCGGAg                                                      >  2:77615/1‑149 (MQ=255)                                                                                                 aGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAggg                                                    >  2:16100/1‑149 (MQ=255)                                                                                                  gcgcCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTACAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAggg                                                    >  1:94524/1‑149 (MQ=255)                                                                                                      cATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCACCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGcc                                                >  2:90024/1‑148 (MQ=255)                                                                                                           aCGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGc                                          <  4:67993/149‑1 (MQ=255)                                                                                                              gTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTc                                       <  4:254797/149‑1 (MQ=255)                                                                                                                tATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGAc                                             >  1:145565/1‑141 (MQ=255)                                                                                                                tATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGAc                                             <  2:145565/141‑1 (MQ=255)                                                                                                                    ccTTGCGCCT‑GCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATg                                 <  6:46513/148‑1 (MQ=255)                                                                                                                     cTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATg                                 >  4:181644/1‑148 (MQ=255)                                                                                                                     cTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATg                                 >  8:129331/1‑148 (MQ=255)                                                                                                                        gcgcCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAg                              <  7:60004/148‑1 (MQ=255)                                                                                                                          gcCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGa                            >  3:247860/1‑148 (MQ=255)                                                                                                                            cTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGaaaa                         >  4:24595/1‑149 (MQ=255)                                                                                                                            cTGGCGGACTGGCAGGGCCAGCCGCTGCCGGG‑CATCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTTCCTGCGCGCGGATTATTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGaaa                          <  1:90024/149‑1 (MQ=255)                                                                                                                                 tGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCGGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACAcc                     >  7:217600/1‑148 (MQ=255)                                                                                                                                 tGACCGGCAGCGGCTGGCCCTGCCAGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACAcc                     >  2:50868/1‑148 (MQ=255)                                                                                                                                   aCTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGaa                  <  3:240417/149‑1 (MQ=255)                                                                                                                                   aCTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGa                   <  5:305732/148‑1 (MQ=255)                                                                                                                                   aCTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGa                   <  5:50916/148‑1 (MQ=255)                                                                                                                                     tGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAg                 >  1:142746/1‑148 (MQ=255)                                                                                                                                       gCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTAc              >  6:283912/1‑149 (MQ=255)                                                                                                                                       gCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTAc              >  3:117776/1‑149 (MQ=255)                                                                                                                                         aGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTc            >  7:139969/1‑149 (MQ=255)                                                                                                                                          gggCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCg           <  8:139969/149‑1 (MQ=255)                                                                                                                                               aGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGTCTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGAt      <  1:16100/149‑1 (MQ=255)                                                                                                                                                 ccGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGATgg    >  6:33466/1‑149 (MQ=255)                                                                                                                                                   gCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGATGGcc  <  5:33466/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                       ATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGATGGCC  >  NC_002947/2421806‑2422100

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