Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,669,975	G→A	7.0%	R781R (CGC→CGT)	acnB ←	aconitate hydratase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,669,975	0	G	A	7.0%	167.7 / 2.9	57	R781R (CGC→CGT)	acnB	aconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base A (2/2);  total (28/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 = 9.97e-01

CATGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCAC  >  NC_002947/2669828‑2670123                                                                                                                                                    |                                                                                                                                                      cATGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCAcgc                                                                                                                                                      >  3:110796/1‑149 (MQ=255)  aTGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCAcgc‑g                                                                                                                                                    >  7:60081/1‑149 (MQ=255)  aTGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCAcgc‑g                                                                                                                                                    >  8:341175/1‑149 (MQ=255)      aCTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGTCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGcc                                                                                                                                                >  3:298896/1‑149 (MQ=255)          ttCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGcc                                                                                                                                                >  7:287529/1‑145 (MQ=255)          ttCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGcc                                                                                                                                                <  8:287529/145‑1 (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:69115‑M1/24‑1 (MQ=255)                 gtcaagggtggcattccgacccgtctgtggctggctccgccaaccaagatggacgctcaccagctgaccgaagaaggctactacggcatctacggcaaggccggtgcacgtatggaaatgccaggCTGCTCACTGTGCCTGGTTGCCCa                                                                                                                                      <  4:64336‑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                                                                                                                     >  5:320255‑M1/104‑143 (MQ=255)                                       gtctgtggctggctccgccaaccaagatggacgctcaccagctgaccgaagaaggctactacggcatctacggcaaggccggtgcacgtatggaaatgccaggCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCt                                                                                                                     <  6:320255‑M1/40‑1 (MQ=255)                                           gACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGc                                                                                                                              >  7:107702/1‑130 (MQ=255)                                           gACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGc                                                                                                                              <  8:107702/130‑1 (MQ=255)                                              aGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATAc                                                                                                        <  8:60081/149‑1 (MQ=255)                                                cGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGt                                                                                                      <  6:212730/149‑1 (MQ=255)                                                cGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGt                                                                                                      >  2:236090/1‑149 (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)                                                                                                                                                   cgc‑gTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGGATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCa   <  3:209984/149‑1 (MQ=255)                                                                                                                                                   cgc‑gTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAAt       >  7:76516/1‑145 (MQ=255)                                                                                                                                                   cgc‑gTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAAt       <  8:76516/145‑1 (MQ=255)                                                                                                                                                    gc‑gTGCCTGGTTGCCCCTGCACAGCGAGCAGCCTGGCATTTCCATACGGGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCAc  <  4:194626/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CATGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCAC  >  NC_002947/2669828‑2670123

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