Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,339,148 C→T 7.4% G618S (GGT→AGT)  gmr ← cyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of Rnase II stability

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,339,1480CT 7.4% 828.2 / 36.0 241G618S (GGT→AGT) gmrcyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of Rnase II stability
Reads supporting (aligned to +/- strand):  ref base C (115/108);  new base T (6/12);  total (121/120)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

ATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAG  >  CP009273/1339000‑1339295
                                                                                                                                                    |                                                                                                                                                   
atatatTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCt                                                                                                                                                    >  2:1083745/1‑150 (MQ=255)
   tatTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCAACACCTTCg                                                                                                                                                 >  1:1284327/1‑150 (MQ=255)
    atTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCgg                                                                                                                                                <  2:3865658/150‑1 (MQ=255)
    atTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCgg                                                                                                                                                >  1:1915197/1‑150 (MQ=255)
      tCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCg                                                                                                                                              >  1:178119/1‑150 (MQ=255)
            aCTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCAcc                                                                                                                                        >  1:438203/1‑150 (MQ=255)
            aCTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCAcc                                                                                                                                        >  2:4688094/1‑150 (MQ=255)
              ttATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTg                                                                                                                                      >  2:1811756/1‑150 (MQ=255)
              ttATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTg                                                                                                                                      >  2:4197077/1‑150 (MQ=255)
               tATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGa                                                                                                                                     <  2:585671/150‑1 (MQ=255)
                  gcgcgcTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGa                                                                                                                                  <  1:3531789/150‑1 (MQ=255)
                  gcgcgcTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGa                                                                                                                                  <  1:3531911/150‑1 (MQ=255)
                   cgcgcTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGAt                                                                                                                                 >  2:777187/1‑150 (MQ=255)
                    gcgcTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGAtt                                                                                                                                >  1:3665868/1‑150 (MQ=255)
                     cgcTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTc                                                                                                                               <  2:4626043/150‑1 (MQ=255)
                       cTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCaa                                                                                                                             >  2:2670180/1‑150 (MQ=255)
                              tAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTgg                                                                                                                      <  1:117030/150‑1 (MQ=255)
                               aGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTggg                                                                                                                     >  1:3929425/1‑150 (MQ=255)
                                  gTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCca                                                                                                                  >  2:652122/1‑150 (MQ=255)
                                  gTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCca                                                                                                                  >  2:2165769/1‑150 (MQ=255)
                                    ttATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCcaca                                                                                                                <  2:1433376/150‑1 (MQ=255)
                                    ttATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCcaca                                                                                                                <  1:1732838/150‑1 (MQ=255)
                                    ttATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCcaca                                                                                                                <  2:200230/150‑1 (MQ=255)
                                      ataCCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGc                                                                                                              <  2:2848773/150‑1 (MQ=255)
                                         ccAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGcgac                                                                                                           >  1:4592302/1‑150 (MQ=255)
                                           aGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACTTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGcgacga                                                                                                         <  2:198175/150‑1 (MQ=255)
                                             cGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATc                                                                                                       <  1:266244/150‑1 (MQ=255)
                                              gTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCg                                                                                                      >  1:1091084/1‑150 (MQ=255)
                                                 cGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGccc                                                                                                   <  1:1018060/150‑1 (MQ=255)
                                                 cGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGccc                                                                                                   <  1:2702040/150‑1 (MQ=255)
                                                   aaGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCgg                                                                                                 >  1:3978976/1‑150 (MQ=255)
                                                     ggCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGAc                                                                                               <  2:2431635/150‑1 (MQ=255)
                                                      gCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACTTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGAcc                                                                                              <  1:755264/150‑1 (MQ=255)
                                                       cGACGGCGGGCATCGGTTTGGCAAACAACCATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCGGAAGATTCCATGCCTGGGCCCCAGCGACGATCGCCCGGACCa                                                                                             <  1:1543337/150‑1 (MQ=255)
                                                       cGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCa                                                                                             <  1:3263020/150‑1 (MQ=255)
                                                       cGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCa                                                                                             >  2:3938275/1‑150 (MQ=255)
                                                          cggcggGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTg                                                                                          <  2:1653001/150‑1 (MQ=255)
                                                              gggCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTg                                                                                      >  1:461406/1‑150 (MQ=255)
                                                               ggCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGCTCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGc                                                                                     <  2:3749426/150‑1 (MQ=255)
                                                                   tCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACTTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCgaga                                                                                 <  2:2929131/150‑1 (MQ=255)
                                                                   tCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACTTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCgaga                                                                                 <  2:2470329/150‑1 (MQ=255)
                                                                    cGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGAc                                                                                >  2:3754708/1‑150 (MQ=255)
                                                                     ggTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACa                                                                               <  2:3274879/150‑1 (MQ=255)
                                                                      gTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAg                                                                              >  2:4426505/1‑150 (MQ=255)
                                                                       tttGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAgg                                                                             <  1:1650173/150‑1 (MQ=255)
                                                                       tttGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAgg                                                                             <  2:4533506/150‑1 (MQ=255)
                                                                       tttGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAgg                                                                             <  1:601865/150‑1 (MQ=255)
                                                                        ttGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGgt                                                                            <  2:2900592/150‑1 (MQ=255)
                                                                         tGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGgtt                                                                           <  1:1285631/150‑1 (MQ=255)
                                                                              aacaaaAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGttgtttg                                                                      <  2:184045/150‑1 (MQ=255)
                                                                              aacaaaAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGttgtttg                                                                      <  1:3550736/150‑1 (MQ=255)
                                                                                 aaaaaTCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGa                                                                   >  2:3167826/1‑150 (MQ=255)
                                                                                  aaaaTCCTTGCCGCTCAGTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACTTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGaa                                                                  >  2:3689119/1‑150 (MQ=255)
                                                                                   aaaTCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACTTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAat                                                                 >  1:3126120/1‑150 (MQ=255)
                                                                                   aaaTCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAat                                                                 >  2:3226421/1‑150 (MQ=255)
                                                                                     aTCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAatat                                                               >  2:1185358/1‑150 (MQ=255)
                                                                                      tCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATAtc                                                              >  2:910047/1‑150 (MQ=255)
                                                                                      tCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATAtc                                                              <  1:4711561/150‑1 (MQ=255)
                                                                                        cgtGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTGCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAAGGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATAtctc                                                            <  2:3170061/148‑1 (MQ=255)
                                                                                          tGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGa                                                          >  2:4729095/1‑150 (MQ=255)
                                                                                           gCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGaa                                                         <  2:2763340/150‑1 (MQ=255)
                                                                                              gCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCAACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACaa                                                      <  2:1284327/150‑1 (MQ=255)
                                                                                               cTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACaaa                                                     <  1:194355/150‑1 (MQ=255)
                                                                                                tCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACaaaa                                                    >  2:4032554/1‑150 (MQ=255)
                                                                                                tCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACaaaa                                                    <  1:3482261/150‑1 (MQ=255)
                                                                                                  aTTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAAc                                                  >  1:2391596/1‑150 (MQ=255)
                                                                                                  aTTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAAc                                                  <  1:2795704/150‑1 (MQ=255)
                                                                                                   ttGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAAcc                                                 >  1:476143/1‑150 (MQ=255)
                                                                                                    tGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCt                                                <  1:2006117/150‑1 (MQ=255)
                                                                                                    tGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCt                                                <  1:3414603/150‑1 (MQ=255)
                                                                                                    tGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCt                                                <  2:2083921/150‑1 (MQ=255)
                                                                                                     gATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTg                                               <  2:2408576/150‑1 (MQ=255)
                                                                                                         ccGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCa                                           >  1:2307499/1‑150 (MQ=255)
                                                                                                         ccGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCa                                           >  1:2930616/1‑150 (MQ=255)
                                                                                                          cGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCaa                                          >  2:2750261/1‑150 (MQ=255)
                                                                                                          cGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCaa                                          >  1:299926/1‑150 (MQ=255)
                                                                                                            ttCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGt                                        >  2:842555/1‑150 (MQ=255)
                                                                                                            ttCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGt                                        >  1:2467515/1‑150 (MQ=255)
                                                                                                              cTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGttt                                      <  2:489668/150‑1 (MQ=255)
                                                                                                              cTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGttt                                      >  1:3107223/1‑150 (MQ=255)
                                                                                                                 ggTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGAt                                   >  2:3709838/1‑150 (MQ=255)
                                                                                                                  gTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATg                                  <  2:1501668/150‑1 (MQ=255)
                                                                                                                    tAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGc                                <  2:703784/150‑1 (MQ=255)
                                                                                                                    tAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGc                                >  1:4512353/1‑150 (MQ=255)
                                                                                                                          aGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCatcgat                          <  1:3034677/150‑1 (MQ=255)
                                                                                                                          aGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCatcgat                          <  1:2165769/150‑1 (MQ=255)
                                                                                                                            cATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCatcgatcg                        >  1:4462954/1‑150 (MQ=255)
                                                                                                                            cATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCatcgatcg                        >  1:4462655/1‑150 (MQ=255)
                                                                                                                              tCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGa                      <  1:4319506/150‑1 (MQ=255)
                                                                                                                               cTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGaa                     <  1:848556/150‑1 (MQ=255)
                                                                                                                                  ccTTTGCACTCTCTACACTTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAgc                  <  1:736852/150‑1 (MQ=255)
                                                                                                                                    tttGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAgcgc                >  2:821066/1‑150 (MQ=255)
                                                                                                                                     ttGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAgcgcg               <  2:4469741/150‑1 (MQ=255)
                                                                                                                                      tGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAgcgcgc              >  2:4115073/1‑150 (MQ=255)
                                                                                                                                       gCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGcc             <  1:4197077/150‑1 (MQ=255)
                                                                                                                                        cACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCa            <  2:2357966/150‑1 (MQ=255)
                                                                                                                                         aCTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAg           >  1:3758384/1‑150 (MQ=255)
                                                                                                                                             tctACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGc       <  2:1241705/150‑1 (MQ=255)
                                                                                                                                                acacCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGaa    >  2:3530491/1‑150 (MQ=255)
                                                                                                                                                acacCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGaa    >  2:3794373/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTGGCAAACAAAAATCCTTGCCGCTCATTGATCCCGTTCTTGGTTAAAAAAGCATCTTCCTTTGCACTCTCTACACCTTCGGCGATCACCTGAAGATTCAATGCCTGGGCCACAGCGACGATCGCCCGGACCAGTGACTGCGAGACAGGTTGTTTGTGAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAG  >  CP009273/1339000‑1339295

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