Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,249,777 C→G 6.4% S596T (AGC→ACC) ‡ 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,7770CG 6.4% 835.5 / 34.3 250S596T (AGC→ACC) ‡PP_1983PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (117/117);  new base G (9/7);  total (126/124)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 275 total aligned reads displayed.

CGCGGGCGGTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCGGTGATGGTGATGGCGTCGGG  >  NC_002947/2249628‑2249925
                                                                                                                                                     |                                                                                                                                                    
cgcgGGCGGTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAgcgcgc                                                                                                                                                      >  2:7792365/1‑150 (MQ=255)
 gcgggcggTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCgcgct                                                                                                                                                     <  1:6787248/150‑1 (MQ=255)
   gggcggTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCgcgctgc                                                                                                                                                   <  1:9084172/150‑1 (MQ=255)
   gggcggTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCgcgctgc                                                                                                                                                   <  1:8047019/150‑1 (MQ=255)
     gcggTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCgcgctgcgc                                                                                                                                                 <  1:5713485/150‑1 (MQ=255)
     gcggTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCgcgctgcgc                                                                                                                                                 <  2:4444439/150‑1 (MQ=255)
       ggTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGcgctgcgctg                                                                                                                                               >  1:2813337/1‑150 (MQ=255)
         tCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCgctgct                                                                                                                                             <  2:4690294/150‑1 (MQ=255)
         tCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCgctgct                                                                                                                                             >  2:743092/1‑150 (MQ=255)
            gcagcagGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCAGCGCgctgctggc                                                                                                                                          <  1:6939729/150‑4 (MQ=255)
             cagcagGCACTCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCCTCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTgtc                                                                                                                                         <  2:4048720/150‑1 (MQ=255)
             cagcagGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTgtc                                                                                                                                         >  1:5286322/1‑150 (MQ=255)
             cagcagGCAATCCACTTCGTTGAGGGTGATCGGTCCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGACAAGCAGCGCAGCGCgctgctggcc                                                                                                                                         <  2:891304/150‑5 (MQ=255)
               gcagGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTgtcga                                                                                                                                       <  2:5006117/150‑2 (MQ=255)
                 agGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGtcgtcg                                                                                                                                     <  2:6440396/150‑1 (MQ=255)
                     aaTCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGccc                                                                                                                                 >  1:8374740/1‑150 (MQ=255)
                     aaTCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGccc                                                                                                                                 >  1:5774621/1‑150 (MQ=255)
                      aTCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCa                                                                                                                                <  2:4335352/150‑1 (MQ=255)
                         cACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGAt                                                                                                                             <  2:7360735/150‑1 (MQ=255)
                         cACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGAt                                                                                                                             >  2:6189700/1‑150 (MQ=255)
                           cTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGAGCGTGTGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGc                                                                                                                           >  2:5422994/1‑150 (MQ=255)
                           cTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGc                                                                                                                           <  2:1802016/150‑1 (MQ=255)
                            ttCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCa                                                                                                                                <  2:5085684/144‑1 (MQ=255)
                                ttGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGAtc                                                                                                                      <  1:8479834/150‑1 (MQ=255)
                                       tGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATaca                                                                                                               >  1:7546903/1‑150 (MQ=255)
                                           cGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCgt                                                                                                           <  2:7897520/150‑1 (MQ=255)
                                             gTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCgtgt                                                                                                         >  2:290658/1‑150 (MQ=255)
                                              ttCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGtgtg                                                                                                        >  1:1245964/1‑150 (MQ=255)
                                              ttCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGtgtg                                                                                                        >  1:8729827/1‑150 (MQ=255)
                                               tCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTgg                                                                                                       <  2:9002497/150‑1 (MQ=255)
                                                  aCCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCg                                                                                                    >  2:8452954/1‑150 (MQ=255)
                                                  aCCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCg                                                                                                    >  1:2489666/1‑150 (MQ=255)
                                                    cGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGAt                                                                                                  >  1:1485130/1‑150 (MQ=255)
                                                     gAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATc                                                                                                 >  1:7888326/1‑150 (MQ=255)
                                                      aaaCCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCa                                                                                                <  2:6965252/150‑1 (MQ=255)
                                                        cccTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACt                                                                                              <  1:305608/149‑1 (MQ=255)
                                                         ccTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACtt                                                                                             <  2:1285331/150‑1 (MQ=255)
                                                           tCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTtca                                                                                           >  2:1112675/1‑150 (MQ=255)
                                                            cGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTtcat                                                                                          <  2:523694/150‑1 (MQ=255)
                                                             gACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTtcatc                                                                                         >  1:7731989/1‑150 (MQ=255)
                                                              aCCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTtcatca                                                                                        <  2:6932087/150‑1 (MQ=255)
                                                              aCCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTtcatca                                                                                        <  2:2354144/150‑1 (MQ=255)
                                                              aCCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTtcatca                                                                                        <  2:660499/150‑1 (MQ=255)
                                                                cGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGt                                                                                      <  1:4853829/150‑1 (MQ=255)
                                                                   gAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCt                                                                                   <  2:6599991/150‑1 (MQ=255)
                                                                      gATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGATCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTa                                                                                <  1:1023848/150‑1 (MQ=255)
                                                                        tatCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGc                                                                              <  1:8276385/150‑1 (MQ=255)
                                                                             ccACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGc                                                                                    <  1:6273720/139‑1 (MQ=255)
                                                                              cACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCtt                                                                                                                                            <  2:1944543/82‑1 (MQ=255)
                                                                              cACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCa                                                                        <  1:7319218/150‑1 (MQ=255)
                                                                              cACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCa                                                                        >  2:7240692/1‑150 (MQ=255)
                                                                               aCGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAg                                                                       >  1:6387829/1‑150 (MQ=255)
                                                                               aCGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGCACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAg                                                                       >  1:7898234/1‑150 (MQ=255)
                                                                                cGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAgg                                                                      <  1:4087867/150‑1 (MQ=255)
                                                                                 gCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAggc                                                                     <  1:8667565/150‑1 (MQ=255)
                                                                                   ttgttgCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAggcgg                                                                   <  2:8692581/150‑1 (MQ=255)
                                                                                     gttgCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCa                                                                 <  1:4499883/150‑1 (MQ=255)
                                                                                     gttgCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGCGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCa                                                                 <  2:1485228/150‑1 (MQ=255)
                                                                                        gcgcCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGaa                                                              >  1:8921676/1‑150 (MQ=255)
                                                                                        gcgcCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGaa                                                              <  2:6000721/150‑1 (MQ=255)
                                                                                            ccGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGAcg                                                          >  1:4266758/1‑150 (MQ=255)
                                                                                              gAGCATCTCGCGGTGGTGGACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCcgtc                                                        <  1:3662175/150‑1 (MQ=255)
                                                                                                   tctcGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGa                                                   >  1:5984303/1‑150 (MQ=255)
                                                                                                   tcacgcgGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGa                                                   >  1:3744161/4‑150 (MQ=255)
                                                                                                    ctcGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAc                                                  <  1:4641336/150‑1 (MQ=255)
                                                                                                    ctcGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAc                                                  <  2:6848945/150‑1 (MQ=255)
                                                                                                     tcGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCGGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAcc                                                 <  2:6741495/150‑1 (MQ=255)
                                                                                                     tcGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAcc                                                 >  2:3786390/1‑150 (MQ=255)
                                                                                                     tcGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAcc                                                 >  2:3962938/1‑150 (MQ=255)
                                                                                                     gcgcgGTGGTGGACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGAcc                                                 <  2:6290876/149‑1 (MQ=255)
                                                                                                      cgcgGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCt                                                >  1:6421165/1‑150 (MQ=255)
                                                                                                       gcgGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTc                                               >  2:8921453/1‑150 (MQ=255)
                                                                                                             gtgtACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTa                                         >  2:3413238/1‑150 (MQ=255)
                                                                                                               gtACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGc                                       >  1:1704649/1‑150 (MQ=255)
                                                                                                               gtACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGc                                       <  1:8189599/150‑1 (MQ=255)
                                                                                                                 tgctggCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCgg                                     >  1:7095931/5‑150 (MQ=255)
                                                                                                                   ccGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGcc                                   >  1:2576027/1‑150 (MQ=255)
                                                                                                                   ccGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGcc                                   >  2:6143773/1‑150 (MQ=255)
                                                                                                                   ccGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGcc                                   >  1:4954288/1‑150 (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)
                                                                                                                               gcgcAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCCCTGCCTCAGTGCGGTAGCCGGTCAGGCGGCAGAAGCCGGCGTTGACCGCCAGGTAGCGGCCGCTGTGGCgggc                       >  2:8865840/1‑147 (MQ=255)
                                                                                                                               gcgcAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACGTCCTCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTCGCGGCCGCTGTGGCGTTc                       >  2:7144821/1‑150 (MQ=255)
                                                                                                                                 gcAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCgg                     <  1:5517283/150‑1 (MQ=255)
                                                                                                                                 gcAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCgg                     <  2:8374740/150‑1 (MQ=255)
                                                                                                                                   aaCTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtg                   >  1:2834221/1‑150 (MQ=255)
                                                                                                                                    aCTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtga                  <  2:2482060/150‑1 (MQ=255)
                                                                                                                                     cTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgat                 >  1:122549/1‑150 (MQ=255)
                                                                                                                                     cTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgat                 >  2:7653593/1‑150 (MQ=255)
                                                                                                                                         acaagcagcGCAGCGCGCTGCTGGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtg             >  1:4146246/4‑150 (MQ=255)
                                                                                                                                          cgagttgcgcgcTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtga            <  2:3501132/144‑1 (MQ=255)
                                                                                                                                          ccAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtga            <  2:5538364/150‑1 (MQ=255)
                                                                                                                                          ccAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtga            <  2:1220636/150‑1 (MQ=255)
                                                                                                                                           cagcagCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtgat           <  2:7862142/150‑1 (MQ=255)
                                                                                                                                           aagcagcGCAGCGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtgat           >  2:819165/2‑150 (MQ=255)
                                                                                                                                           aagcagcGCAGCGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCggtgatggtgat           >  2:8332392/2‑150 (MQ=255)
                                                                                                                                            agcagcGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGGTCGgtgatggtgatg          >  2:3044427/1‑150 (MQ=255)
                                                                                                                                             gcagcgcagcGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCGGtgatggtgatgg         >  2:2136090/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCGGGCGGTCAGCAGCAGGCAATCCACTTCGTTGAGGGTGATCGGTTCCACCGAAACCTCGACCGTGAGGATATCGCCACGCTTGTTGCGCCCGAGCATCTCGCGGTGGTGTACCCGGCCGCGCTCGCGCAACTCGGCCAGCAGCGCGCTGCGCTGCTTGTCGTCGGCCCAGATGCCGATCTCATACACCGTGTGGCCGATCACTTCATCAGTGCTGTAGCCGGTCAGGCGGCAGAAGCCGTCGTTGACCTCCAGGTAGCGGCCGCTGTGGCGTTCGGTGATGGTGATGGCGTCGGG  >  NC_002947/2249628‑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: