Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,249,776 G→A 6.8% S596S (AGC→AGT) ‡ PP_1983 ← PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,249,7760GA 6.8% 770.9 / 29.9 250S596S (AGC→AGT) ‡PP_1983PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (117/116);  new base A (9/8);  total (126/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 278 total aligned reads displayed.

TCGCGGGCGGTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCGGTGATGGTGATGGCGTCGGG  >  NC_002947/2249627‑2249925
                                                                                                                                                     |                                                                                                                                                     
tCGCGGGCGGTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAgcgcg                                                                                                                                                       >  1:523694/1‑150 (MQ=255)
tCGCGGGCGGTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAgcgcg                                                                                                                                                       >  2:8276385/1‑150 (MQ=255)
 cgcgGGCGGTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAgcgcgc                                                                                                                                                      >  2:7792365/1‑150 (MQ=255)
    gggcggTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCgcgctgc                                                                                                                                                   <  1:6575492/150‑1 (MQ=255)
      gcggTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCgcgctgcgc                                                                                                                                                 <  1:8619887/150‑1 (MQ=255)
          tCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCgctgct                                                                                                                                             <  2:4690294/150‑1 (MQ=255)
           cagcagcagGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCtt                                                                                                                                            <  1:6851054/150‑1 (MQ=255)
           cagcagcagGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCtt                                                                                                                                            <  2:878624/150‑1 (MQ=255)
           cagcagcagGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCtt                                                                                                                                            <  1:7885701/150‑1 (MQ=255)
              cagcagGCACTCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCCTCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTgtc                                                                                                                                         <  2:4048720/150‑1 (MQ=255)
              cagcagGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCAAGCAGCGCAGCGCgctgctggcc                                                                                                                                         <  2:1213439/150‑5 (MQ=255)
              cagcagGCAATCCACTTCGTTGAGGGTGATCGGTCCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGACAAGCAGCGCAGCGCgctgctggcc                                                                                                                                         <  2:891304/150‑5 (MQ=255)
                gcagGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTgtcga                                                                                                                                       <  2:5006117/150‑2 (MQ=255)
                  agGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGtcgtcg                                                                                                                                     <  2:6440396/150‑1 (MQ=255)
                       aTCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCa                                                                                                                                <  2:4335352/150‑1 (MQ=255)
                        tCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGACCAg                                                                                                                               >  1:806157/1‑150 (MQ=255)
                          cACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGAt                                                                                                                             <  2:7360735/150‑1 (MQ=255)
                          cACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGAt                                                                                                                             >  2:6189700/1‑150 (MQ=255)
                            cTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGAGCGTGTGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGc                                                                                                                           >  2:5422994/1‑150 (MQ=255)
                             ttCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCa                                                                                                                                <  2:5085684/144‑1 (MQ=255)
                             ttCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCa                                                                                                                                >  1:5085684/1‑144 (MQ=255)
                                 ttGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCCGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGAtc                                                                                                                      >  2:6622871/1‑150 (MQ=255)
                                 ttGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGAtc                                                                                                                      <  1:3874978/150‑1 (MQ=255)
                                               ttCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGtgtg                                                                                                        >  1:1245964/1‑150 (MQ=255)
                                                tCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTgg                                                                                                       <  2:9002497/150‑1 (MQ=255)
                                                   aCCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCg                                                                                                    >  1:6245706/1‑150 (MQ=255)
                                                   aCCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCg                                                                                                    >  1:2489666/1‑150 (MQ=255)
                                                   aCCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCg                                                                                                    >  2:8452954/1‑150 (MQ=255)
                                                    ccGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGa                                                                                                   >  1:5977683/1‑150 (MQ=255)
                                                     cGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGAt                                                                                                  >  1:1485130/1‑150 (MQ=255)
                                                       aaaCCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCa                                                                                                >  1:4625219/1‑150 (MQ=255)
                                                        aaCCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCAc                                                                                               <  1:9047649/150‑1 (MQ=255)
                                                         cccTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACt                                                                                              <  1:305608/149‑1 (MQ=255)
                                                          ccTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACtt                                                                                             <  2:1285331/150‑1 (MQ=255)
                                                           cTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTtc                                                                                            >  1:2816673/1‑150 (MQ=255)
                                                             cGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTtcat                                                                                          <  2:523694/150‑1 (MQ=255)
                                                             cGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTtcat                                                                                          <  1:8457192/150‑1 (MQ=255)
                                                               aCCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTtcatca                                                                                        <  2:2354144/150‑1 (MQ=255)
                                                               aCCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTtcatca                                                                                        <  2:6932087/150‑1 (MQ=255)
                                                                 cGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGt                                                                                      <  1:4853829/150‑1 (MQ=255)
                                                                   tGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGc                                                                                    <  1:7487386/150‑1 (MQ=255)
                                                                    gAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCt                                                                                   <  2:6599991/150‑1 (MQ=255)
                                                                      ggATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGt                                                                                 >  1:8298909/1‑150 (MQ=255)
                                                                       gATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGATCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTa                                                                                <  1:1023848/150‑1 (MQ=255)
                                                                       gATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGAGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTa                                                                                <  2:2385125/150‑1 (MQ=255)
                                                                              ccACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGc                                                                                    <  1:6273778/139‑1 (MQ=255)
                                                                               cACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCa                                                                        >  2:7240692/1‑150 (MQ=255)
                                                                               cACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCa                                                                        <  1:5727926/150‑1 (MQ=255)
                                                                                aCGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAg                                                                       >  1:7898258/1‑150 (MQ=255)
                                                                                 cGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAgg                                                                      <  1:4087867/150‑1 (MQ=255)
                                                                                    ttgttgCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAggcgg                                                                   >  1:7898050/1‑150 (MQ=255)
                                                                                    ttgttgCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAggcgg                                                                   <  2:8692581/150‑1 (MQ=255)
                                                                                     tgttgCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAggcggc                                                                  <  1:8710332/150‑1 (MQ=255)
                                                                                      gttgCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGCGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCa                                                                 <  2:1485228/150‑1 (MQ=255)
                                                                                         gcgcCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGaa                                                              <  2:6000721/150‑1 (MQ=255)
                                                                                         gcgcCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGaa                                                              >  1:8921676/1‑150 (MQ=255)
                                                                                             ccGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCcg                                                          >  2:1579548/1‑150 (MQ=255)
                                                                                             ccGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGAcg                                                          >  1:4266758/1‑150 (MQ=255)
                                                                                              cGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCcgt                                                         <  2:372621/150‑1 (MQ=255)
                                                                                               gAGCATCTCGCGGTGGTGGACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCcgtc                                                        <  1:3662175/150‑1 (MQ=255)
                                                                                                aGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCcgtcg                                                       <  2:1889066/150‑1 (MQ=255)
                                                                                                  cATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGtt                                                     <  1:7673431/150‑1 (MQ=255)
                                                                                                     ctcGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAc                                                  <  2:6848945/150‑1 (MQ=255)
                                                                                                     ctcGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAc                                                  <  1:465988/150‑1 (MQ=255)
                                                                                                      tcGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCGGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAcc                                                 <  2:6741495/150‑1 (MQ=255)
                                                                                                      tcGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAcc                                                 >  1:6956284/1‑150 (MQ=255)
                                                                                                      tcGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAcc                                                 >  2:3962938/1‑150 (MQ=255)
                                                                                                      gcgcgGTGGTGGACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAcc                                                 <  2:6290876/149‑1 (MQ=255)
                                                                                                        gcgGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTc                                               >  2:8051773/1‑150 (MQ=255)
                                                                                                        gcgGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTc                                               >  2:8921453/1‑150 (MQ=255)
                                                                                                              gtgtACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTa                                         >  2:3413238/1‑150 (MQ=255)
                                                                                                                gtACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGc                                       <  2:4391167/150‑1 (MQ=255)
                                                                                                                gtACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGc                                       >  1:1704649/1‑150 (MQ=255)
                                                                                                                  aCCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCgg                                     >  1:8224688/1‑150 (MQ=255)
                                                                                                                    ccGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGcc                                   >  1:2576027/1‑150 (MQ=255)
                                                                                                                    ccGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGcc                                   >  2:6143773/1‑150 (MQ=255)
                                                                                                                    ccGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGcc                                   >  1:6416158/1‑150 (MQ=255)
                                                                                                                         cgcgCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCtgt                              <  2:2816673/150‑1 (MQ=255)
                                                                                                                           cgcTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTGGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGGTGACCTCCAGGTAGCGGCCGCTGTgg                            >  1:4745060/1‑150 (MQ=255)
                                                                                                                            gcTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGc                           >  1:1813830/1‑150 (MQ=255)
                                                                                                                             cTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCg                          >  1:2057495/1‑150 (MQ=255)
                                                                                                                              tCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGGCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTCGCGGCCGCTGTGGCGt                         >  2:7688648/1‑150 (MQ=255)
                                                                                                                                gcgcAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACGTCCTCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTCGCGGCCGCTGTGGCGTTc                       >  2:7144821/1‑150 (MQ=255)
                                                                                                                                  gcAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCgg                     <  2:8374740/150‑1 (MQ=255)
                                                                                                                                    aaCTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtg                   >  1:2834221/1‑150 (MQ=255)
                                                                                                                                      cTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgat                 >  1:122549/1‑150 (MQ=255)
                                                                                                                                      cTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgat                 >  2:7664986/1‑150 (MQ=255)
                                                                                                                                      cTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgat                 <  1:7241217/150‑1 (MQ=255)
                                                                                                                                          gCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtg             >  1:5576929/1‑150 (MQ=255)
                                                                                                                                          acaagcagcGCAGCGCGCTGCTGGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtg             >  1:4146246/4‑150 (MQ=255)
                                                                                                                                           cgagttgcgcgcTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtga            <  2:3501132/144‑1 (MQ=255)
                                                                                                                                           ccAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtga            >  2:7006/1‑150 (MQ=255)
                                                                                                                                           ccAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtga            <  2:5538364/150‑1 (MQ=255)
                                                                                                                                            cagcagCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtgat           <  2:7862184/150‑1 (MQ=255)
                                                                                                                                            cagcagCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtgat           >  2:1104784/1‑150 (MQ=255)
                                                                                                                                            aagcagcGCAGCGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtgat           >  2:819165/2‑150 (MQ=255)
                                                                                                                                            aagcagcGCAGCGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtgat           >  2:7499077/2‑150 (MQ=255)
                                                                                                                                             agcagcGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGGTCGgtgatggtgatg          >  2:3044427/1‑150 (MQ=255)
                                                                                                                                              gcagcgcagcGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCGGtgatggtgatgg         >  2:2136090/1‑150 (MQ=255)
                                                                                                                                                    cgctgcgctgCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCGGTGATGGTGATGGCGTCgg   <  1:5347540/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCGGGCGGTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCGGTGATGGTGATGGCGTCGGG  >  NC_002947/2249627‑2249925

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: