Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	41,101	G→A	6.5%	P211P (CCG→CCA)	oprP →	porin P

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	41,101	0	G	A	6.5%	174.8 / 2.0	61	P211P (CCG→CCA)	oprP	porin P
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base A (2/2);  total (32/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 = 5.14e-01

CAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACG  >  NC_002947/40957‑41246                                                                                                                                                 |                                                                                                                                                  cAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGGGAAACCAATAACAACTCGGATGGCGACTGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGc                                                                                                                                               >  4:177664/1‑149 (MQ=255)   aCGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGca                                                                                                                                              >  1:11377/1‑148 (MQ=255)    cGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGca                                                                                                                                              >  1:138875/1‑147 (MQ=255)     gACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAg                                                                                                                                           >  1:154194/1‑149 (MQ=255)          cGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGGGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATg                                                                                                                                      <  1:224360/149‑1 (MQ=255)          cGCAGGCATGGGTGTTCAAGCCTCCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATCACAACTCGGATGGCGAGAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATg                                                                                                                                      <  3:13956/149‑1 (MQ=255)           gCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATgg                                                                                                                                     >  7:218583/1‑149 (MQ=255)                          cAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCACTGCACAGTGATGGCAACGTCCTGCACGt                                                                                                                      >  1:85967/1‑149 (MQ=255)                           aaGCCACCAGCACCATCGCTGACATGGCTTATGGTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGt                                                                                                                      <  2:287091/148‑1 (MQ=255)                               caccaGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCACTGCACAGTGATGGCAACGTCCTGCACGTCGGTg                                                                                                                 <  2:85967/149‑1 (MQ=255)                                accaGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCCCCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGc                                                                                                                >  7:142275/1‑149 (MQ=255)                                accaGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTcc                                                                                                                             <  6:228163/136‑1 (MQ=255)                                accaGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTcc                                                                                                                             >  5:228163/1‑136 (MQ=255)                                 ccaGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCt                                                                                                               >  4:142796/1‑149 (MQ=255)                                   aGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCa                                                                                                             <  8:218583/149‑1 (MQ=255)                                   aGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAAAAACTCGGAATGCGACAGCACCAAGCGCTACAACGTCCGTGGCGTGTTCTCAACGCTGCACACTGATGGCAACGTCCTGCACGTCGGTGCt                                                                                                               >  5:237221/1‑147 (MQ=255)                                   aGCACCATCGCTGACATGGCTTATGTTACCGGTACCGATTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACCTGGCCTGTTCGCACCGCTCCACAGAGATGGCAACGTCCTGCACGTCGGTGCt                                                                                                               >  5:237222/1‑147 (MQ=255)                                       ccATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTAc                                                                                                         >  5:212504/1‑149 (MQ=255)                                         aTCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGc                                                                                                       <  3:171805/149‑1 (MQ=255)                                         aTCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGc                                                                                                       <  2:154194/149‑1 (MQ=255)                                               gACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACcg                                                                                                 >  8:154511/1‑149 (MQ=255)                                                aCATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATg                                                                                                                                      >  3:280326/1‑111 (MQ=255)                                                aCATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATg                                                                                                                                      <  4:280326/111‑1 (MQ=255)                                                  aTGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGa                                                                                              <  1:203610/149‑1 (MQ=255)                                                    ggCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGAcc                                                                                            >  6:250487/1‑149 (MQ=255)                                                     gCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGAcc                                                                                            <  2:11377/148‑1 (MQ=255)                                                              tCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACa                                                                                   <  3:139115/148‑1 (MQ=255)                                                                  gTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCg                                                                              <  3:142796/149‑1 (MQ=255)                                                                  gTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCg                                                                              >  4:202579/1‑149 (MQ=255)                                                                     gCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGCGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTc                                                                           <  6:11224/149‑1 (MQ=255)                                                                       gTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGa                                                                         >  1:137938/1‑149 (MQ=255)                                                                          ttCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATa                                                                       <  7:154511/148‑1 (MQ=255)                                                                            cAGCTAAACCAATAACAACTCGGATGGCGTCAGCACCAAGCGCTTCAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAAGGTCGTGCACGTCGGTGCTCAGTACGCTTAGCGCGACCTCAAAAACAGTGCGGTGGATAcc                                                                     <  6:237222/148‑1 (MQ=255)                                                                            cAGCGAAAACAATAATACTTCGTATGGCGACAGCACAAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCGGCACAGTGATGTAAACGTCCTGCACGTCGGTGATCAGTACGATTACCGCGACCTCAAAAACAGTTCGGTAGATAcc                                                                     <  6:237221/148‑1 (MQ=255)                                                                              gCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGCCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCgg                                                                  <  2:138875/149‑1 (MQ=255)                                                                                      aaTAACAACTTGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCACTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGAttcgttc                                                          >  4:145318/1‑149 (MQ=255)                                                                                            aacTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCt                                                    >  2:315084/1‑149 (MQ=255)                                                                                             acTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTg                                                   <  5:51093/149‑1 (MQ=255)                                                                                                 ggATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCg                                               >  8:65939/1‑149 (MQ=255)                                                                                                       cGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCcgcg                                         <  3:202579/149‑1 (MQ=255)                                                                                                       cGACAGCACCAAGCGCTACAACGTACGT‑GCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCgg                                        >  6:31249/1‑149 (MQ=255)                                                                                                         aCAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGgt                                       >  6:231702/1‑149 (MQ=255)                                                                                                          cAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAATGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGgtg                                      <  5:250487/149‑1 (MQ=255)                                                                                                           aGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGgtgt                                     <  6:212504/149‑1 (MQ=255)                                                                                                               ccAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGAACTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGa                                  >  8:178589/1‑148 (MQ=255)                                                                                                                cAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATa                                >  3:60409/1‑149 (MQ=255)                                                                                                                cAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATa                                >  5:22145/1‑149 (MQ=255)                                                                                                                cAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCGCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCTTTCGCGTCTGGGCGTCCGCGGTGTAGATa                                <  8:142275/149‑1 (MQ=255)                                                                                                                cAAGCGCTACAACGTACGTGGCGTGTTAGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATa                                >  1:244642/1‑149 (MQ=255)                                                                                                                  aGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATAc                               >  6:44317/1‑148 (MQ=255)                                                                                                                    cgcTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCaa                            <  4:60409/149‑1 (MQ=255)                                                                                                                      cTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAg                                   >  5:180209/1‑140 (MQ=255)                                                                                                                      cTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAg                                   <  6:180209/140‑1 (MQ=255)                                                                                                                       tACAACGTACGT‑GCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATggc                        >  2:18901/1‑149 (MQ=255)                                                                                                                         cAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATggc                        <  2:137938/148‑1 (MQ=255)                                                                                                                                cgtGGCGTGTTCGCACCACTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTTCCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGa                <  3:145318/149‑1 (MQ=255)                                                                                                                                     cGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATggcggcgg                   >  7:219048/1‑141 (MQ=255)                                                                                                                                     cGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATggcggcgg                   <  8:219048/141‑1 (MQ=255)                                                                                                                                         ttCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCAc       <  5:44317/149‑1 (MQ=255)                                                                                                                                              aCCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACg  >  8:33258/1‑149 (MQ=255)                                                                                                                                              aCCGCTGCACAGTGATCGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACACGGATTCGTTCGCGTCTGGGAGTCCGCGGTGTAGATAACAATGGCGGCGGCGATGCCGGCACCAACg  >  4:266297/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  CAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACG  >  NC_002947/40957‑41246

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