Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	727,935	C→A	5.8%	E145* (GAG→TAG)	bamD ←	outer membrane protein assembly factor BamD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	727,935	0	C	A	5.8%	217.9 / 3.0	69	E145* (GAG→TAG)	bamD	outer membrane protein assembly factor BamD
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base A (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

GGCAGCCGCTACATAAGCCTGGCGGCTCAGGTAGTAGTCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGCGTAATCGACGTTCGGGTGCTGCGGGTGCAGG  >  NC_002947/727786‑728080                                                                                                                                                      |                                                                                                                                                  ggCAGCCGCTACATAAGCCTGGCGGCTCAGGTAGTAGTCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTc                                                                                                                                                   <  2:1052889/150‑1 (MQ=255)   cAGCCGCTACATAAGCCTGGCGGCTCAGGTAGTAGTCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGt                                                                                                                                                 >  1:1195925/1‑150 (MQ=255)         cTACATAAGCCTGGCGGCTCAGGTAGTAGTCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACg                                                                                                                                           >  1:388800/1‑150 (MQ=255)          tACATAAGCCTGGCGGCTCAGGTAGTAGTCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTGGATCATGCGCTGCTTGGCGTCTGGCGCGTGGCGGCTGTTGGGGAAGCGGCTGGTAAGCTGGGCGAACTCGTTGTACGa                                                                                                                                          >  2:58262/1‑150 (MQ=255)          tACATAAGCCTGGCGGCTCAGGTAGTAGTCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGAGCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGa                                                                                                                                          >  2:58183/1‑150 (MQ=255)                aGCCTGGCGGCTCAGGTAGTAGTCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTcgcg                                                                                                                                    <  2:24893/150‑1 (MQ=255)                   cTGGCGGCTCAGGTAGTAGTCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCgcgggc                                                                                                                                 <  2:479034/150‑1 (MQ=255)                   cTGGCGGCTCAGGTAGTAGTCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCgcgggc                                                                                                                                 >  2:697984/1‑150 (MQ=255)                    tGGCGGCTCAGGTAATAGTCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCgcgggcg                                                                                                                                <  1:339343/150‑1 (MQ=255)                      gcggcTCAGGTAGTAGTCGGCCACACGATCCACGTACGAAGCCAGCAGGTTGAGCAGGTAGATAATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGATGGGGAAGCGGCTGGGCAGCTGGGCGAACTCGTTGTACGAGTCGCGggcggc                                                                                                                              >  1:31900/1‑150 (MQ=255)                                      tCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGCCGCCCgcgac                                                                                                                     <  2:552850/143‑5 (MQ=255)                                      tCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGCCGCCCgcgac                                                                                                                     >  1:552850/1‑139 (MQ=255)                                           cACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGa                                                                                                                                          >  1:90522/1‑117 (MQ=255)                                           cACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGa                                                                                                                                          <  2:90522/117‑1 (MQ=255)                                           cACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCAt                                                                                                         <  1:38872/150‑1 (MQ=255)                                           cACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCAt                                                                                                         >  1:911697/1‑150 (MQ=255)                                           cACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCAt                                                                                                         <  2:388800/150‑1 (MQ=255)                                                    ttCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCgg                                                                                                <  1:637387/150‑1 (MQ=255)                                                     tCGTACGAATCCAGCAGGTTGCGCAGGTAGAGCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTGGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGc                                                                                               <  1:1153004/150‑1 (MQ=255)                                                          ggAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGaa                                                                                          <  1:107880/149‑1 (MQ=255)                                                            aaGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCgcgggc                                                                                                                                 >  2:80724/1‑109 (MQ=255)                                                            aaGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCgcgggc                                                                                                                                 <  1:80724/109‑1 (MQ=255)                                                              gCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAgcgc                                                                                      <  2:1116553/150‑1 (MQ=255)                                                                cagcagGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAgcgcgc                                                                                    >  1:913732/1‑150 (MQ=255)                                                                 agcagGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTAGCGGGCGGCGCCCGGGTCACGCTTGGTCAt                                                                                                         >  1:60491/1‑128 (MQ=255)                                                                 agcagGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTAGCGGGCGGCGCCCGGGTCACGCTTGGTCAt                                                                                                         <  2:60491/128‑1 (MQ=255)                                                                   cagGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCcag                                                                                 <  2:1195925/150‑1 (MQ=255)                                                                           gcAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGTGGCAGTAGGGGTAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCAc                                                                         >  1:210479/1‑150 (MQ=255)                                                                             aGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTAc                                                                                                                                            <  2:294644/81‑1 (MQ=255)                                                                             aGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTAc                                                                                                                                            >  1:294644/1‑81 (MQ=255)                                                                             aGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACgg                                                                       <  2:913732/150‑1 (MQ=255)                                                                               gTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTc                                                                     >  2:157589/1‑150 (MQ=255)                                                                                          gctgctTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACg                                                          >  2:563783/1‑150 (MQ=255)                                                                                            tgctTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGt                                                                                                                      >  2:407059/1‑88 (MQ=255)                                                                                            tgctTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGt                                                                                                                      <  1:407059/88‑1 (MQ=255)                                                                                            tgctTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAg                                                        <  2:911697/150‑1 (MQ=255)                                                                                                   cGTCTGGCGCATAGCGGCTGTTGGGGAAGCGGCTGGTCAGTTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTg                                                                                                              >  1:1149747/1‑89 (MQ=255)                                                                                                   cGTCTGGCGCATAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTg                                                                                                              <  2:1149747/89‑1 (MQ=255)                                                                                                      gtGGCCCGTAGCGGCTGTTGGGTAAGCGGCTGGTCAGCGGGGTGAACTCGGTGTACGAGTCGCGGGCGGCGCTCGGGTCACGCTTGGTCGTGTCCAGCGGCAGGAAGCGCGCCAGCACGCCACGGTCCTGGTCGAACGAGGTCAGGCCCt                                              <  2:61699/149‑1 (MQ=255)                                                                                                          cgcgTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCCTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAAGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAg                                          <  1:697984/150‑1 (MQ=255)                                                                                                           gcgTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAg                                          <  1:231635/149‑1 (MQ=255)                                                                                                           gcgTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAg                                          >  2:231635/1‑149 (MQ=255)                                                                                                                gCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGgtagta                                    <  2:762352/150‑1 (MQ=255)                                                                                                                gCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGgtagta                                    >  1:90168/1‑150 (MQ=255)                                                                                                                gCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGgtagta                                    >  1:806031/1‑150 (MQ=255)                                                                                                                 cGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGTCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGACCACCGCTAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGtagtag                                   <  1:633590/150‑1 (MQ=255)                                                                                                                 cGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGAAGGCAGCGCGCAACCAGCACAAGATGCTCGGCTAACAAGGTCATGCCCTTGAGGtagtag                                   >  2:132403/1‑150 (MQ=255)                                                                                                                   ggtggtGGGGAAGGGCCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCAAGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGc                                 <  1:58262/145‑1 (MQ=255)                                                                                                                   gCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTcc                                                                    >  1:852071/1‑115 (MQ=255)                                                                                                                   gCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTcc                                                                    <  2:852071/115‑1 (MQ=255)                                                                                                                   gCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGc                                 <  1:58183/150‑1 (MQ=255)                                                                                                                     tGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCAt                                                                                                         >  2:208681/1‑76 (MQ=255)                                                                                                                     tGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCAt                                                                                                         <  1:208681/76‑1 (MQ=255)                                                                                                                     tGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGCgg                               >  1:505598/1‑149 (MQ=255)                                                                                                                          gggAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGTGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAGGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGAATGAGGCGTAATCg                          >  1:548273/1‑150 (MQ=255)                                                                                                                          gggAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGCGTAATCg                          >  2:341930/1‑150 (MQ=255)                                                                                                                            gAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGTTCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGCGTAATCGAc                        >  1:1181721/1‑150 (MQ=255)                                                                                                                                   ttGGTCAGCGGGCCGAACTCGTTGTACGAGTCGCGTGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGCGTAATCGACGTTCggg                 <  2:548273/149‑1 (MQ=255)                                                                                                                                   cTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGCGTAATCGACGTTCggg                 >  1:856931/1‑150 (MQ=255)                                                                                                                                   cTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGCGTAATCGACGTTCcgg                 <  2:806031/150‑1 (MQ=255)                                                                                                                                    tGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGCGTAATCGACGTTCGGGt                <  2:210479/150‑1 (MQ=255)                                                                                                                                         aGCTGGGCGAACTAGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAg                                          <  1:88238/119‑1 (MQ=255)                                                                                                                                         aGCTGGGCGAACTAGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAg                                          >  2:88238/1‑119 (MQ=255)                                                                                                                                         aGCTGGGCGAACTAGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAg                                          >  2:87904/1‑119 (MQ=255)                                                                                                                                         aGCTGGGCGAACTAGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAg                                          <  1:87904/119‑1 (MQ=255)                                                                                                                                           cTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGGAGTAGGTGTAATCGACGCTCGGGTGCTGCgg         >  2:152400/1‑150 (MQ=255)                                                                                                                                           cTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTACCACATGGTAATGTCCAGGGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGCGTAATCGACGTTCGGGTGCTGCgg         <  2:31900/150‑1 (MQ=255)                                                                                                                                            tGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGCGTAATCGACGTTCGGGTGCTGCggg        >  2:102634/1‑150 (MQ=255)                                                                                                                                             gggCGAACTCGTTGTCCGAGTCGCGGGCGGCGCCCGGGTGACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGCGTAATCGACGTTCGGGTGCTGCGGGt       <  2:89926/150‑1 (MQ=255)                                                                                                                                                  taCTCGTTGTCCGAGTCGCGGGCGGTGCACGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGTCCCTTGGGGTAGTAGGCGTAATCGACGTTCGGGTGCTGCGGGTGCAgg  <  2:505598/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  GGCAGCCGCTACATAAGCCTGGCGGCTCAGGTAGTAGTCGGCCACGTGGATTTCGTACGAAGCCAGCAGGTTGCGCAGGTAGATCATGCGCTGCTTGGCGTCTGGCGCGTAGCGGCTGTTGGGGAAGCGGCTGGTCAGCTGGGCGAACTCGTTGTACGAGTCGCGGGCGGCGCCCGGGTCACGCTTGGTCATGTCCAGCGGCAGGAAGCGCGCCAGCAGGCCACGGTCCTGGTCGAACGAGGTCAGGCCCTTGAGGTAGTAGGCGTAATCGACGTTCGGGTGCTGCGGGTGCAGG  >  NC_002947/727786‑728080

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