Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,495,409 A→G 10.7% *206W (TGA→TGG PP_2188 → tRNA 2‑methylthio‑N6‑isopentenyl adenosine(37) hydroxylase MiaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,495,4090AG10.7% 191.6 / 12.4 75*206W (TGA→TGGPP_2188tRNA 2‑methylthio‑N6‑isopentenyl adenosine(37) hydroxylase MiaE
Reads supporting (aligned to +/- strand):  ref base A (34/32);  new base G (4/4);  total (39/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 = 8.21e-01

CTACCTGAAACTGGCGCACAACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCTTG  >  NC_002947/2495262‑2495552
                                                                                                                                                   |                                                                                                                                                
cTACCTGAAACTGGCGCACAACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGAt                                                                                                                                                 <  2:258197/149‑1 (MQ=255)
   ccTGAAACTGGCGCACAACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATggg                                                                                                                                              <  3:148100/149‑1 (MQ=255)
       aaaCTGGCGCACAACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCAc                                                                                                                                           >  2:351969/1‑148 (MQ=255)
                             atgaAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGcc                                                                                                                           >  1:238061/1‑141 (MQ=255)
                             atgaAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGcc                                                                                                                           <  2:238061/141‑1 (MQ=255)
                                  gCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCgg                                                                                                              <  1:152654/149‑1 (MQ=255)
                                    ggACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGcc                                                                                                            <  8:232679/149‑1 (MQ=255)
                                        aTCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGTTGGGCACCAGGCGGCGCAGGc                                                                                                                            >  5:313713/1‑129 (MQ=255)
                                        aTCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGc                                                                                                                            <  6:313713/129‑1 (MQ=255)
                                        aTCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCAACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAg                                                                                                         >  6:363084/1‑148 (MQ=255)
                                         tcgctcgcTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAg                                                                                                         >  1:297296/1‑147 (MQ=255)
                                         gcgcacgctGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAg                                                                                                         >  2:34791/2‑147 (MQ=255)
                                           gctcgctGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTaa                                                                                                                       >  6:163899/1‑131 (MQ=255)
                                           gctcgctGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTaa                                                                                                                       <  5:163899/131‑1 (MQ=255)
                                            ctcgctGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAg                                                                                                      >  3:336388/1‑147 (MQ=255)
                                               gcttcttCGAGATGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCGTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAg                                                                                                   <  3:312200/141‑1 (MQ=255)
                                                     tCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCg                                                                                           >  2:66453/1‑149 (MQ=255)
                                                       gAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGcc                                                                                         >  6:126525/1‑149 (MQ=255)
                                                            cgcaCGCCCGCCCGAAATGCAGCTGACCCAGGCCCCCGCCCATGCGCTGCCCGTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGccc                                                                                    <  8:304614/146‑1 (MQ=255)
                                                               aCGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAg                                                                                 >  1:27224/1‑149 (MQ=255)
                                                                cgcgcgGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAg                                                                                 <  6:173226/148‑1 (MQ=255)
                                                                  cgcgGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGcc                                                                              <  1:351969/149‑1 (MQ=255)
                                                                     ggcggAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGc                                                                           <  2:313336/149‑1 (MQ=255)
                                                                        ggAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGa                                                                        <  5:260366/149‑1 (MQ=255)
                                                                          aaaTGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCAACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGaaa                                                                      <  5:363084/149‑1 (MQ=255)
                                                                            aTGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGaaa                                                                      >  2:85690/1‑147 (MQ=255)
                                                                            aTGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCATAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATa                                                                    <  4:154012/149‑1 (MQ=255)
                                                                                aGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTAcc                                                                                            >  5:248570/1‑121 (MQ=255)
                                                                                aGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTAcc                                                                                            <  6:248570/121‑1 (MQ=255)
                                                                                 gCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGc                                                                                                     <  8:88005/111‑1 (MQ=255)
                                                                                   tGATCCAGTTACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGaaa                                                                      >  3:329506/1‑140 (MQ=255)
                                                                                   tGATCCAGTTACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGaaa                                                                      <  4:329506/140‑1 (MQ=255)
                                                                                      tCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCt                                                          <  2:297296/149‑1 (MQ=255)
                                                                                        cAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGc                                                        >  2:155044/1‑149 (MQ=255)
                                                                                         aGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAAATGCGCCGGCCCAGAAGCAGGCCTACCGCCAGACCGACGCCGGCAGAAAGAGCACCAGTCTGc                                                        >  8:144798/1‑148 (MQ=255)
                                                                                          gTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACGAGTCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCa                                                      <  7:232829/149‑1 (MQ=255)
                                                                                             accagaccagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGGCGGCAGAAATAGCAACAGTCTGCCAggg                                                   <  4:230989/149‑1 (MQ=255)
                                                                                             accagaccagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAggg                                                   >  3:180291/1‑149 (MQ=255)
                                                                                             accagaccagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAggg                                                   >  1:241847/1‑149 (MQ=255)
                                                                                             accagaccagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAggg                                                   >  3:230989/1‑149 (MQ=255)
                                                                                               cagaccagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCtt                                                                                                                         <  2:34924/77‑1 (MQ=255)
                                                                                               cagaccagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCtt                                                                                                                         >  1:34924/1‑77 (MQ=255)
                                                                                                agaccagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTg                                                >  6:183981/1‑149 (MQ=255)
                                                                                                  accagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGgc                                              >  5:172597/1‑149 (MQ=255)
                                                                                                         gCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGAt                                       <  1:34791/149‑1 (MQ=255)
                                                                                                                ttCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTcc                                >  2:237799/1‑149 (MQ=255)
                                                                                                                ttCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTcc                                >  4:304429/1‑149 (MQ=255)
                                                                                                                   cacaGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTg                              <  1:66453/148‑1 (MQ=255)
                                                                                                                    acaGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTggg                            >  8:10566/1‑149 (MQ=255)
                                                                                                                    acaGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTggg                            >  2:63616/1‑149 (MQ=255)
                                                                                                                        cGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGc                        <  4:180291/149‑1 (MQ=255)
                                                                                                                           cATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTc                     <  4:336388/149‑1 (MQ=255)
                                                                                                                           cATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTc                     <  1:85690/149‑1 (MQ=255)
                                                                                                                             tCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTcca                   >  5:222359/1‑149 (MQ=255)
                                                                                                                             tCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTcca                   >  8:156308/1‑149 (MQ=255)
                                                                                                                                cGCAGGCTTTGGCCGCCTGATGGGCACCAGGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGgccgcc                                          >  1:219688/1‑124 (MQ=255)
                                                                                                                                cGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGgccgcc                                          <  2:219688/123‑1 (MQ=255)
                                                                                                                                cGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTccacca                >  4:197881/1‑149 (MQ=255)
                                                                                                                                  cAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCaccacc              <  5:183981/149‑1 (MQ=255)
                                                                                                                                     gCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAg                                                     >  1:262645/1‑107 (MQ=255)
                                                                                                                                     gCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAg                                                     <  2:262645/107‑1 (MQ=255)
                                                                                                                                     gCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGt           <  7:144798/149‑1 (MQ=255)
                                                                                                                                     gCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGt           >  3:234811/1‑149 (MQ=255)
                                                                                                                                           gccgccTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGcc     >  6:103967/1‑149 (MQ=255)
                                                                                                                                            ccgccTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCt    <  3:304429/149‑1 (MQ=255)
                                                                                                                                            ccgccTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGACCAGTAGCAGGCCTACAGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGACAGGGGTGGCCGCCGATGTAGTCATGGGTGGCGTCCACCAACCGTTGCGACt    >  7:268778/1‑149 (MQ=255)
                                                                                                                                             cgccTGGTGCCCATCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATa                                                                    >  8:60579/1‑84 (MQ=39)
                                                                                                                                             cgccTGGTGCCCATCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATa                                                                    <  7:60579/84‑1 (MQ=39)
                                                                                                                                             cgccTGGTGCCCATCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCt                                                          <  5:96138/94‑1 (MQ=255)
                                                                                                                                             cgccTGGTGCCCATCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCt                                                          >  6:96138/1‑94 (MQ=255)
                                                                                                                                             cgccTGGTGCCCATCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTgcgc      >  8:13665/1‑146 (MQ=255)
                                                                                                                                             cgccTGGTGCCCATCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTgcgc      <  7:13665/146‑1 (MQ=255)
                                                                                                                                             cgccTGGTGCCCATCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCt  >  2:204773/1‑148 (MQ=255)
                                                                                                                                             cgccTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCCCCACCGGTTGCGCCtt   >  3:81461/1‑149 (MQ=255)
                                                                                                                                              gccTGGTGCCCATCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCTtg  <  1:204773/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTACCTGAAACTGGCGCACAACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCTTG  >  NC_002947/2495262‑2495552

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: