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.5% 599.8 / ‑4.0 184P2T (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 (0/1);  new base T (61/122);  total (61/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-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 184 total aligned reads displayed.

AAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACC  >  CP009273/1348629‑1348916
                                                                                                                                               |                                                                                                                                                
  aCCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATAcc                                                                                                                                            >  2:1731713/1‑148 (MQ=255)
     aGCCGCCAAAAAGACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATACAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCAtt                                                                                                                                         <  2:2643598/148‑1 (MQ=255)
      gccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTc                                                                                                                                        >  2:1613453/1‑148 (MQ=255)
      gccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTc                                                                                                                                        <  1:2501626/148‑1 (MQ=255)
      gccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTc                                                                                                                                        <  1:1788320/148‑1 (MQ=255)
       ccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTcc                                                                                                                                       <  2:3095114/148‑1 (MQ=255)
       ccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTcc                                                                                                                                       <  2:3095142/148‑1 (MQ=255)
       ccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTcc                                                                                                                                       <  2:3795398/148‑1 (MQ=255)
        cgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCt                                                                                                                                      <  1:3789409/148‑1 (MQ=255)
        cgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCt                                                                                                                                      <  1:1879474/148‑1 (MQ=255)
         gccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCtt                                                                                                                                     <  2:449534/148‑1 (MQ=255)
         gccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCtt                                                                                                                                     <  1:145519/148‑1 (MQ=255)
         gccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCtt                                                                                                                                     <  2:2063222/148‑1 (MQ=255)
          ccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTa                                                                                                                                    >  2:2274436/1‑148 (MQ=255)
             aaaaTACACAGTACAGCCAGCCCCTCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGt                                                                                                                                 <  2:2547733/148‑1 (MQ=255)
                aTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTc                                                                                                                              >  1:308900/1‑148 (MQ=255)
                aTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTc                                                                                                                              <  1:1554931/148‑1 (MQ=255)
                aTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTc                                                                                                                              <  1:3453177/148‑1 (MQ=255)
                  acacAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCag                                                                                                                            >  1:3188522/1‑148 (MQ=255)
                  acacAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCag                                                                                                                            >  1:1363285/1‑148 (MQ=255)
                   cacaGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCCTTCCTGATGTTTCaga                                                                                                                           >  2:2709342/1‑148 (MQ=255)
                   cacaGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCaga                                                                                                                           <  1:2521079/148‑1 (MQ=255)
                   cacaGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCGTATGTTTCaga                                                                                                                           >  1:2632061/1‑148 (MQ=255)
                   cacaGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCGATCGTAAGTCATACCATTCCTGATGTTTCaga                                                                                                                           >  2:2875907/1‑148 (MQ=255)
                    acaGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCagag                                                                                                                          <  1:485976/148‑1 (MQ=255)
                    acaGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCagag                                                                                                                          <  2:254587/148‑1 (MQ=255)
                    acaGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCagag                                                                                                                          <  1:485899/148‑1 (MQ=255)
                         aCAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTa                                                                                                                     <  2:691544/148‑1 (MQ=255)
                           agccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCaga                                                                                                                           <  1:3572640/140‑1 (MQ=255)
                            gccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGcc                                                                                                                  <  2:1882202/148‑1 (MQ=255)
                            gccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGcc                                                                                                                  <  2:977450/148‑1 (MQ=255)
                             ccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCa                                                                                                                 <  2:1580125/148‑1 (MQ=255)
                             ccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCa                                                                                                                 <  1:220556/148‑1 (MQ=255)
                                gccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATgg                                                                                                              <  2:1189590/148‑1 (MQ=255)
                                gccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATgg                                                                                                              >  2:227324/1‑148 (MQ=255)
                                gccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATgg                                                                                                              <  2:1189016/148‑1 (MQ=255)
                                 ccccGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATAAGCTATCGGAAGTCATACCAATCCTTATGTTTCAGAGGGTTAGCCAAGGc                                                                                                             >  1:34683/1‑148 (MQ=255)
                                    cgcgCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAcc                                                                                                          <  1:2419336/148‑1 (MQ=255)
                                        cAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCaaa                                                                                                      <  1:312346/148‑1 (MQ=255)
                                          gCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCGTATGTTTAAGAGGGTTAGCCAAGGCACCCaaaaa                                                                                                    >  1:2605268/1‑147 (MQ=255)
                                           ccGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAata                                                                                                   >  2:2944484/1‑148 (MQ=255)
                                                 aGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGaa                                                                                             >  2:348980/1‑148 (MQ=255)
                                                   gCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAAt                                                                                           >  2:290047/1‑148 (MQ=255)
                                                    ccGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATc                                                                                          <  1:973708/148‑1 (MQ=255)
                                                       aCCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACg                                                                                       <  1:106563/148‑1 (MQ=255)
                                                       aCCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCAGTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACg                                                                                       >  1:820670/1‑148 (MQ=255)
                                                          aTCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGACATCACGTTa                                                                                    >  2:578372/1‑148 (MQ=255)
                                                          aTCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                                                    <  1:1817263/148‑1 (MQ=255)
                                                           tCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTaa                                                                                   >  1:2735503/1‑148 (MQ=255)
                                                           tCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCGATCGTAAGTCCAACCATTCCTTATGTTTCAGAGGGTTAGCCCTGGCACCCCAAATATCAGAAAACACGTTaa                                                                                   >  1:376516/1‑148 (MQ=255)
                                                              cAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaa                                                                                <  2:3454776/148‑1 (MQ=255)
                                                                gagGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCTTACGTTTTGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaata                                                                              <  1:3747094/148‑1 (MQ=255)
                                                                gagGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaata                                                                              <  2:660321/148‑1 (MQ=255)
                                                                gagGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaata                                                                              <  1:3475991/148‑1 (MQ=255)
                                                                 agGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaataa                                                                             >  2:303413/1‑148 (MQ=255)
                                                                 agGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaataa                                                                             <  2:2871910/148‑1 (MQ=255)
                                                                 agGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaataa                                                                             <  1:1277899/148‑1 (MQ=255)
                                                                 agGCGTCACTATAAAATTTGCGCCAGGCAGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaataa                                                                             <  1:691653/148‑1 (MQ=255)
                                                                   gCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCGTATGTGTCAGAGGGTTAGCCATGGCACCCACAATATCAGAAATCACGTTAACAATAATg                                                                           >  2:3303969/1‑148 (MQ=255)
                                                                      tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGAcc                                                                        <  2:491234/148‑1 (MQ=255)
                                                                      tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATAAGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGAcc                                                                        <  2:1244712/148‑1 (MQ=255)
                                                                          tataAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTg                                                                    >  1:2804785/1‑148 (MQ=255)
                                                                          tataAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTg                                                                    >  1:889628/1‑148 (MQ=255)
                                                                            taAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCAATCCTTATGTTTCAGAGGGTTAGCCAAGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAg                                                                  >  1:2966897/1‑148 (MQ=255)
                                                                               aaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCca                                                               <  1:3283526/148‑1 (MQ=255)
                                                                               aaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCca                                                               <  1:3283651/148‑1 (MQ=255)
                                                                               aaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCca                                                               <  2:356056/148‑1 (MQ=255)
                                                                                aTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCcac                                                              >  2:594886/1‑148 (MQ=255)
                                                                                  ttGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCcacac                                                            <  2:3235020/148‑1 (MQ=255)
                                                                                   tGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCcacaca                                                           >  1:1678102/1‑148 (MQ=255)
                                                                                         aGGCTGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATAGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCTCACACCATCa                                                     <  2:2622498/148‑1 (MQ=255)
                                                                                          ggCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCAc                                                    >  2:3621096/1‑148 (MQ=255)
                                                                                            cGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTc                                                  <  1:3418261/148‑1 (MQ=255)
                                                                                               tACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcgg                                               <  1:3468849/148‑1 (MQ=255)
                                                                                               tACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcgg                                               >  2:3599440/1‑148 (MQ=255)
                                                                                               tACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGGTAACAATAATGACCAGTGAGCCACACACCATCACTCcgg                                               >  1:1819344/1‑148 (MQ=255)
                                                                                                    aGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCATTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACACTAATGACCAGTGAGCCACACACCATCACTCCGGCGGaa                                          <  1:2896650/148‑1 (MQ=255)
                                                                                                      cGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAAt                                        <  2:1149175/148‑1 (MQ=255)
                                                                                                      cGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAAt                                        >  2:1519947/1‑148 (MQ=255)
                                                                                                      cGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAAt                                        <  1:385876/148‑1 (MQ=255)
                                                                                                        tGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATgg                                      <  2:3763842/148‑1 (MQ=255)
                                                                                                          cccGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCt                                    <  2:820670/148‑1 (MQ=255)
                                                                                                               cgggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTa                               <  1:3636223/148‑1 (MQ=255)
                                                                                                                gggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTaa                              <  2:2804785/148‑1 (MQ=255)
                                                                                                                gggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTaa                              <  2:889628/148‑1 (MQ=255)
                                                                                                                 ggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCGAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAAt                             <  2:1060088/148‑1 (MQ=255)
                                                                                                                 ggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAAt                             <  1:2128401/148‑1 (MQ=255)
                                                                                                                  gcgcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATc                            <  1:2228114/148‑1 (MQ=255)
                                                                                                                  gcgcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATc                            <  1:2001680/148‑1 (MQ=255)
                                                                                                                     cTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctg                         <  1:658881/148‑1 (MQ=255)
                                                                                                                     cTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctg                         >  2:2383462/1‑148 (MQ=255)
                                                                                                                        ttCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGGACTCctgctg                      >  1:2281217/1‑148 (MQ=255)
                                                                                                                            cTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCgg                  >  2:876636/1‑148 (MQ=255)
                                                                                                                             tGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGa                 <  1:2859427/148‑1 (MQ=255)
                                                                                                                                      tATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGttaa        <  2:863553/148‑1 (MQ=255)
                                                                                                                                       aTCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGttaat       >  2:1969040/1‑148 (MQ=255)
                                                                                                                                         cGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCCTCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGttaatta     >  2:2439295/1‑148 (MQ=255)
                                                                                                                                         cGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGttaatta     >  2:517968/1‑148 (MQ=255)
                                                                                                                                          gTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGttaattaa    <  2:1664066/148‑1 (MQ=255)
                                                                                                                                            aaGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAAcc  <  2:1363285/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
AAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACC  >  CP009273/1348629‑1348916

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: