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,7720GT99.1% 756.9 / ‑2.1 233P2T (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/0);  new base T (92/139);  total (94/139)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 234 total aligned reads displayed.

GCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCG  >  CP009273/1348627‑1348920
                                                                                                                                                 |                                                                                                                                                    
gcAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCAt                                                                                                                                                   >  1:752988/1‑149 (MQ=255)
gcAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCAt                                                                                                                                                   <  2:2622762/149‑1 (MQ=255)
 cAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATa                                                                                                                                                  <  2:1221621/149‑1 (MQ=255)
   aaCCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATAcc                                                                                                                                                >  2:1016769/1‑149 (MQ=255)
    aCCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCa                                                                                                                                               >  2:1244588/1‑149 (MQ=255)
      cAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCAtt                                                                                                                                             >  2:2363988/1‑149 (MQ=255)
       aGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTc                                                                                                                                            <  1:1066987/149‑1 (MQ=255)
       aGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTc                                                                                                                                            <  1:3128826/149‑1 (MQ=255)
           gccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTa                                                                                                                                        <  1:1405516/149‑1 (MQ=255)
           gccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTa                                                                                                                                        <  1:2342677/149‑1 (MQ=255)
               aaaaTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGtt                                                                                                                                    <  1:4082075/149‑1 (MQ=255)
                aaaTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGttt                                                                                                                                   <  1:3043539/149‑1 (MQ=255)
                 aaTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTc                                                                                                                                  <  1:1716130/149‑1 (MQ=255)
                 aaTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTc                                                                                                                                  <  1:1938935/149‑1 (MQ=255)
                  aTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCAATCCTTATGTTTCa                                                                                                                                 >  1:3995845/1‑149 (MQ=255)
                      acaGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAgg                                                                                                                             <  2:2498918/149‑1 (MQ=255)
                        aGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGt                                                                                                                           <  2:3844979/149‑1 (MQ=255)
                        aGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGt                                                                                                                           <  2:1510062/149‑1 (MQ=255)
                         gTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGtt                                                                                                                          <  2:857557/149‑1 (MQ=255)
                         gTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGtt                                                                                                                          <  1:2533316/149‑1 (MQ=255)
                             agccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGcc                                                                                                                      <  2:3351989/149‑1 (MQ=255)
                             agccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGcc                                                                                                                      <  1:1350470/149‑1 (MQ=255)
                              gccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCa                                                                                                                     <  1:1126522/149‑1 (MQ=255)
                              gccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCa                                                                                                                     <  1:1006420/149‑1 (MQ=255)
                              gccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCa                                                                                                                     <  1:1497324/149‑1 (MQ=255)
                               ccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCAt                                                                                                                    <  1:2096304/149‑1 (MQ=255)
                                cagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATg                                                                                                                   >  2:2178448/1‑149 (MQ=255)
                                    cccGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTAAACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAc                                                                                                               >  2:2735330/1‑149 (MQ=255)
                                    cccGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATGTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAc                                                                                                               <  2:358647/149‑1 (MQ=255)
                                      cgcgCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAccc                                                                                                             <  1:206020/149‑1 (MQ=255)
                                       gcgcAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCa                                                                                                            >  1:94011/1‑149 (MQ=255)
                                         gcAGTCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCaaa                                                                                                          <  1:2366786/149‑1 (MQ=255)
                                            gCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAata                                                                                                       >  2:3723403/1‑149 (MQ=255)
                                            gCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAata                                                                                                       <  1:3548968/149‑1 (MQ=255)
                                            gCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAata                                                                                                       <  2:1871940/149‑1 (MQ=255)
                                             ccGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAatat                                                                                                      <  1:374140/149‑1 (MQ=255)
                                             ccGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAatat                                                                                                      <  1:3975780/149‑1 (MQ=255)
                                               gTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCa                                                                                                    >  2:780529/1‑149 (MQ=255)
                                                   aGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGaaa                                                                                                <  2:3893723/149‑1 (MQ=255)
                                                   aGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGaaa                                                                                                <  2:733540/149‑1 (MQ=255)
                                                    ggCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAAt                                                                                               <  2:1710161/149‑1 (MQ=255)
                                                        gACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATAAGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCAAGGCACCCCAAATATCAGAAATCACg                                                                                           >  1:3475437/1‑149 (MQ=255)
                                                         aCCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGt                                                                                          <  2:3673556/149‑1 (MQ=255)
                                                          ccATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGtt                                                                                         >  2:656257/1‑149 (MQ=255)
                                                          ccATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTAGCTGTATAAGCTATCGTAATTCATAACATTACTTATGTTTCAGAGGGTTAGCCATGGGCCCCAAAATATCAGAAATAACGtt                                                                                         >  1:3103781/1‑149 (MQ=255)
                                                           cATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                                                        <  1:3922468/149‑1 (MQ=255)
                                                               gCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaa                                                                                    <  1:2019354/149‑1 (MQ=255)
                                                                 agagGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaata                                                                                  >  2:3212620/1‑149 (MQ=255)
                                                                  gagGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaataa                                                                                 <  1:199116/149‑1 (MQ=255)
                                                                  gagGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaataa                                                                                 <  2:510030/149‑1 (MQ=255)
                                                                  gagGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCGTATGTTTCAGAGGGTTAGCCATGGCCCCCAAAATATCAGAAAACACGTTAACaataa                                                                                 >  1:1561017/1‑149 (MQ=255)
                                                                     gCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGa                                                                              <  1:1840504/149‑1 (MQ=255)
                                                                     gCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGa                                                                              <  1:1818355/149‑1 (MQ=255)
                                                                     gCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGa                                                                              >  1:1639792/1‑149 (MQ=255)
                                                                     gCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCGTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGa                                                                              >  1:915438/1‑149 (MQ=255)
                                                                     gCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATAACATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAAAAATGa                                                                              >  1:2974683/1‑149 (MQ=255)
                                                                         cACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAg                                                                          <  2:925435/149‑1 (MQ=255)
                                                                         cACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAg                                                                          <  1:3765192/149‑1 (MQ=255)
                                                                         cACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAg                                                                          >  1:3954399/1‑149 (MQ=255)
                                                                          aCTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGt                                                                         >  1:3731585/1‑149 (MQ=255)
                                                                           cTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGTGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTg                                                                        <  1:3051718/149‑1 (MQ=255)
                                                                           cTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTg                                                                        >  1:2408373/1‑149 (MQ=255)
                                                                           cTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTg                                                                        <  1:3584849/149‑1 (MQ=255)
                                                                           cTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTg                                                                        <  1:3218809/149‑1 (MQ=255)
                                                                             ataAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAg                                                                      >  2:3402585/1‑149 (MQ=255)
                                                                              taAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATTTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGc                                                                     <  2:2940606/149‑1 (MQ=255)
                                                                              taAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGc                                                                     <  2:1550852/149‑1 (MQ=255)
                                                                                 aaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCcac                                                                  >  1:860835/1‑149 (MQ=255)
                                                                                 aaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTAGCTGTATACGCTATCGTAAGTCATACCATTCCGTATGTTTCAGAGGGTTAGCCATGGGACCCAAAATATCAGAAATAACGTTAAAAATAATGACCAGTGTGCcac                                                                  >  1:2805688/1‑149 (MQ=255)
                                                                                    ttGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCcacaca                                                               <  1:3394622/149‑1 (MQ=255)
                                                                                    gtGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCcacaca                                                               <  2:2986418/148‑1 (MQ=255)
                                                                                       cgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCa                                                            <  1:2959792/149‑1 (MQ=255)
                                                                                           aGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATAACATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAAAAATGACCAGTGAGCCACACACCATCAc                                                        >  1:457977/1‑149 (MQ=255)
                                                                                             gCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTc                                                      >  1:1655524/1‑149 (MQ=255)
                                                                                              cGGTGCGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTcc                                                     <  2:3525214/149‑1 (MQ=255)
                                                                                              cGGTACGCAGCGTGCCGGGCTTTCGATTTTCGCTGTATACGCTATCGTTAGTAATACTATTCATTATGTTTCAGAGGGTTATCCATGGCACCCAAAATATCAGATATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACtac                                                     <  2:2815540/149‑3 (MQ=255)
                                                                                              cGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTcc                                                     <  2:752988/149‑1 (MQ=255)
                                                                                               ggTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcg                                                    <  1:3315132/149‑1 (MQ=255)
                                                                                                 tACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcggc                                                  <  2:3016454/149‑1 (MQ=255)
                                                                                                  aCGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcggcg                                                 <  1:2208626/149‑1 (MQ=255)
                                                                                                   cGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcggcgg                                                <  2:602020/149‑1 (MQ=255)
                                                                                                     caacGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaa                                              >  1:2999224/4‑149 (MQ=255)
                                                                                                        cGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATg                                           <  2:1607143/149‑1 (MQ=255)
                                                                                                          tGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGc                                         <  2:3111388/149‑1 (MQ=255)
                                                                                                           gCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCt                                        <  2:3689363/149‑1 (MQ=255)
                                                                                                                   ggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATc                                >  1:2264121/1‑149 (MQ=255)
                                                                                                                    gcgcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATTATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATcc                               <  2:2231804/149‑1 (MQ=255)
                                                                                                                    gcgcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATcc                               <  1:19500/149‑1 (MQ=255)
                                                                                                                    gcgcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATcc                               >  2:3494759/1‑149 (MQ=255)
                                                                                                                       cTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgc                            <  1:741440/149‑1 (MQ=255)
                                                                                                                                gTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATg                   <  2:2051494/149‑1 (MQ=255)
                                                                                                                                 tataCGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATgg                  <  1:3609364/149‑1 (MQ=255)
                                                                                                                                      gCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGtta             >  1:3007150/1‑149 (MQ=255)
                                                                                                                                      gCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGtta             >  1:3215552/1‑149 (MQ=255)
                                                                                                                                         aTCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGttaatt          >  1:3060706/1‑149 (MQ=255)
                                                                                                                                            gTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAAc       <  2:3954399/149‑1 (MQ=255)
                                                                                                                                             tAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAAcc      >  2:1422240/1‑149 (MQ=255)
                                                                                                                                             tAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAAcc      <  1:2744570/149‑1 (MQ=255)
                                                                                                                                              aaGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCa     >  1:239640/1‑149 (MQ=255)
                                                                                                                                              aaGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCa     >  1:3277398/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCG  >  CP009273/1348627‑1348920

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: