Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 718,686 C→T 5.6% W395* (TGG→TGA kdpD ← fused sensory histidine kinase in two‑component regulatory system with KdpE: signal sensing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273718,6860CT 5.6% 1076.6 / 34.1 301W395* (TGG→TGAkdpDfused sensory histidine kinase in two‑component regulatory system with KdpE: signal sensing protein
Reads supporting (aligned to +/- strand):  ref base C (142/141);  new base T (6/11);  total (149/152)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 307 total aligned reads displayed.

CGTCCATAAAATAGCGCCACCACCACCACGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGCCACCAGCGGCGCGAGGC  >  CP009273/718537‑718835
                                                                                                                                                     |                                                                                                                                                     
  tCCATAAAATAGCGCCACCACCACCACGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCa                                                                                                                                                     <  2:5203201/150‑1 (MQ=255)
      tAAAATAGCGCCACCACCACCACGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGTCACTTg                                                                                                                                                 <  2:2395413/150‑1 (MQ=255)
      tAAAATAGCGCCACCACCACCACGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTg                                                                                                                                                 <  1:3243675/150‑1 (MQ=255)
       aaaaTAGCGCCACCACCACCACGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGt                                                                                                                                                >  1:3076995/1‑150 (MQ=255)
         aaTAGCGCCACCACCACCACGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCt                                                                                                                                              <  1:4657009/150‑1 (MQ=255)
           tAGCGCCACCACCACCACGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCttt                                                                                                                                            >  1:2066240/1‑150 (MQ=255)
             gcgcCACCACCACCACGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTaa                                                                                                                                          <  2:1445575/150‑1 (MQ=255)
              cgcCACCACCACCACGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTaaa                                                                                                                                         >  1:50316/1‑150 (MQ=255)
              cgcCACCACCACCACGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAAGTTGTACACGCCACTTGTCTTTaaa                                                                                                                                         >  2:645224/1‑150 (MQ=255)
                    caccaccacGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCg                                                                                                                                   >  2:1613502/1‑150 (MQ=255)
                     accaccacGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACCGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCgg                                                                                                                                  <  2:4656786/150‑1 (MQ=255)
                          cacGCCAAGCAGATACAGCATCACCTGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATTACGTCGCATAACGCGGCGGCATCCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATc                                                                                                                             <  1:5150940/150‑1 (MQ=255)
                          cacGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATc                                                                                                                             >  2:2527812/1‑150 (MQ=255)
                          cacGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATc                                                                                                                             <  1:2696551/150‑1 (MQ=255)
                            cGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCg                                                                                                                           <  1:225144/150‑1 (MQ=255)
                             gCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCgg                                                                                                                          >  1:3755573/1‑150 (MQ=255)
                              ccAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGgc                                                                                                                         <  2:1005811/150‑1 (MQ=255)
                              ccAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGgc                                                                                                                         >  1:4886825/1‑150 (MQ=255)
                              ccAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGgc                                                                                                                         >  1:441686/1‑150 (MQ=255)
                                 aGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGTCACTTGTCTTTAAAAGAGCGGTTATCCGGcgcg                                                                                                                      >  1:1973565/1‑150 (MQ=255)
                                   cAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCgtt                                                                                                                    <  1:5533162/150‑1 (MQ=255)
                                    aGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCgttg                                                                                                                   >  2:2088772/1‑150 (MQ=255)
                                      aTACAGCATCACCAGGTTGGCGGCATCAATCGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCTCATAACGCGGCGGCAACCACGCATCCCTGATTTTGTACACGCAACTTGTCTTTAAAAGAGCGGTTATCCGGCGCgttgtt                                                                                                                 <  2:1519102/150‑1 (MQ=255)
                                       tACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTa                                                                                                                >  2:4512161/1‑150 (MQ=255)
                                          aGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATc                                                                                                             <  1:2952496/150‑1 (MQ=255)
                                          aGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATc                                                                                                             >  2:3783087/1‑150 (MQ=255)
                                             aTCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTg                                                                                                          <  2:573390/150‑1 (MQ=255)
                                             aTCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTg                                                                                                          >  1:469997/1‑150 (MQ=255)
                                               cACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTgcg                                                                                                        <  2:81973/150‑1 (MQ=255)
                                                aCCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTgcgg                                                                                                       <  1:1028112/150‑1 (MQ=255)
                                                 ccAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTgcggg                                                                                                      >  2:1979950/1‑150 (MQ=255)
                                                   aGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGTCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTgcgggcg                                                                                                    <  1:2819508/150‑1 (MQ=255)
                                                    ggTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTGAATCGTgcgggcgg                                                                                                   >  2:3336390/1‑150 (MQ=255)
                                                      ttGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGggt                                                                                                 <  2:2067677/150‑1 (MQ=255)
                                                        ggcggcATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGggtgg                                                                                               >  1:4708559/1‑150 (MQ=255)
                                                           ggcATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTtc                                                                                            >  2:747821/1‑150 (MQ=255)
                                                             cATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTtcat                                                                                          <  1:3487840/150‑1 (MQ=255)
                                                             cATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTtcat                                                                                          <  1:636811/150‑1 (MQ=255)
                                                              aTCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTGCACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTtcatc                                                                                         >  1:4094520/1‑150 (MQ=255)
                                                               tCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTtcatca                                                                                        >  1:4529202/1‑150 (MQ=255)
                                                               tCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTtcatca                                                                                        >  2:3570446/1‑150 (MQ=255)
                                                                cAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTCAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCaa                                                                                       <  2:439350/150‑1 (MQ=255)
                                                                        aTCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGAcc                                                                               >  2:1403860/1‑150 (MQ=255)
                                                                        aTCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGAcc                                                                               <  1:2336039/150‑1 (MQ=255)
                                                                           aGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGTCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGc                                                                            <  1:530657/150‑1 (MQ=255)
                                                                             ccACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCAc                                                                          >  1:2426483/1‑150 (MQ=255)
                                                                             ccACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCAc                                                                          >  1:5572934/1‑150 (MQ=255)
                                                                             ccACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCAc                                                                          >  1:1682331/1‑150 (MQ=255)
                                                                              cACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCAcc                                                                         >  1:1811924/1‑150 (MQ=255)
                                                                                cTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTg                                                                       >  2:2008786/1‑150 (MQ=255)
                                                                                cTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTg                                                                       >  2:4054388/1‑150 (MQ=255)
                                                                                cTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTg                                                                       <  2:4159495/150‑1 (MQ=255)
                                                                                 tGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGa                                                                      >  1:2191462/1‑150 (MQ=255)
                                                                                 tGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGa                                                                      >  1:2306828/1‑150 (MQ=255)
                                                                                  gCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGTCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGAt                                                                     >  2:2328132/1‑150 (MQ=255)
                                                                                    aTGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCg                                                                   <  2:1235259/150‑1 (MQ=255)
                                                                                    aTGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCg                                                                   >  2:5099536/1‑150 (MQ=255)
                                                                                    aTGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCg                                                                   >  1:1578736/1‑150 (MQ=255)
                                                                                    aTGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCg                                                                   <  2:3878267/150‑1 (MQ=255)
                                                                                       gCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCgaga                                                                >  1:501050/1‑150 (MQ=255)
                                                                                         aaTTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGTCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATc                                                              <  2:3105647/150‑1 (MQ=255)
                                                                                         aaTTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATc                                                              <  1:4747804/150‑1 (MQ=255)
                                                                                         aaTTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATc                                                              <  1:423846/150‑1 (MQ=255)
                                                                                         aaTTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATc                                                              <  1:1915207/150‑1 (MQ=255)
                                                                                         aaTTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATc                                                              <  1:217001/150‑1 (MQ=255)
                                                                                             aaGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCgggg                                                          <  2:3039762/150‑1 (MQ=255)
                                                                                              aGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCggggg                                                         <  1:1793219/150‑1 (MQ=255)
                                                                                                   aTAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATg                                                    >  1:1383345/1‑150 (MQ=255)
                                                                                                    tAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATgc                                                   >  1:1220955/1‑150 (MQ=255)
                                                                                                    tAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATgc                                                   >  1:4968675/1‑150 (MQ=255)
                                                                                                     aaCGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATgcg                                                  >  1:1463533/1‑150 (MQ=255)
                                                                                                         gcgcATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGcc                                              <  2:2489637/150‑1 (MQ=255)
                                                                                                                   cggcggCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAg                                    <  2:5494762/150‑1 (MQ=255)
                                                                                                                    ggcggcAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGc                                   >  1:5602206/1‑150 (MQ=255)
                                                                                                                     gcggcAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCa                                  <  2:3424938/150‑1 (MQ=255)
                                                                                                                        gcAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAAc                               <  2:4813476/150‑1 (MQ=255)
                                                                                                                            ccACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGttt                           >  1:4015707/1‑150 (MQ=255)
                                                                                                                              aCGCATCACTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTcc                         <  2:4792656/150‑1 (MQ=255)
                                                                                                                                gCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCg                       <  2:1113107/150‑1 (MQ=255)
                                                                                                                                  aTCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCACACGTTTCCCGAc                     >  2:3201439/1‑150 (MQ=255)
                                                                                                                                   tCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACg                    >  2:4939872/1‑150 (MQ=255)
                                                                                                                                   tCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACg                    >  1:5488342/1‑150 (MQ=255)
                                                                                                                                   tCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACg                    >  2:534878/1‑150 (MQ=255)
                                                                                                                                   tCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACg                    <  2:1704809/150‑1 (MQ=255)
                                                                                                                                   tCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACg                    >  2:841565/1‑150 (MQ=255)
                                                                                                                                     ccTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGcc                  <  2:614858/150‑1 (MQ=255)
                                                                                                                                     ccTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGcc                  <  2:3467467/150‑1 (MQ=255)
                                                                                                                                       tGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGccac                >  1:3375423/1‑150 (MQ=255)
                                                                                                                                       tGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGccac                >  1:520338/1‑150 (MQ=255)
                                                                                                                                       tGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGccac                >  2:3857181/1‑150 (MQ=255)
                                                                                                                                         aagctcgacacGCCACTTGTCTTCAAAAGAGCGGTGATCCGGCGCGTGGTTAATCGTGCGGGCGGGTGGTTCATCACGCGCGACCAGCACCTGATCGAGATCGGGGGCGAGGCGCGCCAGGCGGTCAGCAAACGTTTCCCGACGccacca              <  2:2936357/143‑1 (MQ=255)
                                                                                                                                          aTTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGCCACCAg             >  1:888065/1‑150 (MQ=255)
                                                                                                                                          aTTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGCCACCAg             <  2:3632462/150‑1 (MQ=255)
                                                                                                                                            ttGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGCCACCAgcg           >  2:3025565/1‑150 (MQ=255)
                                                                                                                                            ttGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGCCACCAgcg           <  2:1943282/150‑1 (MQ=255)
                                                                                                                                                acacGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGCCACCAGCGgcgc       <  1:3008207/150‑1 (MQ=255)
                                                                                                                                                acacGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGCCACCAGCGgcgc       <  2:263268/150‑1 (MQ=255)
                                                                                                                                                 cacGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGCCACCAGCGGcgcg      >  1:2663851/1‑150 (MQ=255)
                                                                                                                                                     ccACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGCCACCAGCGGCGCGAGGc  <  1:365950/150‑1 (MQ=255)
                                                                                                                                                     ccACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGCCACCAGCGGCGCGAGGc  >  2:5316493/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTCCATAAAATAGCGCCACCACCACCACGCCAAGCAGATACAGCATCACCAGGTTGGCGGCATCAAACGCCATCAGCCACTGCATGGCAATTAAGGTGATAACGGCGCATAACGCGGCGGCAACCACGCATCCCTGAATTTGTACACGCCACTTGTCTTTAAAAGAGCGGTTATCCGGCGCGTTGTTAATCGTGCGGGCGGGTGGTTCATCAAGCGCGACCAGCACCTGATCGAGATCGGGGGCGATGCGCGCCAGTCGGTCAGCAAACGTTTCCCGACGCCACCAGCGGCGCGAGGC  >  CP009273/718537‑718835

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: