Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,923,250 C→G 5.6% intergenic (‑35/+83) 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,2500CG 5.6% 666.9 / 22.0 198intergenic (‑35/+83)PP_1718/prcEAL domain‑containing diguanylate phosphodiesterase/tail‑specific protease Prc
Reads supporting (aligned to +/- strand):  ref base C (107/80);  new base G (3/8);  total (110/88)
Fisher's exact test for biased strand distribution p-value = 6.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 202 total aligned reads displayed.

CCTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCGCGAGG  >  NC_002947/1923101‑1923399
                                                                                                                                                     |                                                                                                                                                     
ccTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGc                                                                                                                                                       >  2:3826938/1‑150 (MQ=255)
ccTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGc                                                                                                                                                       <  2:4175611/150‑1 (MQ=255)
ccTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACCGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGc                                                                                                                                                       <  1:3969160/150‑1 (MQ=255)
 cTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCa                                                                                                                                                      <  2:3511459/150‑1 (MQ=255)
 cTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCa                                                                                                                                                      >  2:732648/1‑150 (MQ=255)
      ggCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGccc                                                                                                                                                 <  1:5258202/150‑1 (MQ=255)
         cGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGc                                                                                                                                              <  2:1824708/150‑1 (MQ=255)
          gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGcc                                                                                                                                             <  2:4894428/150‑1 (MQ=255)
          gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGcc                                                                                                                                             <  1:954766/150‑1 (MQ=255)
          gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGcc                                                                                                                                             >  2:4279828/1‑150 (MQ=255)
          gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGcc                                                                                                                                             >  2:837455/1‑150 (MQ=255)
          gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGcc                                                                                                                                             >  1:1853157/1‑150 (MQ=255)
           aaTACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGccc                                                                                                                                            <  1:2410813/150‑1 (MQ=255)
              aCCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCa                                                                                                                                         <  2:2994301/150‑1 (MQ=255)
                 cTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCatat                                                                                                                                      >  2:4066289/1‑150 (MQ=255)
                     gTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGCCCTTACTTTTCa                                                                                                                                  <  1:3695137/150‑1 (MQ=39)
                     gTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGCCCTTACTTTTCa                                                                                                                                  <  2:3780101/150‑1 (MQ=39)
                      ttCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCAc                                                                                                                                 <  1:2131923/150‑1 (MQ=255)
                       tCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACg                                                                                                                                >  2:2547119/1‑150 (MQ=255)
                             aGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCa                                                                                                                          >  1:2891273/1‑150 (MQ=255)
                              ggCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCAc                                                                                                                         >  1:5356053/1‑150 (MQ=255)
                                 acaTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGAt                                                                                                                      <  2:5612766/150‑1 (MQ=255)
                                    tGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGAc                                                                                                                   <  2:2471202/150‑1 (MQ=255)
                                      aTTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTg                                                                                                                 <  1:2567565/150‑1 (MQ=255)
                                      aTTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTg                                                                                                                 >  2:2687148/1‑150 (MQ=255)
                                         ggCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGttt                                                                                                              <  1:1642611/150‑1 (MQ=255)
                                         ggCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGttt                                                                                                              <  1:4066288/150‑1 (MQ=255)
                                            tGGAACAGACTGTGTATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgat                                                                                                           <  2:3724991/150‑1 (MQ=255)
                                            tGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgat                                                                                                           >  2:3611455/1‑150 (MQ=255)
                                            tGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgat                                                                                                           >  2:933566/1‑150 (MQ=255)
                                            tGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgat                                                                                                           >  2:2184436/1‑150 (MQ=255)
                                            tGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgat                                                                                                           >  2:426383/1‑150 (MQ=255)
                                            tGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgat                                                                                                           <  1:1202919/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                                                                                               >  1:5602946/1‑150 (MQ=255)
                                                            tACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACattatt                                                                                           >  2:1915830/1‑150 (MQ=255)
                                                            tACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACattatt                                                                                           >  2:5510991/1‑150 (MQ=255)
                                                              cTGCGCTGAGCCTGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGTATTCCTGAAAAGTATGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCt                                                                                         <  1:3794545/150‑1 (MQ=255)
                                                                gcgcTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGc                                                                                       <  2:313060/150‑1 (MQ=255)
                                                                    tGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCcatc                                                                                   >  2:119231/1‑150 (MQ=255)
                                                                       gCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCAtt                                                                                >  1:4385419/1‑150 (MQ=255)
                                                                        ccAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTg                                                                               >  1:792888/1‑150 (MQ=255)
                                                                          aGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCt                                                                             <  2:5028030/150‑1 (MQ=255)
                                                                               ggCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGt                                                                        <  1:4884787/150‑1 (MQ=255)
                                                                                  aCTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTAGTCTGCCATCATTGCTGAGGTGcc                                                                     >  2:4238989/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:3632774/1‑150 (MQ=255)
                                                                                           aCTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTActtctt                                                            >  1:4188122/1‑150 (MQ=255)
                                                                                               aGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGTTGCCTTACTTCTTggac                                                        >  2:3362206/1‑148 (MQ=255)
                                                                                                   gTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCt                                                    >  2:2287909/1‑150 (MQ=255)
                                                                                                    ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGa                                                                           >  2:4084562/1‑126 (MQ=255)
                                                                                                    ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCtt                                                   >  1:4489591/1‑150 (MQ=255)
                                                                                                    ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCtt                                                   <  2:2512706/150‑1 (MQ=255)
                                                                                                    ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCtt                                                   >  1:716918/1‑150 (MQ=255)
                                                                                                     tCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGTCTTACTTCTTGGCCACCTTg                                                  <  1:5594647/150‑1 (MQ=255)
                                                                                                     tCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTg                                                  <  2:2501467/150‑1 (MQ=255)
                                                                                                      cGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTgg                                                 <  2:2728049/150‑1 (MQ=255)
                                                                                                       ggCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGt                                                <  2:5168592/150‑1 (MQ=255)
                                                                                                          aaCGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCt                                             >  1:1779275/1‑150 (MQ=255)
                                                                                                          aaCGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCt                                             >  1:3501120/1‑150 (MQ=255)
                                                                                                               tCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCAc                                        >  2:2899475/1‑150 (MQ=255)
                                                                                                                  tgatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                     >  1:1108214/1‑150 (MQ=255)
                                                                                                                  tgatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                     <  1:1915830/150‑1 (MQ=255)
                                                                                                                  tgatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                     >  1:882191/1‑150 (MQ=255)
                                                                                                                  tgatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTACCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                     >  2:844583/1‑150 (MQ=255)
                                                                                                                   gatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAg                                    <  2:1658618/150‑1 (MQ=255)
                                                                                                                    atgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACGGAgg                                   >  2:3650116/1‑150 (MQ=255)
                                                                                                                      ggATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTa                                 <  2:2376905/150‑1 (MQ=255)
                                                                                                                           ccTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGtggccaa                            <  1:1254691/150‑7 (MQ=255)
                                                                                                                            cTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGAt                           >  1:4079082/1‑150 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                          >  2:4048187/1‑150 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                          <  2:3186716/150‑1 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                          <  2:1414688/150‑1 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                          >  2:4461303/1‑150 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                          >  1:3382898/1‑150 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                          >  1:4938525/1‑150 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                          >  2:2640239/1‑150 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                          >  2:4932452/1‑150 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTCGGTGCTGGCACTGAGGTAGTCGATc                          >  1:4061805/1‑150 (MQ=255)
                                                                                                                                  aGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGAt                     <  1:2290464/150‑1 (MQ=255)
                                                                                                                                        ggCCGGGCATGCAGATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGAcc               >  2:847289/1‑150 (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)
                                                                                                                                           cGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGt            >  2:3620783/1‑150 (MQ=255)
                                                                                                                                            gggCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGtt           >  1:3321289/1‑150 (MQ=255)
                                                                                                                                                aTCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCcg       >  2:1539779/1‑150 (MQ=255)
                                                                                                                                                aTCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCcg       >  2:5606539/1‑150 (MQ=255)
                                                                                                                                                   tGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCGCGa    <  2:2160151/150‑1 (MQ=255)
                                                                                                                                                    gCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGGGTCcgcggg   >  2:4323227/1‑148 (MQ=255)
                                                                                                                                                     cATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCGCGAgg  <  1:3620783/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCGCGAGG  >  NC_002947/1923101‑1923399

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: