Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,421,948	A→T	5.8%	Q79L (CAG→CTG)	moeA →	molybdenum::molybdopterin ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,421,948	0	A	T	5.8%	215.2 / 2.4	69	Q79L (CAG→CTG)	moeA	molybdenum::molybdopterin ligase
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base T (2/2);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

CCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGAT  >  NC_002947/2421803‑2422096                                                                                                                                                  |                                                                                                                                                       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                                         <  2:145565/141‑1 (MQ=255)                                                                                                                   tATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGAc                                         >  1:145565/1‑141 (MQ=255)                                                                                                                       ccTTGCGCCT‑GCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATg                             <  6:46513/148‑1 (MQ=255)                                                                                                                        cTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATg                             >  8:129331/1‑148 (MQ=255)                                                                                                                        cTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATg                             >  4:181644/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)                                                                                                                                                  |                                                                                                                                                       CCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGT‑CAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGA‑TCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCGAAGAT  >  NC_002947/2421803‑2422096

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