Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 104,582 A→T 17.9% T131S (ACG→TCG)  secM → regulator of secA translation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273104,5820AT17.9% 735.5 / 152.6 302T131S (ACG→TCG) secMregulator of secA translation
Reads supporting (aligned to +/- strand):  ref base A (143/105);  new base T (23/31);  total (166/136)
Fisher's exact test for biased strand distribution p-value = 5.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

Only 100 of 308 total aligned reads displayed.

GGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTACTTCATTTT  >  CP009273/104433‑104729
                                                                                                                                                     |                                                                                                                                                   
ggCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCa                                                                                                                                                     <  2:2350485/150‑1 (MQ=255)
ggCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCa                                                                                                                                                     <  2:330108/150‑1 (MQ=255)
 gCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCAc                                                                                                                                                    <  1:3972743/150‑1 (MQ=255)
   aTCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGc                                                                                                                                                  <  1:517130/150‑1 (MQ=255)
    tCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGcc                                                                                                                                                 >  2:2782904/1‑150 (MQ=255)
     cAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCg                                                                                                                                                >  1:1580445/1‑150 (MQ=255)
      aaCATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGt                                                                                                                                               <  1:4002742/150‑1 (MQ=255)
       aCATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTc                                                                                                                                              <  2:443225/150‑1 (MQ=255)
        cATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCt                                                                                                                                             >  2:638090/1‑150 (MQ=255)
           gCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGaa                                                                                                                                          >  1:1142283/1‑150 (MQ=255)
             cATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGaaaa                                                                                                                                        >  2:1990174/1‑150 (MQ=255)
             cATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGaaaa                                                                                                                                        >  2:4015034/1‑150 (MQ=255)
              aTTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAg                                                                                                                                       <  1:1664103/150‑1 (MQ=255)
                    aCGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTAt                                                                                                                                 <  1:3498432/150‑1 (MQ=255)
                      ggTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCg                                                                                                                               >  2:54218/1‑150 (MQ=255)
                       gTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAATGCTCGCCGTCTGAAAAGGGTTATCGc                                                                                                                              >  2:1219533/1‑150 (MQ=255)
                       gTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGc                                                                                                                              >  2:1219028/1‑150 (MQ=255)
                         aaTCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCAt                                                                                                                            <  2:2009909/150‑1 (MQ=255)
                           tCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTg                                                                                                                          >  1:2016933/1‑150 (MQ=255)
                              gTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGAtt                                                                                                                       >  1:3394075/1‑150 (MQ=255)
                              gTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGAtt                                                                                                                       >  2:2337878/1‑150 (MQ=255)
                                   ctttctttCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATgcg                                                                                                                  >  2:1814639/1‑150 (MQ=255)
                                    tttctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATgcgc                                                                                                                 <  2:3135893/150‑1 (MQ=255)
                                     ttctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCa                                                                                                                >  2:846443/1‑150 (MQ=255)
                                      tctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAt                                                                                                               <  1:2161150/150‑1 (MQ=255)
                                       ctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAtt                                                                                                              <  1:3374205/150‑1 (MQ=255)
                                       ctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACACAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAtt                                                                                                              <  2:4245784/150‑1 (MQ=255)
                                       ctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAtt                                                                                                              >  1:1884571/1‑150 (MQ=255)
                                         ttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAtttt                                                                                                            <  2:3768676/150‑1 (MQ=255)
                                         ttcGCAATGGCACCGCAAAAACTGCCCGTTGCTGAATTATCTTTGCCTCTTCAGGCGCAACATCTTGCATTTCTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAtttt                                                                                                            <  2:693399/150‑1 (MQ=255)
                                           cGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAc                                                                                                          <  2:1415851/150‑1 (MQ=255)
                                            gCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAcc                                                                                                         <  2:2663287/150‑1 (MQ=255)
                                            gCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAcc                                                                                                         <  2:723545/150‑1 (MQ=255)
                                                   cACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACACTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAg                                                                                                  <  1:3074696/150‑1 (MQ=255)
                                                     ccGCAAACACGGCCCGGTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATGACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCa                                                                                                <  2:2712218/150‑1 (MQ=255)
                                                         aaaCACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAAt                                                                                            >  2:2595003/1‑150 (MQ=255)
                                                         aaaCACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAAt                                                                                            <  2:3499122/150‑1 (MQ=255)
                                                         aaaCACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAAt                                                                                            <  2:3499223/150‑1 (MQ=255)
                                                          aaCACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAAtt                                                                                           <  2:1495103/150‑1 (MQ=255)
                                                          aaCACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAAtt                                                                                           >  1:2663193/1‑150 (MQ=255)
                                                           acacTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTc                                                                                          <  2:1935027/150‑1 (MQ=255)
                                                              cTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGc                                                                                       <  1:249395/150‑1 (MQ=255)
                                                               tGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCa                                                                                      >  2:2286669/1‑150 (MQ=255)
                                                               tGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCa                                                                                      <  1:4137076/150‑1 (MQ=255)
                                                                gCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCTCATTTTACCCCACAAGCAAAATTCAGCAc                                                                                     <  2:794447/150‑1 (MQ=255)
                                                                 cccGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACg                                                                                    >  1:199043/1‑150 (MQ=255)
                                                                 cccGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACg                                                                                    >  2:3607154/1‑150 (MQ=255)
                                                                 cccGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACg                                                                                    >  2:3606920/1‑150 (MQ=255)
                                                                          gaagaaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTgg                                                                           >  1:1281656/1‑150 (MQ=255)
                                                                           aagaaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGa                                                                          >  1:724283/1‑150 (MQ=255)
                                                                            agaaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGAt                                                                         >  2:100541/1‑150 (MQ=255)
                                                                              aaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATaa                                                                       >  2:2307597/1‑150 (MQ=255)
                                                                               aTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAg                                                                      >  2:1598732/1‑150 (MQ=255)
                                                                                   gtGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAATAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAg                                                                  <  2:2867545/149‑1 (MQ=255)
                                                                                          ttCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAg                                                           <  1:2213830/150‑1 (MQ=255)
                                                                                              ggCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCAt                                                       <  1:563623/150‑1 (MQ=255)
                                                                                              ggCGCAACATCGTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCAt                                                       <  1:2848491/150‑1 (MQ=255)
                                                                                                 gcAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCg                                                    >  1:496670/1‑150 (MQ=255)
                                                                                                   aaCATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTg                                                  <  1:846434/150‑1 (MQ=255)
                                                                                                   aaCATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTg                                                  <  2:4535376/150‑1 (MQ=255)
                                                                                                    aCATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGc                                                 >  2:802854/1‑150 (MQ=255)
                                                                                                     cctCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCt                                                <  2:3325296/148‑1 (MQ=255)
                                                                                                     cATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCt                                                <  2:1855503/150‑1 (MQ=255)
                                                                                                            cATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTc                                         <  2:1290471/150‑1 (MQ=255)
                                                                                                              ttACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCaa                                       <  2:1912666/150‑1 (MQ=255)
                                                                                                              ttACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCaa                                       <  1:4375861/150‑1 (MQ=255)
                                                                                                              ttACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCaa                                       <  1:2394451/150‑1 (MQ=255)
                                                                                                               tACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAAc                                      <  2:1733619/150‑1 (MQ=255)
                                                                                                                aCTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACg                                     <  2:430808/150‑1 (MQ=255)
                                                                                                                  tGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGcc                                   >  1:369909/1‑150 (MQ=255)
                                                                                                                    gATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTc                                 >  1:329293/1‑150 (MQ=255)
                                                                                                                    gATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTc                                 >  1:378333/1‑150 (MQ=255)
                                                                                                                    gATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTc                                 >  1:1290074/1‑150 (MQ=255)
                                                                                                                      tACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCAc                               >  1:3791059/1‑150 (MQ=255)
                                                                                                                      tACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCAc                               >  1:1739430/1‑150 (MQ=255)
                                                                                                                       aCGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCAcc                              >  1:461175/1‑150 (MQ=255)
                                                                                                                         gCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGTAAAGGGTTATCGCATTGATTATGCGCATTTTTCCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCATGCGCAAGGCATCCGTGCTGGCCCTCTACGCCTCACCTa                            <  2:3045637/150‑1 (MQ=255)
                                                                                                                         gCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTa                            >  1:910069/1‑150 (MQ=255)
                                                                                                                           tCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTaac                          >  1:4514299/1‑150 (MQ=255)
                                                                                                                             aGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTaacaa                        <  1:2544326/150‑1 (MQ=255)
                                                                                                                                tcgcTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAAcaacaa                     >  1:763975/2‑150 (MQ=255)
                                                                                                                                tcgcTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAAcaacaa                     <  2:3743099/149‑1 (MQ=255)
                                                                                                                                 cgcTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAAt                    >  1:1195668/1‑150 (MQ=255)
                                                                                                                                 cgcTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTAAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAAt                    <  1:656192/150‑1 (MQ=255)
                                                                                                                                  gctgctGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATa                   >  1:4264477/1‑150 (MQ=255)
                                                                                                                                    tgctgACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATaaa                 >  1:1021030/1‑150 (MQ=255)
                                                                                                                                    tgctgACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATaaa                 >  2:3074340/1‑150 (MQ=255)
                                                                                                                                     gctgACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAAc                >  1:3817501/1‑150 (MQ=255)
                                                                                                                                     gctgACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAAc                >  1:3817161/1‑150 (MQ=255)
                                                                                                                                       tgACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCt              >  1:3869540/1‑150 (MQ=255)
                                                                                                                                       tgACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCt              >  1:3869493/1‑150 (MQ=255)
                                                                                                                                       tgACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCt              >  2:2272574/1‑150 (MQ=255)
                                                                                                                                       tgACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCt              <  1:3568962/150‑1 (MQ=255)
                                                                                                                                        gACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCtt             <  1:2738827/150‑1 (MQ=255)
                                                                                                                                          cccAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTa           >  1:589107/1‑150 (MQ=255)
                                                                                                                                          cccAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTa           >  1:2875108/1‑150 (MQ=255)
                                                                                                                                            cAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTACt         >  1:1316329/1‑150 (MQ=255)
                                                                                                                                             aGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTACtt        >  2:2622319/1‑150 (MQ=255)
                                                                                                                                                  ggCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTACTTCAttt   <  2:1157085/150‑1 (MQ=255)
                                                                                                                                                   gCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTACTTCAtttt  <  2:2820025/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTACTTCATTTT  >  CP009273/104433‑104729

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: