Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,246,284 C→T 7.3% 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,2840CT 7.3% 319.9 / 11.0 110intergenic (+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 (45/57);  new base T (6/2);  total (51/59)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 111 total aligned reads displayed.

GAAGTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATG  >  NC_002947/2246135‑2246431
                                                                                                                                                       |                                                                                                                                                    
gAAGTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTcc                                                                                                                                                      >  1:2288201/1‑150 (MQ=255)
         gccgccGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGcc                                                                                                                                             <  2:1217991/150‑1 (MQ=255)
                gCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTc                                                                                                                                      >  2:88145/1‑150 (MQ=255)
                gCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGATGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTc                                                                                                                                      <  2:1418570/150‑1 (MQ=255)
                       cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGcc                                                                                                                                             >  1:1786830/1‑136 (MQ=255)
                       cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGcc                                                                                                                                             <  2:1786830/136‑1 (MQ=255)
                       cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCtttt                                                                                                                                        >  1:1884655/1‑141 (MQ=255)
                       cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCtttt                                                                                                                                        <  2:1884655/141‑1 (MQ=255)
                       cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCtt                                                                                                                                          >  2:307125/1‑139 (MQ=255)
                       cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCtt                                                                                                                                          <  1:307125/139‑1 (MQ=255)
                       cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTtc                                                                                                                               >  2:1108629/1‑150 (MQ=255)
                        aaTGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCCtttttttcgt                                                                                                                                  <  2:744576/146‑6 (MQ=255)
                        aaTGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTtct                                                                                                                              >  2:1761998/1‑150 (MQ=255)
                        aaTGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTtct                                                                                                                              >  2:847046/1‑150 (MQ=255)
                        aaTGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATAGC‑GC‑GGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCCtttttttcgt                                                                                                                                  >  1:744576/1‑141 (MQ=39)
                            gCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑Cg                                                                                                                         >  1:76287/1‑150 (MQ=255)
                            gCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑Cg                                                                                                                         >  1:2375603/1‑150 (MQ=255)
                              cTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCa                                                                                                                       >  2:519208/1‑150 (MQ=255)
                              cTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCa                                                                                                                       <  2:2288201/150‑1 (MQ=255)
                                 cTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTc                                                                                                                    <  1:2067536/150‑1 (MQ=255)
                                   cAGCATCGACCGGGCGCGCATCGCGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTccc                                                                                                                  <  1:265862/150‑1 (MQ=255)
                                    aGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAATACGAGCGATG‑CCATAGCTA‑ACTTCCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCGTTTCA‑CGCAGTCCcg                                                                                                                 >  1:1994750/1‑150 (MQ=255)
                                      cATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑c                                                                                                                          >  2:2171439/1‑139 (MQ=255)
                                       aTCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGt                                                                                                              >  2:2417412/1‑150 (MQ=255)
                                        tCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTa                                                                                                             <  2:76287/150‑1 (MQ=255)
                                         cGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAg                                                                                                            <  1:2417412/150‑1 (MQ=255)
                                           aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGt                                                                                                          <  2:2207480/150‑1 (MQ=255)
                                                      aTCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTc                                                                                               >  2:151892/1‑150 (MQ=255)
                                                        cGTGCCTGAAGAAGCCACGGGCGAGCGCGGGCGTAGCAAGATCGACAAACCAGCGGTGCACT‑TTCAACGC‑TTGAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTATGTGCGCGTCCGTCAt                                                                                             <  1:1224016/150‑1 (MQ=255)
                                                             cTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGAt                                                                                        <  2:2375603/150‑1 (MQ=255)
                                                              tGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCGCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGACCGTCATGCGATc                                                                                       >  2:1300436/1‑150 (MQ=255)
                                                                    agCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCATCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAAt                                                                                 <  2:1487886/150‑1 (MQ=255)
                                                                    agCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAAt                                                                                 <  1:460918/150‑1 (MQ=255)
                                                                    agCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAAt                                                                                 <  1:519208/150‑1 (MQ=255)
                                                                      ccaaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAATGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCt                                                                               <  1:948006/150‑1 (MQ=255)
                                                                       caaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTg                                                                              >  1:115039/1‑150 (MQ=255)
                                                                          ggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGCTTATTTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCg                                                                           >  1:747502/1‑150 (MQ=255)
                                                                           gccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTAGTGCTAGTTC‑CCCAGGCCCGCGTAGGTACGCGGCCGTCATGCGATCAGTAATCTGCCgg                                                                          >  2:1244277/1‑150 (MQ=255)
                                                                            ccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgc                                                                         <  2:931061/150‑1 (MQ=255)
                                                                             caagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGc                                                                                           <  2:2442590/131‑1 (MQ=255)
                                                                             caagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGc                                                                                           >  1:2442590/1‑131 (MQ=255)
                                                                                   cATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCa                                                                  >  2:2159761/1‑150 (MQ=255)
                                                                                    aTTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAg                                                                 >  2:276705/1‑150 (MQ=255)
                                                                                     ttCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCGGCGGAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGc                                                                >  1:671229/1‑150 (MQ=255)
                                                                                     ttCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGc                                                                >  2:2025398/1‑150 (MQ=255)
                                                                                       ccGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCaa                                                              <  1:276705/150‑1 (MQ=255)
                                                                                        cGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAAc                                                             <  1:2025398/150‑1 (MQ=255)
                                                                                         gATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACa                                                            >  1:884328/1‑150 (MQ=255)
                                                                                            cAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCt                                                         <  1:2132048/150‑1 (MQ=255)
                                                                                              aTGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGc                                                       <  1:151892/150‑1 (MQ=255)
                                                                                              aTGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGc                                                       >  2:1409490/1‑150 (MQ=255)
                                                                                               tGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGcc                                                      <  2:107743/150‑1 (MQ=255)
                                                                                                ggCGACAATCCAGCTCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCg                                                     <  1:110828/150‑1 (MQ=255)
                                                                                                   gAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGaca                                                  >  2:256176/1‑148 (MQ=255)
                                                                                                   gAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGaaa                                                  <  1:1300436/150‑1 (MQ=255)
                                                                                                   gAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGaaa                                                  >  1:2114181/1‑150 (MQ=255)
                                                                                                      aaaCCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTc                                               <  1:88145/150‑1 (MQ=255)
                                                                                                        aCCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcg                                             <  1:2059245/150‑1 (MQ=255)
                                                                                                          cAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggg                                           >  1:2040151/1‑149 (MQ=255)
                                                                                                              gcTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacg                                       >  1:1036935/1‑150 (MQ=255)
                                                                                                              gcTGAACTGAACGA‑GC‑CGGAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGAGGTGCGTATAGCA‑TCCCGCGTCGGTGCGCGGCCGTCATGAGCTCAGTAAACTGCCGAGTTCATTAGCAACATCTGCCGAAACTCcggtaacg                                       <  1:1751835/150‑1 (MQ=255)
                                                                                                               cTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacgg                                      >  2:1751835/1‑150 (MQ=255)
                                                                                                                tGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacggt                                     >  1:961881/1‑150 (MQ=255)
                                                                                                                  aaCTGATCGA‑GC‑TGCAACGAAAAAAAGTGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCGCCAGCAACATCTGCCCAAACTCcggtaacggtac                                   <  2:271795/150‑2 (MQ=255)
                                                                                                                   aCTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGCTAACGGTAAt                                  <  2:1036935/150‑1 (MQ=255)
                                                                                                                       aTCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGcc                              <  2:884328/150‑1 (MQ=255)
                                                                                                                        tCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGc                                                       <  2:943428/124‑1 (MQ=255)
                                                                                                                        tCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGc                                                       >  1:943428/1‑124 (MQ=255)
                                                                                                                              c‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGGCTTCGCGTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCtgt                        <  1:61389/150‑1 (MQ=255)
                                                                                                                                tGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCtgtg                       >  1:1146901/1‑150 (MQ=255)
                                                                                                                                  cAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATg                                                                                            <  2:1659787/79‑1 (MQ=255)
                                                                                                                                  cAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATg                                                                                            >  1:1659787/1‑79 (MQ=255)
                                                                                                                                  cAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCa                     <  2:2081626/150‑1 (MQ=255)
                                                                                                                                  cAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCa                     <  2:1661753/150‑1 (MQ=255)
                                                                                                                                   aaCGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGc                                                                                                                                  >  2:76351/1‑41 (MQ=255)
                                                                                                                                   aaCGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGc                                                                                                                                  <  1:76351/41‑1 (MQ=255)
                                                                                                                                   aaCGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTTCAg                    >  1:1797689/1‑150 (MQ=255)
                                                                                                                                   aaCGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAg                    >  2:1804384/1‑150 (MQ=255)
                                                                                                                                       aaaaaaaGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAAt                                  <  2:2150513/132‑1 (MQ=255)
                                                                                                                                       aaaaaaaGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAAt                                  >  1:2150513/1‑132 (MQ=255)
                                                                                                                                       aaaaaaaGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATaa                <  2:115039/150‑1 (MQ=255)
                                                                                                                                          aaaaGGGCACCTCTGGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAAccc             >  1:614909/1‑150 (MQ=255)
                                                                                                                                          aaaaGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTACCGAAAATGTTCTGCTGTGCAGATAAccc             <  2:2040151/150‑1 (MQ=255)
                                                                                                                                           aaaGGGCACCTCTGGAGGTGCCTTTTTTCGTGCCTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCt            >  1:1000125/1‑150 (MQ=255)
                                                                                                                                             aGGGCACCTCTGGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACTGTAATGGCCGGCTGTGCAGATAACCCTgg          >  2:2184639/1‑150 (MQ=255)
                                                                                                                                             aGGGCACCTCTGGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTgg          >  1:603382/1‑150 (MQ=255)
                                                                                                                                             aGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTgg          <  1:1098185/150‑1 (MQ=255)
                                                                                                                                              gggCACCTCTGGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGt         >  1:1455236/1‑150 (MQ=255)
                                                                                                                                               ggCACCTCTGGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTa        <  2:603382/150‑1 (MQ=255)
                                                                                                                                               ggCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTa        <  2:728460/150‑1 (MQ=255)
                                                                                                                                                   ccTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCCGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGa    <  1:1108629/150‑1 (MQ=255)
                                                                                                                                                   ccTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGa    <  1:1761998/150‑1 (MQ=255)
                                                                                                                                                    cTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGc                                                                <  1:1093409/89‑1 (MQ=255)
                                                                                                                                                    cTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGc                                                                >  2:1093409/1‑89 (MQ=255)
                                                                                                                                                    cTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGAt   <  2:1993339/150‑1 (MQ=255)
                                                                                                                                                    cTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGAt   <  2:1481993/150‑1 (MQ=255)
                                                                                                                                                    cTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGAt   <  1:2159761/150‑1 (MQ=255)
                                                                                                                                                    cTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGAt   >  1:2184419/1‑150 (MQ=255)
                                                                                                                                                     tCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTgg          >  1:1643573/1‑142 (MQ=255)
                                                                                                                                                     tCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATg  <  2:1797689/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GAAGTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGA‑GC‑TGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCA‑CGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATG  >  NC_002947/2246135‑2246431

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