Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,495,413 G→C 10.8% intergenic (+4/+19) PP_2188 → / ← PP_2189 tRNA 2‑methylthio‑N6‑isopentenyl adenosine(37) hydroxylase MiaE/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,495,4130GC10.8% 216.9 / 14.6 74intergenic (+4/+19)PP_2188/PP_2189tRNA 2‑methylthio‑N6‑isopentenyl adenosine(37) hydroxylase MiaE/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (36/30);  new base C (4/4);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

CCTGAAACTGGCGCACAACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCTTGCTGCGCACG  >  NC_002947/2495265‑2495561
                                                                                                                                                    |                                                                                                                                                     
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                                                                               <  4:329506/140‑1 (MQ=255)
                                                                                tGATCCAGTTACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGaaa                                                                               >  3:329506/1‑140 (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:230989/1‑149 (MQ=255)
                                                                                          accagaccagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAggg                                                            >  1:241847/1‑149 (MQ=255)
                                                                                          accagaccagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAggg                                                            >  3:180291/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                                         >  4:304429/1‑149 (MQ=255)
                                                                                                             ttCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTcc                                         >  2:237799/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                            >  8:156308/1‑149 (MQ=255)
                                                                                                                          tCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTcca                            >  5:222359/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                    >  3:234811/1‑149 (MQ=255)
                                                                                                                                  gCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGt                    <  7:144798/149‑1 (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)
                                                                                                                                          cgccTGGTGCCCATCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCtg            >  2:204773/1‑148 (MQ=255)
                                                                                                                                          cgccTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCCCCACCGGTTGCGCCtt            >  3:81461/1‑149 (MQ=255)
                                                                                                                                           gccTGGTGCCCATCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCTtg           <  1:204773/149‑1 (MQ=255)
                                                                                                                                                 tGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCTTGCTgcg      >  5:166574/1‑148 (MQ=255)
                                                                                                                                                    gCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCTTGCTGCGCACg  >  5:39427/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTGAAACTGGCGCACAACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCTTGCTGCGCACG  >  NC_002947/2495265‑2495561

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: