Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,421,950	C→G	6.0%	P80A (CCG→GCG) ‡	moeA →	molybdenum::molybdopterin ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,421,950	0	C	G	6.0%	217.3 / 3.9	67	P80A (CCG→GCG) ‡	moeA	molybdenum::molybdopterin ligase
Reads supporting (aligned to +/- strand):  ref base C (36/27);  new base G (2/2);  total (38/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 = 3.36e-01

CCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGATGG  >  NC_002947/2421803‑2422098                                                                                                                                                    |                                                                                                                                                       ccGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGcc                                                                                                                                                       <  5:87608/149‑1 (MQ=255)    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                                                                                                                             <  6:302478/143‑1 (MQ=255)                                ccGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTc                                                                                                                             >  5:302478/1‑143 (MQ=255)                                ccGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGc                                                                                                                       >  8:239602/1‑149 (MQ=255)                                ccGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGc                                                                                                                       >  6:134050/1‑149 (MQ=255)                                 cGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGc                                                                                                                       >  4:100675/1‑148 (MQ=255)                                 cGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGc                                                                                                                       >  2:206928/1‑148 (MQ=255)                                  gAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGT‑ca                                                                                                                                           <  1:250570/127‑1 (MQ=38)                                  gAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGT‑ca                                                                                                                                           >  2:250570/1‑127 (MQ=38)                                        gAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGcc                                                                                                                                                       <  4:15214/110‑1 (MQ=255)                                        gAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGcc                                                                                                                                                       >  3:15214/1‑110 (MQ=255)                                          gggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                                  <  7:106845/113‑1 (MQ=255)                                          gggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                                  >  8:106845/1‑113 (MQ=255)                                               gcgTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCg                                                                                                        >  1:205246/1‑149 (MQ=255)                                                 gTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                                      >  7:252231/1‑149 (MQ=255)                                                 gTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                                      <  4:68539/149‑1 (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                                                                                              <  8:62276/149‑1 (MQ=255)                                                         cACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCgg                                                                                              <  3:216273/149‑1 (MQ=255)                                                         cACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCgg                                                                                              <  8:16666/149‑1 (MQ=255)                                                            cGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                                      <  5:30505/138‑1 (MQ=255)                                                            cGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                                      >  6:30505/1‑138 (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:50916/148‑1 (MQ=255)                                                                                                                                      aCTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGa                 <  5:305732/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)                                                                                                                                                    |                                                                                                                                                       CCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGATGG  >  NC_002947/2421803‑2422098

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