Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 909,694 C→T 6.2% R68C (CGT→TGT)  fruB → PTS fructose transporter subunit EI/HPr/IIA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947909,6940CT 6.2% 192.8 / 2.4 64R68C (CGT→TGT) fruBPTS fructose transporter subunit EI/HPr/IIA
Reads supporting (aligned to +/- strand):  ref base C (29/31);  new base T (2/2);  total (31/33)
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.03e-01

GGCCGAGGCGTTACGCCTTCTGGCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCTCCGACGAACACCTGCGCCTGCTGCAGCTGT  >  NC_002947/909546‑909829
                                                                                                                                                    |                                                                                                                                       
ggCCGAGGCGTTACGCCTTCTGGCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAacac                                                                                                                                         <  7:198620/149‑1 (MQ=255)
   cGAGGCGTTACGCCTTCTGGCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTg                                                                                                                                      <  8:323633/149‑1 (MQ=255)
            aCGCCTTCTGGCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGt                                                                                                                                       <  8:233745/139‑1 (MQ=255)
            aCGCCTTCTGGCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGt                                                                                                                                       >  7:233745/1‑139 (MQ=255)
              gCCTTCTGGCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGt                                                                                                                                       <  5:231247/137‑1 (MQ=255)
              gCCTTCTGGCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGt                                                                                                                                       >  6:231247/1‑137 (MQ=255)
              gCCTTCTGGAGGATCGTCTTGTTGCGTACGGCCTGGTCGCCGAGGGTTCCCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTACACTTTCCTTGGCCAGGGCATCGCCATTCACCATGGTACGCCGCAAACACGTGACCTGGTGTAt                                                                                                                           <  4:14464/149‑1 (MQ=255)
               ctttCTGGCGTATGGGCTGGTTGCGGACGGCCTGGTCGCCTAGGTTTACGTGAAAGGGCTGAAGGCCGGGGAGGAAAAAGGCTCCATTTTCCTTGGCCAGGGCATCGCAATACCCCATGGTACGCCGCAAACACGTGACCTGGTGTATg                                                                                                                          <  1:234660/147‑1 (MQ=255)
                cTTCTGGCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGc                                                                                                                         >  7:136746/1‑149 (MQ=255)
                   cTGGCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGgtgt                                                                                                                             >  5:82871/1‑142 (MQ=255)
                   cTGGCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGgtgt                                                                                                                             <  6:82871/142‑1 (MQ=255)
                     ggCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCg                                                                                                                    >  4:285149/1‑149 (MQ=255)
                     ggCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCg                                                                                                                    >  6:170161/1‑149 (MQ=255)
                      gCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCAcc                                                                                                                     <  5:91183/147‑1 (MQ=255)
                         gATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTAt                                                                                                                           <  5:161002/138‑1 (MQ=255)
                         gATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTAt                                                                                                                           >  6:161002/1‑138 (MQ=255)
                             ggCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTAACTGCAAGGGCTGCAGGGCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGc                                                                                                            <  8:297014/149‑1 (MQ=255)
                                tGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCctg                                                                                                         <  8:136746/149‑1 (MQ=255)
                                       ggACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGt                                                                                                  >  6:33837/1‑149 (MQ=255)
                                            gcctGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGttccgc                                                                                             >  4:59851‑M2/5‑147 (MQ=255)
                                             cctGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTccc                                                                                              <  3:43193‑M2/144‑1 (MQ=255)
                                              ctGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACATGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTcccc                                                                                             >  8:248233‑M2/3‑147 (MQ=255)
                                              ctGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACATGGGACCTGGTGTATGCCACCGGCGTACGCCTGATGCAGTTcccc                                                                                             >  8:243516‑M2/3‑147 (MQ=255)
                                                    gCCGATGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGt                                                                                     <  7:313549/149‑1 (MQ=255)
                                                     ccGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTa                                                                                    <  5:170161/149‑1 (MQ=255)
                                                           gTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTgg                                                                              <  8:102249/149‑1 (MQ=255)
                                                               ccTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGt                                                                           >  1:165624/1‑148 (MQ=255)
                                                                cTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGa                                                                         >  6:127218/1‑149 (MQ=255)
                                                                  gCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACg                                                                       >  6:227680/1‑149 (MQ=255)
                                                                           gCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAg                                                                                          <  8:170524/121‑1 (MQ=255)
                                                                           gCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAg                                                                                          >  7:170524/1‑121 (MQ=255)
                                                                            cAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACATGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGt                                                             <  7:248233/149‑1 (MQ=255)
                                                                            cAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACATGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGt                                                             <  7:243516/149‑1 (MQ=255)
                                                                            cAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACgg                                                                      >  1:49077/1‑140 (MQ=255)
                                                                            cAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACgg                                                                      <  2:49077/140‑1 (MQ=255)
                                                                             aGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTc                                                            <  3:285149/149‑1 (MQ=255)
                                                                             aGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCAATCCCCATGGTACGCCGCAAACACGTGAGCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTc                                                            >  8:190198/1‑149 (MQ=255)
                                                                                  cGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGc                                                                    >  6:153551/1‑136 (MQ=255)
                                                                                  cGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGc                                                                    <  5:153551/136‑1 (MQ=255)
                                                                                        gCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCAt                                                 <  5:33837/149‑1 (MQ=255)
                                                                                         cacaAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCGGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCAtt                                                <  1:243279/149‑1 (MQ=255)
                                                                                          acaAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTg                                               <  2:165624/149‑1 (MQ=255)
                                                                                                cTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGTAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTg                                         >  2:27545/1‑149 (MQ=255)
                                                                                                 tCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTggg                                                                             >  2:74160/1‑112 (MQ=255)
                                                                                                 tCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTggg                                                                             <  1:74160/112‑1 (MQ=255)
                                                                                                 tCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTgc                                        >  2:119014/1‑149 (MQ=255)
                                                                                                  ccACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTgcc                                       >  5:234674/1‑149 (MQ=255)
                                                                                                    aCTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTgccgc                                     >  5:201887/1‑149 (MQ=255)
                                                                                                         ccTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCt                                >  7:196959/1‑149 (MQ=255)
                                                                                                          cTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGTAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCTc                               <  1:27545/149‑1 (MQ=255)
                                                                                                              gCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCTCcgac                           >  4:269167/1‑149 (MQ=255)
                                                                                                                       tCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCTCCGACGAACACCTg                  >  5:999/1‑149 (MQ=255)
                                                                                                                         gCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCTCCGACGAACACCTgc                 <  6:234674/148‑1 (MQ=255)
                                                                                                                          ccATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCTCCGACGAACACCTgcgc               >  7:209191/1‑149 (MQ=255)
                                                                                                                           cATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTggg                                                                             <  8:178478/86‑1 (MQ=255)
                                                                                                                           cATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTggg                                                                             >  7:178478/1‑86 (MQ=255)
                                                                                                                                cccATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTgg                                                     <  5:28165/105‑1 (MQ=255)
                                                                                                                                cccATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTgg                                                     >  6:28165/1‑105 (MQ=255)
                                                                                                                                 ccATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTAc                                                         >  1:119539/1‑100 (MQ=255)
                                                                                                                                 ccATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTAc                                                         <  2:119539/100‑1 (MQ=255)
                                                                                                                                  cATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCTCCGACGAACACCTGCGCCTGCTGCa       <  6:999/149‑1 (MQ=255)
                                                                                                                                      gTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCTCCGACGAACACCTGCGCCTGCTGCAGCt    <  5:227680/148‑1 (MQ=255)
                                                                                                                                       tACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCTCCGACGAACACCTGCGCCTGCTGCAGCTGt  <  6:201887/149‑1 (MQ=255)
                                                                                                                                               aaaCACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCTCcgac                           <  5:265890/116‑1 (MQ=255)
                                                                                                                                               aaaCACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCTCcgac                           >  6:265890/1‑116 (MQ=255)
                                                                                                                                                    |                                                                                                                                       
GGCCGAGGCGTTACGCCTTCTGGCGGATCGGCTGGTTGCGGACGGCCTGGTCGCCGAGGGTTACCTGCAAGGGCTGCAGGCCCGGGAGGCACAAGGCTCCACTTTCCTTGGCCAGGGCATCGCCATTCCCCATGGTACGCCGCAAACACGTGACCTGGTGTATGCCACCGGCGTACGCCTGCTGCAGTTCCCCGAGGGGGTAGACTGGGGTGACGGGCAAATGGTCTACCTGGCCATTGGTATTGCCGCCCGCTCCGACGAACACCTGCGCCTGCTGCAGCTGT  >  NC_002947/909546‑909829

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: ‑