Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,923,247 C→G 5.5% intergenic (‑32/+86) PP_1718 ← / ← prc EAL domain‑containing diguanylate phosphodiesterase/tail‑specific protease Prc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,923,2470CG 5.5% 673.9 / 22.0 200intergenic (‑32/+86)PP_1718/prcEAL domain‑containing diguanylate phosphodiesterase/tail‑specific protease Prc
Reads supporting (aligned to +/- strand):  ref base C (107/82);  new base G (3/8);  total (110/90)
Fisher's exact test for biased strand distribution p-value = 6.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

Only 100 of 201 total aligned reads displayed.

GCGCCTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCG  >  NC_002947/1923098‑1923394
                                                                                                                                                     |                                                                                                                                                   
 cgcCTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCt                                                                                                                                                    <  2:1297753/150‑1 (MQ=255)
    cTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCa                                                                                                                                                 >  2:732648/1‑150 (MQ=255)
         ggCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGccc                                                                                                                                            <  1:5258202/150‑1 (MQ=255)
             gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGcc                                                                                                                                        <  1:954766/150‑1 (MQ=255)
             gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGcc                                                                                                                                        <  2:4894428/150‑1 (MQ=255)
             gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGcc                                                                                                                                        >  2:837455/1‑150 (MQ=255)
             gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGcc                                                                                                                                        >  1:1853157/1‑150 (MQ=255)
             gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGcc                                                                                                                                        >  2:2290464/1‑150 (MQ=255)
              aaTACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGccc                                                                                                                                       <  1:2410813/150‑1 (MQ=255)
                   ccTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCata                                                                                                                                  >  2:4066288/1‑150 (MQ=255)
                   ccTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCata                                                                                                                                  >  2:3595000/1‑150 (MQ=255)
                   ccTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCata                                                                                                                                  >  2:2506973/1‑150 (MQ=255)
                        gTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGCCCTTACTTTTCa                                                                                                                             <  1:367687/150‑1 (MQ=39)
                        gTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGCCCTTACTTTTCa                                                                                                                             <  2:3011103/150‑1 (MQ=39)
                         ttCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCAc                                                                                                                            <  1:2131923/150‑1 (MQ=255)
                             cTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACa                                                                                                                        <  2:1549406/150‑1 (MQ=255)
                             cTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACa                                                                                                                        <  2:3641294/150‑1 (MQ=255)
                                aGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCa                                                                                                                     >  1:2891273/1‑150 (MQ=255)
                                 ggCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCAc                                                                                                                    >  1:5356053/1‑150 (MQ=255)
                                    acaTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGAt                                                                                                                 <  2:5612766/150‑1 (MQ=255)
                                      aTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGa                                                                                                               <  2:3772599/150‑1 (MQ=255)
                                      aTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGa                                                                                                               <  1:3765856/150‑1 (MQ=255)
                                      aTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGACAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGa                                                                                                               <  1:4309583/150‑1 (MQ=255)
                                       tGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGAc                                                                                                              <  2:2471128/150‑1 (MQ=255)
                                        gATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACt                                                                                                             >  2:26652/1‑150 (MQ=255)
                                         aTTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTg                                                                                                            <  1:2567565/150‑1 (MQ=255)
                                          ttGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGt                                                                                                           >  2:4162612/1‑150 (MQ=255)
                                            ggCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGttt                                                                                                         <  1:4066289/150‑1 (MQ=255)
                                            ggCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGttt                                                                                                         <  1:1642611/150‑1 (MQ=255)
                                               tgtaaCAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgat                                                                                                      <  2:4276074/147‑1 (MQ=255)
                                               tGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgat                                                                                                      <  1:1202919/150‑1 (MQ=255)
                                               tGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgat                                                                                                      >  2:426332/1‑150 (MQ=255)
                                               tGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgat                                                                                                      >  2:1490550/1‑150 (MQ=255)
                                               tGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgat                                                                                                      >  2:933566/1‑150 (MQ=255)
                                                 gAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgatga                                                                                                    <  2:1869411/150‑1 (MQ=255)
                                                  aaCAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGAc                                                                                                   <  1:5515690/150‑1 (MQ=255)
                                                   aCAGACTGTGACTACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACa                                                                                                  <  2:5456913/150‑1 (MQ=255)
                                                   aCAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACa                                                                                                  <  1:5178490/150‑1 (MQ=255)
                                                        cTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATca                                                                                             >  1:1414688/1‑150 (MQ=255)
                                                         tgtgAATACTGCGCTGAGCCAGGAGGGCACTCAAGGCACTGAGCTGTTCGGCACCGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATcac                                                                                            <  2:4844849/150‑1 (MQ=255)
                                                           tgAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACat                                                                                          >  2:4022889/1‑150 (MQ=255)
                                                           tgAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACat                                                                                          >  1:5602946/1‑150 (MQ=255)
                                                              aTACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACattat                                                                                       <  1:4430776/150‑1 (MQ=255)
                                                               tACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACattatt                                                                                      >  2:5510991/1‑150 (MQ=255)
                                                                     gcTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCca                                                                                <  2:2941660/150‑1 (MQ=255)
                                                                       tGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCcatc                                                                              >  2:119231/1‑150 (MQ=255)
                                                                       tGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCcatc                                                                              >  2:4179447/1‑150 (MQ=255)
                                                                           ccAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTg                                                                          >  1:792888/1‑150 (MQ=255)
                                                                             aGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCt                                                                        <  2:2011024/150‑1 (MQ=255)
                                                                             aGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCt                                                                        <  2:5028030/150‑1 (MQ=255)
                                                                                  ggCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGt                                                                   <  1:4884787/150‑1 (MQ=255)
                                                                                  ggCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGt                                                                   <  1:4003612/150‑1 (MQ=255)
                                                                                     aCTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGcc                                                                >  2:3142013/1‑150 (MQ=255)
                                                                                         aaGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTAc                                                            >  1:806494/1‑150 (MQ=255)
                                                                                          aGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCAGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTAct                                                           >  1:4558085/1‑150 (MQ=255)
                                                                                            gCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTActtc                                                         >  1:5098382/1‑150 (MQ=255)
                                                                                              aCTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTActtctt                                                       <  2:5542960/150‑1 (MQ=255)
                                                                                              aCTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTActtctt                                                       >  1:3632772/1‑150 (MQ=255)
                                                                                               cTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTg                                                      <  2:3687136/150‑1 (MQ=255)
                                                                                                       ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGa                                                                      <  1:4084562/126‑1 (MQ=255)
                                                                                                       ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCtt                                              >  1:4489591/1‑150 (MQ=255)
                                                                                                       ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCtt                                              <  2:4074067/150‑1 (MQ=255)
                                                                                                       ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCtt                                              >  1:716918/1‑150 (MQ=255)
                                                                                                        tCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGTCTTACTTCTTGGCCACCTTg                                             <  1:5594647/150‑1 (MQ=255)
                                                                                                        tCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTg                                             <  2:2501417/150‑1 (MQ=255)
                                                                                                          ggCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGt                                           <  2:5168592/150‑1 (MQ=255)
                                                                                                             aaCGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCt                                        >  1:3501120/1‑150 (MQ=255)
                                                                                                             aaCGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCt                                        >  1:1779275/1‑150 (MQ=255)
                                                                                                              aCGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTg                                       >  2:2599065/1‑150 (MQ=255)
                                                                                                                     tgatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                <  1:1915830/150‑1 (MQ=255)
                                                                                                                     tgatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                >  1:1108214/1‑150 (MQ=255)
                                                                                                                     tgatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                >  2:34481/1‑150 (MQ=255)
                                                                                                                     tgatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                >  1:882191/1‑150 (MQ=255)
                                                                                                                     tgatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTACCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                >  2:844583/1‑150 (MQ=255)
                                                                                                                     tgatgGATTCCTGAAAAGTAAGGGCCGGGCAACTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                >  2:254451/1‑150 (MQ=255)
                                                                                                                         ggATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTCAGTGCCAGCACCAAGGtg                            <  1:3860047/150‑2 (MQ=255)
                                                                                                                              ccTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGtggccaa                       <  1:1254691/150‑7 (MQ=255)
                                                                                                                                tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                     >  2:2269505/1‑150 (MQ=255)
                                                                                                                                tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                     >  1:3382898/1‑150 (MQ=255)
                                                                                                                                tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                     >  2:4932452/1‑150 (MQ=255)
                                                                                                                                tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                     >  2:4461303/1‑150 (MQ=255)
                                                                                                                                tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                     >  1:4938525/1‑150 (MQ=255)
                                                                                                                                tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                     <  2:2698708/150‑1 (MQ=255)
                                                                                                                                tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTCGGTGCTGGCACTGAGGTAGTCGATc                     >  1:4062095/1‑150 (MQ=255)
                                                                                                                                  aaaaGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCCCTGAGGTAGTCGATCAg                   >  2:3615363/1‑150 (MQ=255)
                                                                                                                                    aaGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGa                 <  2:4300074/150‑1 (MQ=255)
                                                                                                                                     aGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGAt                <  2:2891273/150‑1 (MQ=255)
                                                                                                                                     aGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGAt                <  1:2290464/150‑1 (MQ=255)
                                                                                                                                      atgaGGGCCGGGCATGCAGATGCCCGGCCCTTATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATg               >  2:3210556/4‑150 (MQ=255)
                                                                                                                                       tAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGc              >  2:3797364/1‑150 (MQ=255)
                                                                                                                                       tAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTAACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGc              >  2:1659276/1‑150 (MQ=255)
                                                                                                                                           ggCCGGGCATGCAGATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGAcc          >  2:847289/1‑150 (MQ=255)
                                                                                                                                           ggCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAg                               >  2:2166433/1‑129 (MQ=255)
                                                                                                                                           ggCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAg                               <  1:2166433/129‑1 (MQ=255)
                                                                                                                                           ggCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGt                          <  2:5332252/134‑1 (MQ=255)
                                                                                                                                           ggCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGAcc          >  1:659969/1‑150 (MQ=255)
                                                                                                                                            gCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCg         >  2:4648830/1‑150 (MQ=255)
                                                                                                                                              cGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGt       >  2:4578512/1‑150 (MQ=255)
                                                                                                                                               gggCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGtt      >  1:3321289/1‑150 (MQ=255)
                                                                                                                                                   aTCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCcg  >  2:5606539/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCGCCTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCG  >  NC_002947/1923098‑1923394

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: