Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,049,293 T→C 6.3% V52A (GTG→GCG) ‡ 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,2930TC 6.3% 234.3 / 4.5 79V52A (GTG→GCG) ‡spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base C (2/3);  total (39/40)
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.96e-01

ATAGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGC  >  NC_002947/6049148‑6049441
                                                                                                                                                 |                                                                                                                                                    
ataGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGc                                                                                                                                                   >  2:72456/1‑149 (MQ=255)
   gAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCa                                                                                                                                                 <  3:83851/148‑1 (MQ=255)
    aaGCCTTGGCCGAACGGCTTTCGCCCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGAGTGCCTATTTCGAAGCAGAACAGTCCCGCGATGGGCAACGCCCCCGCAGAGGCGAGCCCTAAGTGCCCCACCCGCTGGCCGTGGccccc                                                                                                                                               >  4:194139/1‑146 (MQ=255)
    aaGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAg                                                                                                                                                >  3:170476/1‑148 (MQ=255)
      gCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCAt                                                                                                                                             <  5:94520/149‑1 (MQ=255)
        cTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgggt                                                                                                                                           <  5:266807/149‑5 (MQ=255)
        cTTGGCATAACGTCTTTCGACATATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCTAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgggt                                                                                                                                           <  1:54863/149‑5 (MQ=255)
             ccGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGGGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGcc                                                                                                                                      <  4:219913/149‑1 (MQ=255)
                   ggCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATg                                                                                                                                <  2:243162/149‑1 (MQ=255)
                      tttCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGca                                                                                                                              >  3:207307/1‑148 (MQ=255)
                      tttCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGca                                                                                                                              >  4:92332/1‑148 (MQ=255)
                          gACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGGCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATgg                                                                                                                         <  4:207307/149‑1 (MQ=255)
                          gACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATgg                                                                                                                         <  2:73378/149‑1 (MQ=255)
                          gACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATgg                                                                                                                         <  3:220010/149‑1 (MQ=255)
                          gACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATgg                                                                                                                         <  3:92332/149‑1 (MQ=255)
                          gACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATg                                                                                                                          <  4:212192/148‑1 (MQ=255)
                           aCCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGa                                                                                                                        >  5:330957/1‑149 (MQ=255)
                            ccTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGAc                                                                                                                       <  4:75205/149‑1 (MQ=255)
                                  ttGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCAt                                                                                                                    >  6:68497/1‑146 (MQ=255)
                                     gCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGc                                                                                                                                               <  5:325174/116‑1 (MQ=255)
                                     gCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGc                                                                                                              >  1:217782/1‑149 (MQ=255)
                                       cccGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCt                                                                                                            >  3:5872/1‑149 (MQ=255)
                                       cccGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCt                                                                                                            >  4:203649/1‑149 (MQ=255)
                                       cccGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCt                                                                                                            >  7:260958/1‑149 (MQ=255)
                                              aGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATggcg                                                                                                     <  8:187654/149‑1 (MQ=255)
                                               ggTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATggcgg                                                                                                    <  3:136691/149‑1 (MQ=255)
                                                  cAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCa                                                                                                 >  4:236240/1‑149 (MQ=255)
                                                   aaCCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCAt                                                                                                >  8:293682/1‑149 (MQ=255)
                                                        ggTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCAt                                                                                                                                             >  3:256130/1‑99 (MQ=255)
                                                        ggTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCAt                                                                                                                                             <  4:256130/99‑1 (MQ=255)
                                                        ggTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCAtgctg                                                                                            >  2:206652/1‑148 (MQ=255)
                                                        ggTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCAtgctg                                                                                            >  3:78111/1‑148 (MQ=255)
                                                          ttCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCAc                                                                                         >  6:44362/1‑149 (MQ=255)
                                                             ggCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATCGACCATCAGAGCCTGATGGCGGCCATGCTGCACGa                                                                                       <  4:78111/148‑1 (MQ=255)
                                                                   ccTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc                                                                                <  7:293682/149‑1 (MQ=255)
                                                                   ccTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc                                                                                >  7:60941/1‑149 (MQ=255)
                                                                   ccTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc                                                                                <  8:60941/149‑1 (MQ=255)
                                                                   ccTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc                                                                                <  3:203649/149‑1 (MQ=255)
                                                                    cTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc                                                                                >  4:168467/1‑148 (MQ=255)
                                                                    cTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc                                                                                >  6:173096/1‑148 (MQ=255)
                                                                       tttCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGaa                                                                             <  8:260958/148‑1 (MQ=255)
                                                                          cTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtg                                                                                     >  7:37003/1‑137 (MQ=255)
                                                                          cTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtg                                                                                     <  8:37003/137‑1 (MQ=255)
                                                                             cGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGacac                                                                        <  5:173096/147‑1 (MQ=255)
                                                                               ccGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCgg                                                                     >  4:175963/1‑148 (MQ=255)
                                                                                 gAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCAt                                                                  >  2:62462/1‑149 (MQ=255)
                                                                                  aaCAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATc                                                                 >  5:231410/1‑149 (MQ=255)
                                                                                  aaCAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATc                                                                 <  5:44362/149‑1 (MQ=255)
                                                                                  aaCAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCATCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGCCCATGCTGCACGATGTGATCGAAGACACCGGCATc                                                                 >  2:100958/1‑149 (MQ=255)
                                                                                    cAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGc                                                               <  1:62462/149‑1 (MQ=255)
                                                                                       gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATc                                                                                                                   <  3:155828/94‑1 (MQ=255)
                                                                                       gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATc                                                                                                                   >  4:155828/1‑94 (MQ=255)
                                                                                       gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGcc                                                              >  6:11510/1‑147 (MQ=255)
                                                                                       gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGcc                                                              >  4:11903/1‑147 (MQ=255)
                                                                                       gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCaa                                                            <  6:231410/149‑1 (MQ=255)
                                                                                            cGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACCATGTGATCGAAGACACCGGCATCGCCAAGGAAg                                                       <  1:206652/149‑1 (MQ=255)
                                                                                              aTGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGcc                                                     <  3:175963/149‑1 (MQ=255)
                                                                                                  gCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCt                                                   >  3:162470/1‑147 (MQ=255)
                                                                                                  gCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGACGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCt                                                   >  7:6442/1‑147 (MQ=255)
                                                                                                   cAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTcag                                                >  8:70076/1‑149 (MQ=255)
                                                                                                    aaCGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTcagc                                               <  5:11510/149‑1 (MQ=255)
                                                                                                             gcAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGttt                                      >  7:287399/1‑149 (MQ=255)
                                                                                                              cAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTg                                     <  3:168467/149‑1 (MQ=255)
                                                                                                              cAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTg                                     <  8:287399/149‑1 (MQ=255)
                                                                                                                   gCGAGCCCTACGTGACCCCCCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGCTGTGATCGAAGACACCGGCAACGCAAATGCAGCCCTCAGCCAGCAGTGTGGc                                   >  8:228762/1‑146 (MQ=255)
                                                                                                                        cccTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGt                           >  4:195856/1‑149 (MQ=255)
                                                                                                                         ccTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACGCCGGCATCGCCAAGGAAGCCCTca                                                 <  8:295369/126‑1 (MQ=255)
                                                                                                                         ccTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACGCCGGCATCGCCAAGGAAGCCCTca                                                 >  7:295369/1‑126 (MQ=255)
                                                                                                                          cTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTgg                         <  4:162470/149‑1 (MQ=255)
                                                                                                                            aCGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGcc                       >  7:180075/1‑149 (MQ=255)
                                                                                                                             cGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGc                        <  3:11903/147‑1 (MQ=255)
                                                                                                                             cGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCg                      <  3:236240/149‑1 (MQ=255)
                                                                                                                               tGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGaa                    <  3:61963/149‑1 (MQ=255)
                                                                                                                                 acccacccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGTGAGACCGTGGCCGAAt                   <  4:170476/148‑1 (MQ=255)
                                                                                                                                    cacccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGt               >  5:137766/1‑149 (MQ=255)
                                                                                                                                      gatgCTGGCCACGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAAt                   <  8:146826/140‑1 (MQ=255)
                                                                                                                                      gatgCTGGCCACGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAAt                   >  7:146826/4‑143 (MQ=255)
                                                                                                                                      gatgCTGGCCACGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTc              >  4:225426/4‑148 (MQ=255)
                                                                                                                                          cTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATgg         <  2:18696/149‑1 (MQ=255)
                                                                                                                                                 tGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGACGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGc  <  8:6442/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ATAGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGC  >  NC_002947/6049148‑6049441

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