Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,266,665 G→T 13.7% D423Y (GAT→TAT)  yjcC → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,266,6650GT13.7% 977.9 / 87.2 289D423Y (GAT→TAT) yjcCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (121/128);  new base T (24/15);  total (146/143)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 308 total aligned reads displayed.

TAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCT  >  CP009273/4266516‑4266813
                                                                                                                                                     |                                                                                                                                                    
tAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATtgatg                                                                                                                                                      <  2:3813589/150‑1 (MQ=255)
  gcgcGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTgattat                                                                                                                                                    >  1:98126/1‑150 (MQ=255)
     cgAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGAtttt                                                                                                                                                 >  2:3226527/1‑150 (MQ=255)
      gAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTg                                                                                                                                                >  2:40278/1‑150 (MQ=255)
      gAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTg                                                                                                                                                >  1:4008903/1‑150 (MQ=255)
      gAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTg                                                                                                                                                <  1:2681675/150‑1 (MQ=255)
        atcaatcaGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAATTGGCAATTGATGATTTTGGt                                                                                                                                              >  2:1173035/1‑150 (MQ=255)
           aatcaGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTAtt                                                                                                                                           <  1:2744738/150‑1 (MQ=255)
               aGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCt                                                                                                                                       >  1:4201613/1‑150 (MQ=255)
               aGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCt                                                                                                                                       <  2:1875407/150‑1 (MQ=255)
               aGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCt                                                                                                                                       >  1:2303076/1‑150 (MQ=255)
                 aaaaCAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTAc                                                                                                                                     >  2:232332/1‑150 (MQ=255)
                 aaaaCAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAAGTGATGATTTTGGTATTGGCTAc                                                                                                                                     >  1:2593448/1‑150 (MQ=255)
                    aCAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTActct                                                                                                                                  <  1:1116291/150‑1 (MQ=255)
                    aCAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTActct                                                                                                                                  <  1:3439849/150‑1 (MQ=255)
                          cAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTg                                                                                                                            <  2:4370755/150‑1 (MQ=255)
                          cAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCAATTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTACGCCAGGCAGGTTACGAAGAGGCAATTGATGATTTTGGTATTTGCTACTCTAACTTg                                                                                                                            >  1:1792197/1‑150 (MQ=255)
                           aaTACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGc                                                                                                                           <  1:145660/150‑1 (MQ=255)
                           aaTACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGc                                                                                                                           >  2:1810419/1‑150 (MQ=255)
                            aTACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCa                                                                                                                          >  1:2053433/1‑150 (MQ=255)
                                  ggTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCt                                                                                                                    <  1:2632358/150‑1 (MQ=255)
                                  ggTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCt                                                                                                                    <  2:3652358/150‑1 (MQ=255)
                                   gTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCtt                                                                                                                   <  2:3738632/150‑1 (MQ=255)
                                      cGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTaaa                                                                                                                <  1:732645/150‑1 (MQ=255)
                                       gTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAAt                                                                                                               >  1:1751777/1‑150 (MQ=255)
                                       gTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAAt                                                                                                               >  1:3372048/1‑150 (MQ=255)
                                       gTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAAt                                                                                                               >  1:535721/1‑150 (MQ=255)
                                        tCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATc                                                                                                              <  2:3438218/150‑1 (MQ=255)
                                          cGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCAt                                                                                                            >  2:4119839/1‑150 (MQ=255)
                                           gcagcaAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCAtt                                                                                                           <  1:1553626/150‑1 (MQ=255)
                                            cagcaAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTg                                                                                                          >  1:150822/1‑150 (MQ=255)
                                             agcaAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGa                                                                                                         >  2:28208/1‑150 (MQ=255)
                                             agcaAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGa                                                                                                         <  1:4618715/150‑1 (MQ=255)
                                             agcaAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGa                                                                                                         >  2:705506/1‑150 (MQ=255)
                                                aaattaaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATg                                                                                                      >  2:2432799/1‑150 (MQ=255)
                                                aaattaaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATg                                                                                                      >  1:3936409/1‑150 (MQ=255)
                                                 aattaaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGt                                                                                                     >  1:3198583/1‑150 (MQ=255)
                                                   ttaaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCg                                                                                                   <  2:3589980/150‑1 (MQ=255)
                                                      aattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGata                                                                                                >  2:1080968/1‑150 (MQ=255)
                                                          tGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATAtttt                                                                                            >  1:3297526/1‑150 (MQ=255)
                                                            aaGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGa                                                                                          >  1:3649739/1‑150 (MQ=255)
                                                             aGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGaa                                                                                         >  2:3342564/1‑150 (MQ=255)
                                                             aGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGaa                                                                                         >  1:388870/1‑150 (MQ=255)
                                                              gTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGaaa                                                                                        <  2:448943/150‑1 (MQ=255)
                                                                 aCTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTTAATCTTTGAATGTCGATATTTTGAAAATc                                                                                     <  2:1496124/150‑1 (MQ=255)
                                                                 aCTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATc                                                                                     >  2:3482360/1‑150 (MQ=255)
                                                                  cTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCg                                                                                    <  1:3986999/150‑1 (MQ=255)
                                                                  cTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCg                                                                                    <  2:2713515/150‑1 (MQ=255)
                                                                  cTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCg                                                                                    <  2:4238424/150‑1 (MQ=255)
                                                                    gAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGAc                                                                                  >  1:3752606/1‑150 (MQ=255)
                                                                    gAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGAc                                                                                  <  1:4348890/150‑1 (MQ=255)
                                                                       catgcatTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACaaa                                                                               >  2:540759/1‑150 (MQ=255)
                                                                           catTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGt                                                                           >  2:1304766/1‑150 (MQ=255)
                                                                           catTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGt                                                                           <  1:4381410/150‑1 (MQ=255)
                                                                            atTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGtt                                                                          <  1:4017112/150‑1 (MQ=255)
                                                                            atTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGtt                                                                          >  2:3144738/1‑150 (MQ=255)
                                                                            atTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGtt                                                                          <  2:3053431/150‑1 (MQ=255)
                                                                               tCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTttgt                                                                       <  2:646321/150‑1 (MQ=255)
                                                                               tCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTttgt                                                                       <  2:3872894/150‑1 (MQ=255)
                                                                                cTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTttgtt                                                                      >  2:1991453/1‑150 (MQ=255)
                                                                                  tGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGa                                                                    <  2:4441081/150‑1 (MQ=255)
                                                                                  tGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGa                                                                    <  1:2326403/150‑1 (MQ=255)
                                                                                   gATGTTGACACAATGACGCCGCTTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATGGGCGACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGaa                                                                   <  2:2757799/150‑1 (MQ=255)
                                                                                   gATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGaa                                                                   >  2:1599458/1‑150 (MQ=255)
                                                                                    aTGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGaaa                                                                  <  2:2123947/150‑1 (MQ=255)
                                                                                     tGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAAc                                                                 <  1:1226151/150‑1 (MQ=255)
                                                                                                 gACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGAccaccca                                                     <  2:35652/150‑1 (MQ=255)
                                                                                                  aCGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGAccacccac                                                    <  2:298529/150‑1 (MQ=255)
                                                                                                      cGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACaaaa                                                >  1:3077423/1‑150 (MQ=255)
                                                                                                        attattCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAAcc                                              >  1:458576/1‑150 (MQ=255)
                                                                                                             tCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCa                                         >  1:2442252/1‑150 (MQ=255)
                                                                                                             tCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCa                                         <  1:649595/150‑1 (MQ=255)
                                                                                                               tGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCAtt                                       >  1:4524751/1‑150 (MQ=255)
                                                                                                               tGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCAtt                                       >  1:3568144/1‑150 (MQ=255)
                                                                                                               tGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCAtt                                       >  2:77455/1‑150 (MQ=255)
                                                                                                               tGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCAtt                                       <  2:1180210/150‑1 (MQ=255)
                                                                                                                 gCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTg                                     <  1:1725310/150‑1 (MQ=255)
                                                                                                                 gCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTg                                     <  1:1581863/150‑1 (MQ=255)
                                                                                                                 gCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTg                                     <  2:1752327/150‑1 (MQ=255)
                                                                                                                  cTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGa                                    >  1:2221624/1‑150 (MQ=255)
                                                                                                                   tttCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGAt                                   >  1:4184931/1‑150 (MQ=255)
                                                                                                                    ttCAGCCATACAGGTTACGAAGTGGCAATACATAAATATGGTATTCGCTACTCTAACATGCATAACCTTAAATCAATGAAAGACGAAATTTTCAAAATCGACAAAGCCTATGATGAAACCCTGACCACCCACAAAACCAGTCaatagagt                                  >  2:2239606/1‑143 (MQ=255)
                                                                                                                      ccGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGc                                >  1:2140662/1‑150 (MQ=255)
                                                                                                                       cGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCg                               >  2:2534954/1‑150 (MQ=255)
                                                                                                                         ccAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGa                             >  2:1623017/1‑150 (MQ=255)
                                                                                                                            ggcaggTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAaca                          <  2:1731568/150‑1 (MQ=255)
                                                                                                                             gcaggTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAacac                         <  2:4201613/150‑1 (MQ=255)
                                                                                                                             gcaggTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAacac                         >  2:4289516/1‑150 (MQ=255)
                                                                                                                              caggTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAAcaca                        <  2:2391923/150‑1 (MQ=255)
                                                                                                                                   tACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACatcatc                   <  1:605924/150‑1 (MQ=255)
                                                                                                                                      gAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAg                <  1:4104509/150‑1 (MQ=255)
                                                                                                                                      gAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAg                <  1:2957725/150‑1 (MQ=255)
                                                                                                                                           ggCAATTGATGATTTTGGTATTGGCTACTCTATCTTGCATATCCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGATACGCTGACCACCCACTAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGgc           <  2:3290716/150‑1 (MQ=255)
                                                                                                                                           ggCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGgc           <  1:1852596/150‑1 (MQ=255)
                                                                                                                                           ggCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGgc           <  2:3311782/150‑1 (MQ=255)
                                                                                                                                            gCAATTTATGATTTTGGTAATGGCAACTCTAACTTTCATAAACTTAAATCATTGAATGTCGATATTTTGAAAATCGACCAATCGTATGATGAAACGCTGACCACCCACAAAACCAGTCATTAGATTGCGGAACACATCATCTATCTGgct          >  2:399705/1‑149 (MQ=255)
                                                                                                                                                ttGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGcaca      >  1:169516/1‑150 (MQ=255)
                                                                                                                                                 tgatgaTTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAg     >  1:2491306/1‑150 (MQ=255)
                                                                                                                                                   atgatTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGcc   >  2:666077/1‑150 (MQ=255)
                                                                                                                                                    tgatTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCt  >  2:40545/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCT  >  CP009273/4266516‑4266813

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