Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,266,665 G→T 19.0% 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,6650GT19.0% 844.9 / 134.0 289D423Y (GAT→TAT) yjcCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (107/127);  new base T (31/24);  total (138/151)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 300 total aligned reads displayed.

TAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCT  >  CP009273/4266516‑4266813
                                                                                                                                                     |                                                                                                                                                    
     cgAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGAtttt                                                                                                                                                 >  1:3645678/1‑150 (MQ=255)
      gAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTg                                                                                                                                                <  1:4497823/150‑1 (MQ=255)
          caatcaGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTAt                                                                                                                                            <  1:76121/150‑1 (MQ=255)
              cttaataCAGAGCTATTGGCGGTGCGTCCGCAGCAATTTATATTTGTAGTGTCTGTGCTTGCATTTCTTGATGTTGTCATTATGTCGCCGATTATTGTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCATTTGATGATTTTGGTATTggg                                                                                                                                        <  2:1437836/144‑2 (MQ=255)
              caGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGc                                                                                                                                        >  1:2490609/1‑150 (MQ=255)
               aGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCt                                                                                                                                       >  2:582689/1‑150 (MQ=255)
                gAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTa                                                                                                                                      <  2:1836919/150‑1 (MQ=255)
                  aaaCAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTAct                                                                                                                                    <  1:4376796/150‑1 (MQ=255)
                  aaaCAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTAct                                                                                                                                    <  2:2968244/150‑1 (MQ=255)
                   aaCAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTActc                                                                                                                                   <  1:3906482/150‑1 (MQ=255)
                    aCAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTActct                                                                                                                                  <  1:3822967/150‑1 (MQ=255)
                         gCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACtt                                                                                                                             <  2:419388/150‑1 (MQ=255)
                         gCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACtt                                                                                                                             <  2:4023394/150‑1 (MQ=255)
                         gCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACtt                                                                                                                             <  2:1236259/150‑1 (MQ=255)
                          cAATCCGCGGTGCGTCCGCAGAAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGTCAGGTTACGAAGTGGCAATTTATTATTTTGGTATTGGCTACTCTAACTTg                                                                                                                            <  1:1024644/150‑1 (MQ=255)
                                gcgGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAAc                                                                                                                      <  1:1067160/150‑1 (MQ=255)
                                 cgGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAAcc                                                                                                                     >  2:689233/1‑150 (MQ=255)
                                  ggTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCt                                                                                                                    >  1:999686/1‑150 (MQ=255)
                                      cGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTaaa                                                                                                                >  1:2533071/1‑150 (MQ=255)
                                      cGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTaaa                                                                                                                <  2:117103/150‑1 (MQ=255)
                                       gTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAAt                                                                                                               >  2:3438981/1‑150 (MQ=255)
                                        tCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATc                                                                                                              >  2:4462317/1‑150 (MQ=255)
                                        tCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATc                                                                                                              <  1:3312010/150‑1 (MQ=255)
                                         ccGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCa                                                                                                             <  1:3365055/150‑1 (MQ=255)
                                          cGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCAt                                                                                                            >  1:467550/1‑150 (MQ=255)
                                          cGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCAt                                                                                                            <  2:2693442/150‑1 (MQ=255)
                                           gcagcaAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCAtt                                                                                                           <  2:736266/150‑1 (MQ=255)
                                            cagcaAATTAAATTTGAAGTGACTGAGCATGCATCTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTg                                                                                                          <  1:1946568/150‑1 (MQ=255)
                                               caAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAAt                                                                                                       >  1:1248523/1‑150 (MQ=255)
                                                aaattaaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATg                                                                                                      >  2:2364310/1‑150 (MQ=255)
                                                aaattaaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATg                                                                                                      >  1:357922/1‑150 (MQ=255)
                                                 aattaaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTCAATCATTGAATGt                                                                                                     >  2:4371828/1‑150 (MQ=255)
                                                 aattaaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGt                                                                                                     <  1:428944/150‑1 (MQ=255)
                                                  attaaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTc                                                                                                    >  1:836736/1‑150 (MQ=255)
                                                   ttaaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCg                                                                                                   <  1:2858248/150‑1 (MQ=255)
                                                   ttaaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCg                                                                                                   <  1:546422/150‑1 (MQ=255)
                                                    taaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGa                                                                                                  >  1:2120520/1‑150 (MQ=255)
                                                    taaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGa                                                                                                  >  1:666289/1‑150 (MQ=255)
                                                     aaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGat                                                                                                 >  1:2920418/1‑150 (MQ=255)
                                                     aaattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGat                                                                                                 >  1:2192714/1‑150 (MQ=255)
                                                      aattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTAAGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGata                                                                                                >  1:896149/1‑150 (MQ=255)
                                                       attTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGatat                                                                                               <  1:2354662/150‑1 (MQ=255)
                                                        tttGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATAtt                                                                                              <  1:4483169/150‑1 (MQ=255)
                                                        tttGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATAtt                                                                                              <  1:2855294/150‑1 (MQ=255)
                                                        tttGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATAtt                                                                                              <  2:4258270/150‑1 (MQ=255)
                                                        tttGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATAtt                                                                                              <  2:4419407/150‑1 (MQ=255)
                                                        tttGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATAtt                                                                                              <  2:4073/150‑1 (MQ=255)
                                                          tGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATAtttt                                                                                            <  2:1795406/150‑1 (MQ=255)
                                                           gAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTg                                                                                           >  2:2501748/1‑150 (MQ=255)
                                                            aaGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGa                                                                                          >  1:4220018/1‑150 (MQ=255)
                                                            aaGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGa                                                                                          >  2:2451247/1‑150 (MQ=255)
                                                             aGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGaa                                                                                         >  1:4151538/1‑150 (MQ=255)
                                                                     aGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACa                                                                                 <  1:3046307/150‑1 (MQ=255)
                                                                     aGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACa                                                                                 <  1:2832933/150‑1 (MQ=255)
                                                                        atgcatTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAAt                                                                              >  2:4508590/1‑150 (MQ=255)
                                                                         tccatTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATc                                                                             <  2:2952303/148‑1 (MQ=255)
                                                                          gcatTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCg                                                                            >  1:2555423/1‑150 (MQ=255)
                                                                          gcatTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCg                                                                            <  1:2719594/150‑1 (MQ=255)
                                                                           catTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGt                                                                           >  2:1920224/1‑150 (MQ=255)
                                                                           catTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGt                                                                           >  2:2998515/1‑150 (MQ=255)
                                                                           catTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGt                                                                           <  2:3035807/150‑1 (MQ=255)
                                                                           catTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGt                                                                           >  2:3839392/1‑150 (MQ=255)
                                                                            atTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATATTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGACGATATTTTGAAAATCGACAAATCgat                                                                          >  1:2435768/1‑148 (MQ=255)
                                                                             tttCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGttt                                                                         <  1:1132259/150‑1 (MQ=255)
                                                                                 ttGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTttgttg                                                                     <  1:942414/150‑1 (MQ=255)
                                                                                   gATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGaa                                                                   >  2:2594663/1‑150 (MQ=255)
                                                                                    aTGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGaaa                                                                  <  2:2938892/150‑1 (MQ=255)
                                                                                     tGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAAc                                                                 <  2:3758851/150‑1 (MQ=255)
                                                                                        tGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCt                                                              <  2:3249044/150‑1 (MQ=255)
                                                                                                tGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGAccaccc                                                      <  1:2004792/150‑1 (MQ=255)
                                                                                                   cGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCcaca                                                   <  1:2996648/150‑1 (MQ=255)
                                                                                                   cGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCcaca                                                   <  2:3960777/150‑1 (MQ=255)
                                                                                                    gCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACaa                                                  >  2:1653295/1‑150 (MQ=255)
                                                                                                    gCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACaa                                                  >  2:4346663/1‑150 (MQ=255)
                                                                                                           attCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGt                                           >  2:2603709/1‑150 (MQ=255)
                                                                                                            ttCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTc                                          <  2:4333115/150‑1 (MQ=255)
                                                                                                              cTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCAt                                        >  2:158333/1‑150 (MQ=255)
                                                                                                              cTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCAt                                        <  1:4533239/150‑1 (MQ=255)
                                                                                                              cTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCAt                                        >  1:149684/1‑150 (MQ=255)
                                                                                                               tGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCAtt                                       >  2:133452/1‑150 (MQ=255)
                                                                                                                ggCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCAttt                                      <  2:2747674/150‑1 (MQ=255)
                                                                                                                ggCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCAgtt                                      >  1:157651/1‑150 (MQ=255)
                                                                                                                   tttCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGAt                                   >  1:3977327/1‑150 (MQ=255)
                                                                                                                   tttCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGAt                                   <  2:557661/150‑1 (MQ=255)
                                                                                                                    ttCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGAtt                                  <  2:475449/150‑1 (MQ=255)
                                                                                                                       cGCCAGGCAGGTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCg                               >  1:4168621/1‑150 (MQ=255)
                                                                                                                       cGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCg                               >  2:2064325/1‑150 (MQ=255)
                                                                                                                          caggcaggTTACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGaa                            >  2:4129763/1‑150 (MQ=255)
                                                                                                                            ggcaggTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAaca                          <  1:4445307/150‑1 (MQ=255)
                                                                                                                             gcaggTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAacac                         <  2:3216074/150‑1 (MQ=255)
                                                                                                                                ggTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACAcatc                      >  1:959371/1‑150 (MQ=255)
                                                                                                                                  ttACGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACAcatcat                    >  1:139418/1‑150 (MQ=255)
                                                                                                                                  ttACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACAcatcat                    <  2:1341997/150‑1 (MQ=255)
                                                                                                                                     cGAAGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGa                 >  1:3200622/1‑150 (MQ=255)
                                                                                                                                      gAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAg                >  1:1625498/1‑150 (MQ=255)
                                                                                                                                       aaGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGc               >  1:1507095/1‑150 (MQ=255)
                                                                                                                                        aGTGGCAATTGATTATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCt              >  2:4309078/1‑150 (MQ=255)
                                                                                                                                             cAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGgcgc         >  2:1511391/1‑150 (MQ=255)
                                                                                                                                                  gatgatTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGc    >  2:3771188/1‑150 (MQ=255)
                                                                                                                                                    tgatTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCt  <  1:260915/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCT  >  CP009273/4266516‑4266813

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: