Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,669,973	A→T	12.7%	V782E (GTG→GAG)	acnB ←	aconitate hydratase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,669,973	0	A	T	12.7%	158.9 / 9.8	61	V782E (GTG→GAG)	acnB	aconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base T (4/4);  total (31/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 = 1.00e+00

CATGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGC  >  NC_002947/2669828‑2670120                                                                                                                                                  |                                                                                                                                                     cATGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCAcgc                                                                                                                                                   >  3:110796/1‑149 (MQ=255) cATGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCAc                                                                                                                                                     <  4:187200/147‑1 (MQ=255)  aTGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCAcgc‑g                                                                                                                                                 >  8:341175/1‑149 (MQ=255)  aTGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCAcgc‑g                                                                                                                                                 >  7:60081/1‑149 (MQ=255)      aCTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGTCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGcc                                                                                                                                             >  3:298896/1‑149 (MQ=255)          ttCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGcc                                                                                                                                             <  8:287529/145‑1 (MQ=255)          ttCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGcc                                                                                                                                             >  7:287529/1‑145 (MQ=255)             gACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCTCGC‑GTGCCTGGTTGc                                                                                                                                      <  2:336545/149‑1 (MQ=255)               cGGTCGGCAGTTTGCCGTTGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGTTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGccc                                                                                                                                    <  1:312974/149‑1 (MQ=255)               cGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGccc                                                                                                                                    <  4:315404/149‑1 (MQ=255)                 gtcaagggtggcattccgacccgtctgtggctggctccgccaaccaagatggacgctcaccagctgaccgaagaaggctactacggcatctacggcaaggccggtgcacgtatggaaatgccaggCTGCTCACTGTGCCTGGTTGCCCa                                                                                                                                   <  4:64336‑M1/24‑1 (MQ=255)                 gtcaagggtggcattccgacccgtctgtggctggctccgccaaccaagatggacgctcaccagctgaccgaagaaggctactacggcatctacggcaaggccggtgcacgtatggaaatgccaggCTGCTCACTGTGCCTGGTTGCCCa                                                                                                                                   <  4:69115‑M1/24‑1 (MQ=255)                   cGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATg                                                                                                                                 >  1:335614/1‑148 (MQ=255)                    ggCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGca                                                                                                                               >  6:75788/1‑149 (MQ=255)                      cAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGcaca                                                                                                                             <  7:341175/149‑1 (MQ=255)                       ggtggcattccgacccgtctgtggctggctccgccaaccaagatggacgctcaccagctgaccgaagaaggctactacggcatctacggcaaggccggtgcacgtatggaaatgccaggCTGCTCGCTGTGCCTGGTTGCCCATGcaca                                                                                                                             >  4:271276‑M1/120‑149 (MQ=255)                       ggtggcattccgacccgtctgtggctggctccgccaaccaagatggacgctcaccagctgaccgaagaaggctactacggcatctacggcaaggccggtgcacgtatggaaatgccaggCTGCTCACTGTGCCTGGTTGCCCATGcaca                                                                                                                             >  3:64336‑M1/120‑149 (MQ=255)                       ggtggcattccgacccgtctgtggctggctccgccaaccaagatggacgctcaccagctgaccgaagaaggctactacggcatctacggcaaggccggtgcacgtatggaaatgccaggCTGCTCACTGTGCCTGGTTGCCCATGcaca                                                                                                                             >  3:69115‑M1/120‑149 (MQ=255)                        gTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGc                                                                                                                           <  6:193165/149‑1 (MQ=255)                           tGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGTCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGa                                                                                                                         <  4:298896/148‑1 (MQ=255)                           tGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGa                                                                                                                         <  3:327255/148‑1 (MQ=255)                             attccgacccgtctgtggctggctccgccaaccaagatggacgctcaccagctgaccgaagaaggctactacggcatctacggcaaggccggtgcacgtatggaaatgccaggCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGagca                                                                                                                      >  7:223069‑M1/114‑149 (MQ=255)                                   atCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTgg                                                                                                                >  2:176078/1‑149 (MQ=255)                                   atCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTg                                                                                                                 <  2:335614/148‑1 (MQ=255)                                     cGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCa                                                                                                              >  3:36527/1‑149 (MQ=255)                                       gtctgtggctggctccgccaaccaagatggacgctcaccagctgaccgaagaaggctactacggcatctacggcaaggccggtgcacgtatggaaatgccaggCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCt                                                                                                                  <  6:320255‑M1/40‑1 (MQ=255)                                       gtctgtggctggctccgccaaccaagatggacgctcaccagctgaccgaagaaggctactacggcatctacggcaaggccggtgcacgtatggaaatgccaggCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCt                                                                                                                  >  5:320255‑M1/104‑143 (MQ=255)                                           gACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGc                                                                                                                           <  8:107702/130‑1 (MQ=255)                                           gACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGc                                                                                                                           >  7:107702/1‑130 (MQ=255)                                              aGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATAc                                                                                                     <  8:60081/149‑1 (MQ=255)                                                cGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGt                                                                                                   >  2:236090/1‑149 (MQ=255)                                                cGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGt                                                                                                   <  6:212730/149‑1 (MQ=255)                                                  gCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTg                                                                                                  >  6:155515/1‑148 (MQ=255)                                                      gCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCg                                                                                             >  6:254956/1‑149 (MQ=255)                                                           gCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCtt                                                                                        <  1:176078/149‑1 (MQ=255)                                                            ccGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTg                                                                                       >  7:148940/1‑149 (MQ=255)                                                              gATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGcc                                                                                     <  8:149314/149‑1 (MQ=255)                                                                  ccAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTa                                                                                  <  4:36527/148‑1 (MQ=255)                                                                   cAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGGATTTCCATACGTGCACCGGCCTTGCCg                                                                                    >  4:209984/1‑145 (MQ=255)                                                                   cAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCg                                                                                    >  3:311653/1‑145 (MQ=255)                                                                     ggTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATACACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATg                                                                              <  4:110796/149‑1 (MQ=255)                                                                      gTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGc                                                                             >  3:111375/1‑149 (MQ=255)                                                                          acGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATg                                                                                                                                 >  8:280815/1‑93 (MQ=255)                                                                           cGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCgtag                                                                        <  4:311653/149‑1 (MQ=255)                                                                                gtcgcgtcgcCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGcc                                                                   >  7:320026/1‑149 (MQ=255)                                                                                     gtcgcCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGAACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCTTAGTAGCcttctt                                                              >  4:53075/1‑149 (MQ=255)                                                                                      tcgcCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCttcttc                                                             <  5:75788/149‑1 (MQ=255)                                                                                                  ttCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCcttctt                                                              >  4:281169/1‑136 (MQ=255)                                                                                                  ttCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCcttctt                                                              <  3:281169/136‑1 (MQ=255)                                                                                                         aGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGGTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTcc                                          <  5:254956/149‑1 (MQ=255)                                                                                                             aCGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCt                                      >  6:243008/1‑149 (MQ=255)                                                                                                               gggTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCttg                                    <  1:236090/149‑1 (MQ=255)                                                                                                                ggTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCttgg                                   <  7:327053/149‑1 (MQ=255)                                                                                                                ggTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCttgg                                   <  8:320026/149‑1 (MQ=255)                                                                                                                             cgtatggaaatgccaggCTGCTCGCTGTGCATGGGCAACCAGGCAC‑GCGTGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTGTTCGGTCAGGTGGTGAGCGTCCATCTTGGTTGGCGGagccag                      <  3:15205‑M1/132‑1 (MQ=255)                                                                                                                             cgtatggaaatgccaggCTGCTCGCTGTGCATGGGCAACCAGGCAC‑GCGTGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGagc                         >  4:15205‑M1/18‑146 (MQ=255)                                                                                                                             aCCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGagccag                      <  5:243008/149‑1 (MQ=255)                                                                                                                               cACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGagccagc                     >  8:314398/1‑148 (MQ=255)                                                                                                                                        aCCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACggg           >  7:288934/1‑149 (MQ=255)                                                                                                                                         ccGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGt          >  7:125895/1‑149 (MQ=255)                                                                                                                                                 CACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGc  <  2:212505‑M1/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CATGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGC  >  NC_002947/2669828‑2670120

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