Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,246,241 C→T 7.9% intergenic (+30/+61) ibpA → / ← PP_1983 small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,246,2410CT 7.9% 179.3 / 4.4 63intergenic (+30/+61)ibpA/PP_1983small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base T (2/3);  total (30/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.06e-01

GGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATG  >  NC_002947/2246100‑2246388
                                                                                                                                              |                                                                                                                                                    
ggCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTg                                                                                                                                              <  3:215221/149‑1 (MQ=255)
  gcggcgGCCTGGCCATTGGGCTGGTCATCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCTATGCCAAGCGCATTCCTATCAATGGCGAAAAACAAGCGCTGAACTGATCTAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGc                                                                                                                                             <  4:274354/144‑1 (MQ=255)
   cgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGccc                                                                                                                                           <  8:375847/149‑1 (MQ=255)
   cgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGccc                                                                                                                                           >  2:313727/1‑149 (MQ=255)
   cgccgGCCTGGCCAATGGCCGGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGccc                                                                                                                                           <  8:278528/149‑1 (MQ=255)
     ccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                                 <  1:403547/141‑1 (MQ=255)
     ccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                                 >  2:403547/1‑141 (MQ=255)
      cgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCCtt                                                                                                                                         <  5:238151/148‑1 (MQ=255)
      cgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCCtt                                                                                                                                         >  6:238151/1‑148 (MQ=255)
         ccTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTc                                                                                                                                     >  4:99648/1‑149 (MQ=255)
          cTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCGCCTCCAGCGGTGCGCTTTTTCg                                                                                                                                    >  4:324700/1‑149 (MQ=255)
          cTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCg                                                                                                                                    >  4:356524/1‑149 (MQ=255)
             gCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACG‑‑AAAAAGGGCACCTCTGGAGGTGCCCtttttttcgtt                                                                                                                                <  8:73110/148‑7 (MQ=38)
               cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCtt                                                                                                                               <  6:250811/149‑1 (MQ=255)
                  tGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTtctc                                                                                                                            <  5:47750/149‑1 (MQ=255)
                    gCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCAc                                                                                                                          <  6:255136/149‑1 (MQ=255)
                        gCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAg                                                                                                                      >  2:413043/1‑149 (MQ=255)
                            aGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTc                                                                                                                    <  3:356524/147‑1 (MQ=255)
                                tCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGt                                                                                                              <  1:321858/149‑1 (MQ=255)
                                  gACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTa                                                                                                             <  6:271096/148‑1 (MQ=255)
                                  gACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTa                                                                                                             <  1:413043/148‑1 (MQ=255)
                                   aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCATCCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAg                                                                                                            <  1:313727/149‑1 (MQ=255)
                                    ccTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGt                                                                                                          >  7:230293/1‑149 (MQ=255)
                                      tgctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTATTCGTGCTTCTCACGCAGTCCCGCGTAGGTgc                                                                                                        >  3:174276/1‑149 (MQ=255)
                                       gctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCtttt                                                                                                                                       >  2:231387/1‑117 (MQ=255)
                                       gctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCtttt                                                                                                                                       <  1:231387/117‑1 (MQ=255)
                                        ctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTgcgc                                                                                                      >  7:40948/1‑149 (MQ=255)
                                            gcATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGTGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGcc                                                                                                  <  2:27513/149‑1 (MQ=255)
                                                  tGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATg                                                                                            <  5:430936/149‑1 (MQ=255)
                                                    ccTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATTGCGAAAAAACAGCGCTGAACTGATCGAGCTGCAACGAAAAAAATGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACTCAGTCCCGCGTAGGTGCGAGGCAGTCAttc                                                                                           >  4:252887/1‑146 (MQ=255)
                                                     cTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGa                                                                                         >  1:118715/1‑149 (MQ=255)
                                                        aagaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGAt                                                                                        >  3:360430/1‑147 (MQ=255)
                                                         agaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCa                                                                                      <  3:99648/148‑1 (MQ=255)
                                                         agaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAg                                                                                     <  8:40948/149‑1 (MQ=255)
                                                             gccaaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAAt                                                                                 >  2:128833/1‑149 (MQ=255)
                                                              ccaaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATc                                                                                <  4:360430/149‑1 (MQ=255)
                                                               caaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGcgcg                                                                                                     >  3:287240/1‑127 (MQ=255)
                                                               caaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGcgcg                                                                                                     <  4:287240/127‑1 (MQ=255)
                                                               caaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCt                                                                               >  3:349432/1‑149 (MQ=255)
                                                                  ggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACAAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGcc                                                                            >  3:199628/1‑149 (MQ=255)
                                                                      aagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgcg                                                                        >  6:68719/1‑149 (MQ=255)
                                                                             ttCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAg                                                                 >  4:123928/1‑149 (MQ=255)
                                                                             ttCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAg                                                                 >  8:45641/1‑149 (MQ=255)
                                                                                    cAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATc                                                          >  3:216174/1‑149 (MQ=255)
                                                                                       tGGCGAAAAACCAGCGCTGAACTGATCGAGCTTCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGc                                                       <  8:45361/149‑1 (MQ=255)
                                                                                           gAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGa                                                    >  3:37193/1‑148 (MQ=255)
                                                                                                    gcgcTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggt                                           >  2:383123/1‑148 (MQ=255)
                                                                                                       cTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacg                                       >  1:418508/1‑149 (MQ=255)
                                                                                                          aaCTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacggt                                     <  8:230293/148‑1 (MQ=255)
                                                                                                               aTCGAGCTTCAACGAAAAAAACGGCACGTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATg                                 <  4:174276/147‑1 (MQ=255)
                                                                                                               aTCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGc                               >  2:102856/1‑149 (MQ=255)
                                                                                                               aTCGAGCTGCAACGAAAAAAACGCCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCTCTTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATg                                 <  3:252887/147‑1 (MQ=255)
                                                                                                                    gCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCt                          <  4:216174/149‑1 (MQ=255)
                                                                                                                       gCAACGAAAAAAAGGGCACCTCCAGAGGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCtgt                        >  1:19448/1‑149 (MQ=255)
                                                                                                                        caacaaAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCtgtg                       <  4:199628/148‑1 (MQ=255)
                                                                                                                        cAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGc                      <  2:118715/149‑1 (MQ=255)
                                                                                                                            gAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGa                   >  3:364675/1‑148 (MQ=255)
                                                                                                                            acgaaaaaGGGCACCTCTGGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGa                   >  4:48762/4‑148 (MQ=255)
                                                                                                                              aaaaaaGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATaa                <  4:349432/149‑1 (MQ=255)
                                                                                                                                   agcgCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTg           <  3:324700/146‑1 (MQ=255)
                                                                                                                                   aGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTg           <  1:128833/149‑1 (MQ=255)
                                                                                                                                    gggCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTgg          >  1:306053/1‑149 (MQ=255)
                                                                                                                                           cTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGa    <  4:364675/149‑1 (MQ=255)
                                                                                                                                              cAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATg  <  3:123928/148‑1 (MQ=255)
                                                                                                                                              cAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATg  <  7:45641/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATG  >  NC_002947/2246100‑2246388

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: