Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 104,582 A→T 35.6% 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,5820AT35.6% 318.8 / inf 292T131S (ACG→TCG) secMregulator of secA translation
Reads supporting (aligned to +/- strand):  ref base A (101/87);  new base T (56/48);  total (157/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

Only 100 of 295 total aligned reads displayed.

GGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTACTTCATTTTA  >  CP009273/104433‑104730
                                                                                                                                                     |                                                                                                                                                    
ggCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCa                                                                                                                                                      <  2:2891458/150‑1 (MQ=255)
 gCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGTCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCAc                                                                                                                                                     <  1:3326403/150‑1 (MQ=255)
      aaCATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGt                                                                                                                                                >  2:1416626/1‑150 (MQ=255)
        cATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCt                                                                                                                                              <  2:223069/150‑1 (MQ=255)
         aTGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTg                                                                                                                                             >  1:630044/1‑150 (MQ=255)
             cATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGaaaa                                                                                                                                         >  1:1681552/1‑150 (MQ=255)
             cATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGaaaa                                                                                                                                         <  1:4655330/150‑1 (MQ=255)
              aTTCGCTCGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGACCGTTGCTGAAGAATCTTTGCCTCTACAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAg                                                                                                                                        <  2:4465079/150‑1 (MQ=255)
              aTTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAg                                                                                                                                        >  2:3483397/1‑150 (MQ=255)
               ttCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAgg                                                                                                                                       <  2:4764277/150‑1 (MQ=255)
               ttCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAgg                                                                                                                                       <  2:2834373/150‑1 (MQ=255)
                  gCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGtt                                                                                                                                    >  2:2528931/1‑150 (MQ=255)
                   cACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTa                                                                                                                                   >  1:4689546/1‑150 (MQ=255)
                    aCGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTAt                                                                                                                                  <  1:4447799/150‑1 (MQ=255)
                     cGGTAAGCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATc                                                                                                                                 <  2:4038408/150‑1 (MQ=255)
                        tAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCa                                                                                                                              >  1:1559724/1‑150 (MQ=255)
                         aaTCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCAt                                                                                                                             <  2:896239/150‑1 (MQ=255)
                           tCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTg                                                                                                                           <  1:2128127/150‑1 (MQ=255)
                           tCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTg                                                                                                                           <  1:4283889/150‑1 (MQ=255)
                           tCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTg                                                                                                                           <  1:1881038/150‑1 (MQ=255)
                           tCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTg                                                                                                                           <  2:3575319/150‑1 (MQ=255)
                           tCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTg                                                                                                                           <  2:4581712/150‑1 (MQ=255)
                           tCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTg                                                                                                                           <  2:25768/150‑1 (MQ=255)
                             cGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGAt                                                                                                                         <  1:2980317/150‑1 (MQ=255)
                             cGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGAt                                                                                                                         <  2:408867/150‑1 (MQ=255)
                              gTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGAtt                                                                                                                        >  2:2784674/1‑150 (MQ=255)
                              gTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGAtt                                                                                                                        <  2:3078147/150‑1 (MQ=255)
                                cATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTAt                                                                                                                      >  1:1039884/1‑150 (MQ=255)
                                      tctttcGCAATGGCACCGCATACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAt                                                                                                                <  2:197633/150‑1 (MQ=255)
                                        tttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAttt                                                                                                              >  1:3233439/1‑150 (MQ=255)
                                        tttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAttt                                                                                                              <  1:2169707/150‑1 (MQ=255)
                                        tttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTGTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAttt                                                                                                              <  2:3260400/150‑1 (MQ=255)
                                              aaTGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGAGTATGCGCATTTTAcccc                                                                                                        >  1:1213625/1‑150 (MQ=255)
                                              aaTGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAcccc                                                                                                        >  1:1401061/1‑150 (MQ=255)
                                              aaTGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAcccc                                                                                                        >  2:907704/1‑150 (MQ=255)
                                               aTGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCca                                                                                                       <  2:785961/150‑1 (MQ=255)
                                                 ggCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCcaca                                                                                                     >  1:4630772/1‑150 (MQ=255)
                                                 ggCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCcaca                                                                                                     >  1:2230096/1‑150 (MQ=255)
                                                 ggCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCAGTGAGTCTGCGCATGTTACCCcaca                                                                                                     >  1:2676022/1‑150 (MQ=255)
                                                      cGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCaa                                                                                                <  1:646463/150‑1 (MQ=255)
                                                      cGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCaa                                                                                                >  1:720655/1‑150 (MQ=255)
                                                       gCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCaaa                                                                                               >  1:2010245/1‑150 (MQ=255)
                                                        cAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCaaaa                                                                                              >  1:4020315/1‑150 (MQ=255)
                                                         aaaCACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTCTGCGCATTTTACCCCACAAGCAAAAt                                                                                             >  1:4228101/1‑150 (MQ=255)
                                                          aaCACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAAtt                                                                                            <  1:4792438/150‑1 (MQ=255)
                                                            cacTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCa                                                                                          <  2:120720/150‑1 (MQ=255)
                                                             acTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAg                                                                                         >  1:841832/1‑150 (MQ=255)
                                                              cTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGc                                                                                        >  2:652955/1‑150 (MQ=255)
                                                                   cGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGcc                                                                                   <  1:4523376/150‑1 (MQ=255)
                                                                   cGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGcc                                                                                   >  2:4675845/1‑150 (MQ=255)
                                                                       gCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTc                                                                               >  2:3704931/1‑150 (MQ=255)
                                                                        cTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCt                                                                              <  1:1455887/150‑1 (MQ=255)
                                                                          gaagaaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTgg                                                                            >  2:3795961/1‑150 (MQ=255)
                                                                            agaaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGAt                                                                          >  1:979506/1‑150 (MQ=255)
                                                                              aaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATaa                                                                        >  2:551611/1‑150 (MQ=255)
                                                                              aaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATaa                                                                        <  1:1062998/150‑1 (MQ=255)
                                                                              aaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATaa                                                                        >  2:167794/1‑150 (MQ=255)
                                                                              aaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATaa                                                                        >  2:20705/1‑150 (MQ=255)
                                                                              aaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATaa                                                                        >  2:21651/1‑150 (MQ=255)
                                                                               aTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAg                                                                       >  2:4503363/1‑150 (MQ=255)
                                                                               aTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAACATTCAGCACGCCCGTCTGGATAAg                                                                       >  2:3158318/1‑150 (MQ=255)
                                                                               aTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAg                                                                       >  1:3433689/1‑150 (MQ=255)
                                                                                tCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGc                                                                      >  1:1760163/1‑150 (MQ=255)
                                                                                tCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGc                                                                      >  1:4127006/1‑150 (MQ=255)
                                                                                   ttGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAg                                                                   >  2:3352887/1‑150 (MQ=255)
                                                                                      ccTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGgcg                                                                >  1:3988476/1‑150 (MQ=255)
                                                                                            cAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGc                                                          >  1:969887/1‑150 (MQ=255)
                                                                                             aGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCa                                                         <  1:2784674/150‑1 (MQ=255)
                                                                                             aGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCa                                                         <  1:3377456/150‑1 (MQ=255)
                                                                                             aGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCa                                                         >  2:2463834/1‑150 (MQ=255)
                                                                                                cgcAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATcc                                                      >  2:2825539/1‑150 (MQ=255)
                                                                                                  cAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGt                                                    >  1:3638/1‑150 (MQ=255)
                                                                                                      aTCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTg                                                >  1:2311635/1‑150 (MQ=255)
                                                                                                        cTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGc                                              >  1:2695611/1‑150 (MQ=255)
                                                                                                        cTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGc                                              >  1:4210434/1‑150 (MQ=255)
                                                                                                           gCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCt                                           <  1:3155148/150‑1 (MQ=255)
                                                                                                           gCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCt                                           >  1:1627859/1‑150 (MQ=255)
                                                                                                            cATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTc                                          >  2:1489176/1‑150 (MQ=255)
                                                                                                             aTTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCa                                         >  2:2685515/1‑150 (MQ=255)
                                                                                                              ttACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCaa                                        >  1:793426/1‑150 (MQ=255)
                                                                                                               tACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAAc                                       >  1:3923136/1‑150 (MQ=255)
                                                                                                               tACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAAc                                       <  2:4347972/150‑1 (MQ=255)
                                                                                                                aCTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACg                                      >  1:3964169/1‑150 (MQ=255)
                                                                                                                  tGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGcc                                    <  1:554432/150‑1 (MQ=255)
                                                                                                                  tGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGcc                                    >  2:116165/1‑150 (MQ=255)
                                                                                                                   ggATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCt                                   >  2:1839250/1‑150 (MQ=255)
                                                                                                                      tACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCAc                                >  2:3330168/1‑150 (MQ=255)
                                                                                                                       aCGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCAcc                               <  2:3852966/150‑1 (MQ=255)
                                                                                                                        cGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCt                              <  1:2591559/150‑1 (MQ=255)
                                                                                                                        cGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCt                              <  2:4198166/150‑1 (MQ=255)
                                                                                                                         gCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTa                             <  1:1597198/150‑1 (MQ=255)
                                                                                                                         gCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTa                             <  1:4403558/150‑1 (MQ=255)
                                                                                                                                tcgcTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAAcaacaa                      >  2:3285862/2‑150 (MQ=255)
                                                                                                                                    tgctgACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATaaa                  >  2:1267012/1‑150 (MQ=255)
                                                                                                                                       tgACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCt               >  2:2368380/1‑150 (MQ=255)
                                                                                                                                          cccAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTa            <  2:3791119/150‑1 (MQ=255)
                                                                                                                                          cccAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTa            <  2:1578939/150‑1 (MQ=255)
                                                                                                                                           ccAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTAc           <  2:3987510/150‑1 (MQ=255)
                                                                                                                                             aGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTACtt         >  1:3745736/1‑150 (MQ=255)
                                                                                                                                                   gCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTACTTCAtttt   <  1:2528931/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAACAATAAACCTTTACTTCATTTTA  >  CP009273/104433‑104730

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: