Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,348,772 G→T 100% P2T (CCT→ACT) 
A319D (GCC→GAC) 
sapC ←
sapB ←
antimicrobial peptide transport ABC transporter permease
antimicrobial peptide transport ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,348,7720GT97.0% 1006.9 / 16.5 319P2T (CCT→ACT) 
A319D (GCC→GAC) 
sapC
sapB
antimicrobial peptide transport ABC transporter permease
antimicrobial peptide transport ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (2/8);  new base T (124/185);  total (126/193)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 319 total aligned reads displayed.

GCGCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAG  >  CP009273/1348625‑1348918
                                                                                                                                                   |                                                                                                                                                  
  gcAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATa                                                                                                                                                <  2:224380/150‑1 (MQ=255)
    aaaCCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATAcc                                                                                                                                              >  2:660771/1‑150 (MQ=255)
      aCCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCAt                                                                                                                                            <  1:4861614/150‑1 (MQ=255)
           ccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCtt                                                                                                                                       <  2:1603668/150‑1 (MQ=255)
            cgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTa                                                                                                                                      >  1:4708292/1‑150 (MQ=255)
            cgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTa                                                                                                                                      <  1:722533/150‑1 (MQ=255)
            cgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTa                                                                                                                                      <  2:4585319/150‑1 (MQ=255)
             gccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTAt                                                                                                                                     <  2:841513/150‑1 (MQ=255)
               cAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGt                                                                                                                                   >  1:3153411/1‑150 (MQ=255)
                aaaaaTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGtt                                                                                                                                  <  2:1540575/150‑1 (MQ=255)
                 aaaaTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGttt                                                                                                                                 <  2:289387/150‑1 (MQ=255)
                  aaaTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTc                                                                                                                                >  1:2664088/1‑150 (MQ=255)
                    aTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCag                                                                                                                              <  2:4463848/150‑1 (MQ=255)
                        acaGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCagagtg                                                                                                                          <  1:4218311/150‑3 (MQ=255)
                          aGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGtt                                                                                                                        <  1:1445577/150‑1 (MQ=255)
                          aGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGtt                                                                                                                        >  1:2507421/1‑150 (MQ=255)
                           gTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTa                                                                                                                       <  1:4185078/150‑1 (MQ=255)
                            tACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAg                                                                                                                      >  2:1567034/1‑150 (MQ=255)
                            tACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCGATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAg                                                                                                                      >  1:1856346/1‑150 (MQ=255)
                               agccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCa                                                                                                                   >  2:3323382/1‑150 (MQ=255)
                               agccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCa                                                                                                                   >  2:2466613/1‑150 (MQ=255)
                                gccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCAt                                                                                                                  <  1:1935920/150‑1 (MQ=255)
                                  cagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATgg                                                                                                                >  2:1469004/1‑150 (MQ=255)
                                  cagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATgg                                                                                                                <  1:3623821/150‑1 (MQ=255)
                                  cagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATgg                                                                                                                <  1:199845/150‑1 (MQ=255)
                                   agccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGc                                                                                                               >  1:976241/1‑150 (MQ=255)
                                      cccGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAcc                                                                                                            >  2:3266455/1‑150 (MQ=255)
                                      cccGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAcc                                                                                                            >  1:3253593/1‑150 (MQ=255)
                                         gcgcAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCaa                                                                                                         >  2:717944/1‑150 (MQ=255)
                                          cgcAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCaaa                                                                                                        >  1:3341689/1‑150 (MQ=255)
                                            cAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAat                                                                                                      >  2:2585916/1‑150 (MQ=255)
                                             aGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAata                                                                                                     >  1:480183/1‑150 (MQ=255)
                                              gCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAatat                                                                                                    <  2:4400946/150‑1 (MQ=255)
                                                cGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCa                                                                                                  >  2:4500776/1‑150 (MQ=255)
                                                  tACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCATAGGGTTAGCCATGGCACCCAAAATATCAGa                                                                                                >  1:1523911/1‑150 (MQ=255)
                                                  tACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGa                                                                                                <  1:5085543/150‑1 (MQ=255)
                                                    cAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGaaa                                                                                              >  1:1560682/1‑150 (MQ=255)
                                                    cAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGaaa                                                                                              >  2:4103631/1‑150 (MQ=255)
                                                    cAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGaaa                                                                                              <  1:1725892/150‑1 (MQ=255)
                                                        ccGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCAc                                                                                          <  1:3964147/150‑1 (MQ=255)
                                                         cGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACg                                                                                         <  1:1187846/150‑1 (MQ=255)
                                                           aCCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGGCATGGCACCCAAAATATCAGAAATCACGtt                                                                                       <  2:1645174/150‑1 (MQ=255)
                                                           aCCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGtt                                                                                       <  2:935192/150‑1 (MQ=255)
                                                            ccATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                                                      <  1:1659342/150‑1 (MQ=255)
                                                            ccATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGACATGGCACCCAAAATATCAGAAATCACGTTa                                                                                      <  2:3985021/150‑1 (MQ=255)
                                                             cATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTaa                                                                                     <  2:4904953/150‑1 (MQ=255)
                                                              aTCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAAc                                                                                    <  2:1286023/150‑1 (MQ=255)
                                                               tCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACa                                                                                   >  1:79301/1‑150 (MQ=255)
                                                                 gCAGAGGCTTCAATATAAAATTTGCGCCATGAGGTACGCAGAGTGCCCGGCGGTCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAACCAAAATATCAGAAATCACGTTAACaat                                                                                 <  1:2539230/150‑1 (MQ=255)
                                                                   agagGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaataa                                                                               <  2:4700488/150‑1 (MQ=255)
                                                                     agGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATg                                                                             >  1:3580727/1‑150 (MQ=255)
                                                                         gTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCa                                                                         >  2:691768/1‑150 (MQ=255)
                                                                         gTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCa                                                                         <  1:1213604/150‑1 (MQ=255)
                                                                          tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAg                                                                        <  1:4486192/150‑1 (MQ=255)
                                                                          tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAg                                                                        <  2:90188/150‑1 (MQ=255)
                                                                           cACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGt                                                                       <  1:1619389/150‑1 (MQ=255)
                                                                            aCTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTg                                                                      >  1:406909/1‑150 (MQ=255)
                                                                             cTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGa                                                                     <  1:1259108/150‑1 (MQ=255)
                                                                              tataAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAg                                                                    <  1:3814932/150‑1 (MQ=255)
                                                                              tataAAATATGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCAATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATAACGTTAACAATAATGACCAGTGAg                                                                    >  1:5162212/1‑150 (MQ=255)
                                                                                taAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGcc                                                                  >  1:4595393/1‑150 (MQ=255)
                                                                                 aaaaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCca                                                                 <  2:3594646/150‑1 (MQ=255)
                                                                                 aaaaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCca                                                                 <  2:303825/150‑1 (MQ=255)
                                                                                   aaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCcaca                                                               <  2:3756504/150‑1 (MQ=255)
                                                                                         cggcAGGCGGTACGCAGCGTGCGCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACAATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCAt                                                         <  2:2653518/147‑1 (MQ=255)
                                                                                            cAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCAc                                                      >  2:750956/1‑150 (MQ=255)
                                                                                               gCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTcc                                                   <  1:2265542/150‑1 (MQ=255)
                                                                                                 ggTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcgg                                                 >  2:4264251/1‑150 (MQ=255)
                                                                                                     cGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGa                                             <  1:33663/150‑1 (MQ=255)
                                                                                                      gCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaa                                            <  1:558948/150‑1 (MQ=255)
                                                                                                      gCAGCGTGCCCGGAGGTCGCTTTTAGCTGTATACTCTATCGTAATTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGTAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaa                                            <  2:1195651/150‑1 (MQ=255)
                                                                                                       cAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCGTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaaa                                           <  1:1910475/150‑1 (MQ=255)
                                                                                                           gTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGc                                       <  1:2179787/150‑1 (MQ=255)
                                                                                                            tGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCt                                      <  2:2802369/150‑1 (MQ=255)
                                                                                                            tGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCt                                      <  1:3219750/150‑1 (MQ=255)
                                                                                                            tGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCt                                      <  2:516624/150‑1 (MQ=255)
                                                                                                                   cgggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAAt                               >  1:3723646/1‑150 (MQ=255)
                                                                                                                   cgggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAAt                               >  1:436307/1‑150 (MQ=255)
                                                                                                                   cgggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACACCGGCGGAAATGGCTGCGTAAt                               >  1:1350386/1‑150 (MQ=255)
                                                                                                                    gggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATc                              <  1:2873180/150‑1 (MQ=255)
                                                                                                                    gggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATc                              >  1:1146709/1‑150 (MQ=255)
                                                                                                                     ggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATcc                             <  1:660747/150‑1 (MQ=255)
                                                                                                                     ggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAAGAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATcc                             <  1:3128431/150‑1 (MQ=255)
                                                                                                                       cgcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctg                           >  2:4879435/1‑150 (MQ=255)
                                                                                                                        gcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgc                          <  1:5136596/150‑1 (MQ=255)
                                                                                                                         cTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgct                         <  2:254386/150‑1 (MQ=255)
                                                                                                                         cTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgct                         >  1:398930/1‑150 (MQ=255)
                                                                                                                          ttttCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgctg                        >  2:1121748/1‑150 (MQ=255)
                                                                                                                          ttttCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgctg                        >  2:1393613/1‑150 (MQ=255)
                                                                                                                          ttttCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgctg                        <  2:4797138/150‑1 (MQ=255)
                                                                                                                            ttCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGc                      <  1:2596453/150‑1 (MQ=255)
                                                                                                                            ttCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGc                      <  2:3783804/150‑1 (MQ=255)
                                                                                                                            ttCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGc                      >  1:893502/1‑150 (MQ=255)
                                                                                                                                cTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGAt                  >  2:4738265/1‑150 (MQ=255)
                                                                                                                                   tataCGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGc               <  2:2954746/150‑1 (MQ=255)
                                                                                                                                   tataCGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGc               >  2:2379503/1‑150 (MQ=255)
                                                                                                                                   tataCGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGc               >  1:4123678/1‑150 (MQ=255)
                                                                                                                                   tataCGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGc               >  1:1002434/1‑150 (MQ=255)
                                                                                                                                                aaGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAg  >  2:585473/1‑150 (MQ=255)
                                                                                                                                                aaGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAg  <  2:1916781/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAG  >  CP009273/1348625‑1348918

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: