Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,719,298	Δ1 bp	10.2%	coding (382/585 nt)	nfuA →	Fe/S biogenesis protein NfuA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,719,298	0	A	.	10.2%	173.8 / 14.7	59	coding (382/585 nt)	nfuA	Fe/S biogenesis protein NfuA
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base . (4/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01

TGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCT  >  NC_002947/2719155‑2719446                                                                                                                                                |                                                                                                                                                      tGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATg                                                                                                                                                 <  1:600726/150‑1 (MQ=255)         tACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCCCCCCAGGCCAAGGTGCCGATGGTCAACGAAGAAGGCGCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCC‑GCCATGGCTGGCCAg                                                                                                                                        <  1:570987/150‑1 (MQ=255)            gCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTg                                                                                                                                     >  2:443297/1‑150 (MQ=255)              cACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCC‑GCCATGGCTGGCCa                                                                                                                                         <  2:387635/144‑1 (MQ=255)              cACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCC‑GCCATGGCTGGCCa                                                                                                                                         >  1:387635/1‑144 (MQ=255)                      gCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCtggtgg                                                                                                                           >  1:395217/1‑150 (MQ=255)                               gCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGgcca                                                                                                                                                       <  1:1088949/114‑1 (MQ=255)                               gCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGgcca                                                                                                                                                       >  2:1088949/1‑114 (MQ=255)                                 cAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGAt                                                                                                                >  1:788053/1‑150 (MQ=255)                                   gCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGa                                                                                                              <  2:1198691/150‑1 (MQ=255)                                     tGACCATCAAGGCGCCAAACGCAAAGGTGACCATGGTAAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACg                                                                                                            <  2:800534/150‑1 (MQ=255)                                      gACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACgg                                                                                                           >  1:982088/1‑150 (MQ=255)                                         cATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCtg                                                                                                                               >  2:2689/1‑127 (MQ=255)                                         cATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCtg                                                                                                                               <  1:2689/127‑1 (MQ=255)                                          aTCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTg                                                                                                                      <  2:670485/135‑1 (MQ=255)                                          aTCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTg                                                                                                                      >  1:670485/1‑135 (MQ=255)                                          aTCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCAtt                                                                                                       >  2:1017675/1‑150 (MQ=255)                                                     aaaCGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCtggtg                                                                                                                            <  1:90633/118‑1 (MQ=255)                                                     aaaCGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCtggtg                                                                                                                            >  2:90633/1‑118 (MQ=255)                                                      aaCGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAg                                                                                           >  2:1119461/1‑150 (MQ=255)                                                           cAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGTGGGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCgg                                                                                      >  1:208981/1‑150 (MQ=255)                                                              ggTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGgcca                                                                                                                                                       <  1:1192632/83‑1 (MQ=255)                                                              ggTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGgcca                                                                                                                                                       >  2:1192632/1‑83 (MQ=255)                                                              ggTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGtt                                                                                         >  1:758923/1‑144 (MQ=255)                                                              ggTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGtt                                                                                         <  2:758923/144‑1 (MQ=255)                                                               gTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggtggt                                                                                  <  2:628530/150‑1 (MQ=255)                                                                    gATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGcc                                                                                                    >  2:770697/1‑127 (MQ=255)                                                                    gATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGcc                                                                                                    <  1:770697/127‑1 (MQ=255)                                                                    gATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTg                                                                             >  2:1101073/1‑150 (MQ=255)                                                                         tCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGgtggtg                                                                                 >  2:601021/1‑141 (MQ=255)                                                                         tCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGgtggtg                                                                                 <  1:601021/141‑1 (MQ=255)                                                                           aaCGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGt                                                                      >  2:371408/1‑150 (MQ=255)                                                                            aCGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGtt                                                                     >  2:80099/1‑150 (MQ=255)                                                                             cGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTAGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTg                                                                    >  2:105851/1‑150 (MQ=255)                                                                                aGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATggc‑ggc                                                                                                                                           <  1:692700/76‑1 (MQ=255)                                                                                aGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATggc‑ggc                                                                                                                                           >  2:692700/1‑76 (MQ=255)                                                                                      cccGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGc                                                           <  2:395217/150‑1 (MQ=255)                                                                                      cccGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGc                                                           >  2:908835/1‑150 (MQ=255)                                                                                       ccGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGcc                                                          >  1:964212/1‑150 (MQ=255)                                                                                       ccGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGcc                                                          >  2:518578/1‑150 (MQ=255)                                                                                          aTCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGAt                                                       <  1:1119461/150‑1 (MQ=255)                                                                                          aTCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGAt                                                       <  1:443297/150‑1 (MQ=255)                                                                                            cAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGt                                                     <  1:371408/150‑1 (MQ=255)                                                                                                gaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAg                                                                         >  1:730687/1‑126 (MQ=255)                                                                                                gaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAg                                                                         <  2:730687/126‑1 (MQ=255)                                                                                                     tATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaa                                            <  1:518578/150‑1 (MQ=255)                                                                                                     tATCAAAAACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaa                                            >  2:951732/1‑150 (MQ=255)                                                                                                         aaCTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGCTGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaaggaa                                        >  2:750131/1‑150 (MQ=255)                                                                                                          actactACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaaggaag                                       <  2:964212/150‑1 (MQ=255)                                                                                                          actactACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaaggaag                                       <  1:239596/150‑1 (MQ=255)                                                                                                                  ctacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGACGGAAGGCATCGac                               >  2:1063389/1‑149 (MQ=255)                                                                                                                   tacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAgc                              <  1:951732/150‑1 (MQ=255)                                                                                                                        aCCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCAcc                         >  2:232819/1‑150 (MQ=255)                                                                                                                         ccGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCAccc                        <  1:908835/150‑1 (MQ=255)                                                                                                                              aTCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTc                   <  1:679607/150‑1 (MQ=255)                                                                                                                                  aCCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCAc                                                  <  1:763244/115‑1 (MQ=255)                                                                                                                                  aCCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCAc                                                  >  2:763244/1‑115 (MQ=255)                                                                                                                                   cgccGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCg              >  2:464064/3‑150 (MQ=255)                                                                                                                                      tcaccTGGCC‑GCCATGGCTGGCCAGGCGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTAt           >  1:799744/4‑150 (MQ=255)                                                                                                                                        gCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGGGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTc         <  1:80099/150‑1 (MQ=255)                                                                                                                                         ccTGGCC‑GCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTcc        >  2:215765/1‑150 (MQ=255)                                                                                                                                          cTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCg       <  2:982088/150‑1 (MQ=255)                                                                                                                                          cTGGCC‑GCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCg       >  1:173000/1‑150 (MQ=255)                                                                                                                                            gtccagccaTGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGa     <  1:464064/148‑1 (MQ=255)                                                                                                                                               cagccaTGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGa     <  2:692180/147‑1 (MQ=255)                                                                                                                                               cagccaTGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGa     >  1:692180/1‑147 (MQ=255)                                                                                                                                               cagccaTGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCt  >  1:344526/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                      TGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGC‑GGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCT  >  NC_002947/2719155‑2719446

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑