Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,669,976:1	+T	12.3%	coding (2342/2610 nt)	acnB ←	aconitate hydratase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,669,976	1	.	T	12.3%	189.2 / 26.4	57	coding (2342/2610 nt)	acnB	aconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base . (24/26);  new base T (4/3);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

CATGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCACC  >  NC_002947/2669828‑2670124                                                                                                                                                      |                                                                                                                                                     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)                                                                                                                                                     c‑gTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCAcc  <  8:125895/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CATGTACTCTTCGACGGTCGGCAGTTTGCCGATGATCGAAGCGACAGCGGCCAGCTCGGCCGATGCCAGGTACACGTTGGTCGCGTCGCCCAGACGGTTCGGGAAGTTACGGGTCGAGGTGGAGACCACGGTCGAACCGGTCTGCACGC‑GTGCCTGGTTGCCCATGCACAGCGAGCAGCCTGGCATTTCCATACGTGCACCGGCCTTGCCGTAGATGCCGTAGTAGCCTTCTTCGGTCAGCTGGTGAGCGTCCATCTTGGTTGGCGGAGCCAGCCACAGACGGGTCGGAATGCCACC  >  NC_002947/2669828‑2670124

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