Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,049,292 G→A 6.1% V52M (GTG→ATG)  spoT → bifunctional (p)ppGpp synthase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,049,2920GA 6.1% 193.9 / 2.0 65V52M (GTG→ATG) spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base A (2/2);  total (33/32)
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.06e-01

TATAGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCA  >  NC_002947/6049147‑6049439
                                                                                                                                                 |                                                                                                                                                   
tataGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTgg                                                                                                                                                  <  7:304089/149‑1 (MQ=255)
  taGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGc                                                                                                                                                 <  5:299011/148‑1 (MQ=255)
       gCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCAt                                                                                                                                           >  5:140183/1‑149 (MQ=255)
             gcggAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTgg                                                                                                                                                  <  6:293426/133‑1 (MQ=255)
             gCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTgg                                                                                                                                                  >  5:293426/1‑136 (MQ=255)
              ccGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTg                                                                                                                                                   <  4:50924/134‑1 (MQ=255)
              ccGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTg                                                                                                                                                   >  3:50924/1‑134 (MQ=255)
               cGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgg                                                                                                                                           >  5:155008/1‑139 (MQ=255)
               cGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgg                                                                                                                                           <  6:155008/141‑3 (MQ=255)
                gAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGa                                                                                                                                  >  7:225959/1‑149 (MQ=255)
                            aCCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGc                                                                                                                             >  6:66212/1‑142 (MQ=255)
                            aCCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGGGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGc                                                                                                                             <  5:66212/142‑1 (MQ=255)
                              cTATCTTGGCCCCGAACAGGTCAACGTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgggtgg                                                                                                                                       <  6:287040/130‑7 (MQ=255)
                              cTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGAc                                                                                                                     >  1:344234/1‑148 (MQ=255)
                              cTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgggtgg                                                                                                                                       >  5:287040/1‑124 (MQ=255)
                                   ttGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCa                                                                                                                                            <  4:55421/120‑1 (MQ=255)
                                     ggCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCagag                                                                                                             >  6:193188/1‑149 (MQ=255)
                                       ccccGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGcc                                                                                                                                    <  7:146997/124‑1 (MQ=255)
                                       ccccGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGcc                                                                                                                                    >  8:146997/1‑124 (MQ=255)
                                       ccccGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGcc                                                                                                                                    <  2:128200/124‑1 (MQ=255)
                                       ccccGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGcc                                                                                                                                    >  1:128200/1‑124 (MQ=255)
                                           gAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGAt                                                                                                       >  6:163379/1‑149 (MQ=255)
                                                     aCCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCAtg                                                                                             <  4:234932/149‑1 (MQ=255)
                                                      ccTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATc                                                                                                                 <  7:170814/128‑1 (MQ=255)
                                                      ccTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATc                                                                                                                 >  8:170814/1‑128 (MQ=255)
                                                      ccTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCAtgc                                                                                            >  1:15633/1‑149 (MQ=255)
                                                          gTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCAtgctgc                                                                                         <  8:225959/148‑1 (MQ=255)
                                                               gCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtg                                                                                   <  8:84969/149‑1 (MQ=255)
                                                                cGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtgt                                                                                  <  2:344234/149‑1 (MQ=255)
                                                                     cTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc                                                                              >  6:849/1‑148 (MQ=255)
                                                                     cTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCg                                                                             <  5:163379/149‑1 (MQ=255)
                                                                      tATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCg                                                                             >  5:53017/1‑148 (MQ=255)
                                                                        tttCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGaa                                                                           <  5:193188/148‑1 (MQ=255)
                                                                          tCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGa                                                                         <  5:849/148‑1 (MQ=255)
                                                                              cGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCg                                                                    >  5:176542/1‑149 (MQ=255)
                                                                               gCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCgg                                                                   >  1:219749/1‑149 (MQ=255)
                                                                                      aGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGcc                                                            <  2:15633/149‑1 (MQ=255)
                                                                                        gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATggcg                                                                                                   >  5:173930/1‑108 (MQ=255)
                                                                                        gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATggcg                                                                                                   <  6:173930/108‑1 (MQ=255)
                                                                                        gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCaa                                                          <  2:219749/149‑1 (MQ=255)
                                                                                              gATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGa                                                                                <  6:4620/121‑1 (MQ=255)
                                                                                                 gggCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCGCATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCt                                                 >  4:233341/1‑149 (MQ=255)
                                                                                                  ggCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCt                                                 <  6:140183/148‑1 (MQ=255)
                                                                                                   gCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTca                                               <  6:53017/149‑1 (MQ=255)
                                                                                                   gCAACGCCGCCGCACTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCACCCCGACATGCACATGGACCATCAGACCCTGATGGCGGCCATGCTGCACGATCTGATCCAAGACACCGGCATCGCAAAGGAACCCCt                                                 >  8:234787/1‑147 (MQ=255)
                                                                                                    cAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGCCACCGGCATCGCCAAGGAAGCCCTcag                                              >  5:332050/1‑149 (MQ=255)
                                                                                                             cgcAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTca                                               >  8:196918/1‑139 (MQ=255)
                                                                                                             cgcAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTca                                               <  7:196918/139‑1 (MQ=255)
                                                                                                                aGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTgg                                  <  3:233341/149‑1 (MQ=255)
                                                                                                                 gTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGc                                 >  8:127401/1‑149 (MQ=255)
                                                                                                                    ccGAGACCTACGTCACCCACCCGCTGCCCGTGGCCAGCATCCTCGCCGACATCCACAAGGACCATCACAGCCTGATGGCGGCAATGCTGCACGATGTGATCGAAGACACAGCAATCGCCAAGGAAGCCCTCAGCCAGCAGtttccc                                 >  2:84293/2‑143 (MQ=255)
                                                                                                                         cccTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACTTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGTAGACTGGGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGt                         <  1:84293/149‑1 (MQ=255)
                                                                                                                            tACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGc                      >  7:6266/1‑149 (MQ=255)
                                                                                                                                tGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGaa                  >  3:203175/1‑149 (MQ=255)
                                                                                                                                 gACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAAt                 >  3:334155/1‑149 (MQ=255)
                                                                                                                                   cccacccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTg               <  1:135832/149‑1 (MQ=255)
                                                                                                                                    ccacccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTgg              >  3:168605/1‑149 (MQ=255)
                                                                                                                                    ccacccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTgg              >  4:254808/1‑149 (MQ=255)
                                                                                                                                    ccacccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTgg              >  4:297178/1‑149 (MQ=255)
                                                                                                                                     cacccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGt             <  8:6266/149‑1 (MQ=255)
                                                                                                                                          gCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGc                                                             >  6:241304/1‑96 (MQ=255)
                                                                                                                                          gCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGc                                                             <  5:241304/96‑1 (MQ=255)
                                                                                                                                          gCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATg        >  2:270200/1‑149 (MQ=255)
                                                                                                                                               ccGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTc   <  7:127401/149‑1 (MQ=255)
                                                                                                                                                cGTGGCCAGCATCCTCGCCGACTTGCTCATGGACCATCAGACCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGGATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCa  <  7:234787/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TATAGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCA  >  NC_002947/6049147‑6049439

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