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% 759.1 / ‑2.1 236P2T (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/3);  new base T (96/137);  total (96/140)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 236 total aligned reads displayed.

GCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACC  >  CP009273/1348627‑1348916
                                                                                                                                                 |                                                                                                                                                
gcAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCa                                                                                                                                                <  1:372238/148‑1 (MQ=255)
 cAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCAt                                                                                                                                               <  2:3440418/148‑1 (MQ=255)
   aaCCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATAc                                                                                                                                             <  2:3017872/148‑1 (MQ=255)
      cAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCAt                                                                                                                                          <  2:2869444/148‑1 (MQ=255)
       aGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCAtt                                                                                                                                         <  1:13548/148‑1 (MQ=255)
       aGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCAtt                                                                                                                                         >  2:480050/1‑148 (MQ=255)
       aGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCAtt                                                                                                                                         <  2:396811/148‑1 (MQ=255)
        gccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTc                                                                                                                                        <  1:1423131/148‑1 (MQ=255)
        gccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTc                                                                                                                                        <  1:3144911/148‑1 (MQ=255)
        gccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTc                                                                                                                                        >  2:3063759/1‑148 (MQ=255)
         ccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTcc                                                                                                                                       >  2:1893435/1‑148 (MQ=255)
         ccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTcc                                                                                                                                       >  1:2801210/1‑148 (MQ=255)
         ccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTcc                                                                                                                                       <  1:1512200/148‑1 (MQ=255)
         ccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTcc                                                                                                                                       >  1:1878890/1‑148 (MQ=255)
         ccgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTcc                                                                                                                                       <  1:3203081/148‑1 (MQ=255)
          cgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCt                                                                                                                                      <  2:1043518/148‑1 (MQ=255)
           gccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACAATTCCtt                                                                                                                                     <  2:2705651/148‑1 (MQ=255)
            ccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCCGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTa                                                                                                                                    >  2:325691/1‑148 (MQ=255)
            ccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTa                                                                                                                                    <  1:2799885/148‑1 (MQ=255)
            ccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTa                                                                                                                                    >  2:2094921/1‑148 (MQ=255)
            ccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTa                                                                                                                                    <  2:1302639/148‑1 (MQ=255)
             cAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTAt                                                                                                                                   >  1:3217118/1‑148 (MQ=255)
              aaaaaTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATg                                                                                                                                  <  2:617302/148‑1 (MQ=255)
                aaaTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGtt                                                                                                                                <  2:2914688/148‑1 (MQ=255)
                 aaTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGttt                                                                                                                               <  1:3176564/148‑1 (MQ=255)
                  aTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTc                                                                                                                              <  1:1058462/148‑1 (MQ=255)
                     cacaGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCCGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCaga                                                                                                                           >  2:1945821/1‑148 (MQ=255)
                          tACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGtt                                                                                                                      >  2:1633015/1‑148 (MQ=255)
                          tACAGCCAGCCCCGCGCAGCCGTACAGGCAGACCATCGCTGAGTCGTAACTATAAAATTTGCGCCAGGCGGTACTCAGCGTGCCCGGCGTGCGCTTTTCGCTGTATACTCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGtt                                                                                                                      <  2:687554/148‑1 (MQ=255)
                           aCAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTa                                                                                                                     <  2:779877/148‑1 (MQ=255)
                                 agccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATg                                                                                                               <  2:1819227/148‑1 (MQ=255)
                                  gccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATgg                                                                                                              >  2:2852310/1‑148 (MQ=255)
                                  gccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCGATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGGTAGCCCTgg                                                                                                              >  1:358459/1‑148 (MQ=255)
                                   ccccGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGc                                                                                                             <  1:1133632/148‑1 (MQ=255)
                                      cgcgCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGttt                                                                                                                               >  2:2472177/1‑127 (MQ=255)
                                      cgcgCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGttt                                                                                                                               <  1:2472177/127‑1 (MQ=255)
                                      cgcgCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAcc                                                                                                          <  2:2221295/148‑1 (MQ=255)
                                          cAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACACaaa                                                                                                      >  2:1005961/1‑148 (MQ=255)
                                            gCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCCTACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAat                                                                                                    >  1:1915093/1‑148 (MQ=255)
                                                 aCAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAg                                                                                               <  2:3473352/148‑1 (MQ=255)
                                                    ggCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGaaa                                                                                            >  2:875658/1‑148 (MQ=255)
                                                     gCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAAt                                                                                           <  2:1610994/148‑1 (MQ=255)
                                                      ccGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATc                                                                                          >  1:2631405/1‑148 (MQ=255)
                                                       cGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCa                                                                                         <  2:87623/148‑1 (MQ=255)
                                                          ccATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTGGCTGTATAAGCTATCGTAAGTCATAACAATCCTTATGTTTCAGAGGGTTAGCCAAGGCACCCAAAAAATAAGAAAAAACGt                                                                                      >  1:45586/1‑148 (MQ=255)
                                                          ccATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTGGCTGTATAAGCTATAGTAATTCATACCAATCCTTATTTTTCAGAGGGTTAGCCAAGGCACCCAAAATATAAGAAAAAACGt                                                                                      >  1:45638/1‑148 (MQ=255)
                                                          ccATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGt                                                                                      >  1:810206/1‑148 (MQ=255)
                                                            aTCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                                                    <  1:2497128/148‑1 (MQ=255)
                                                             tCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTaa                                                                                   <  2:2652785/148‑1 (MQ=255)
                                                              cGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTGGCTGTATACGCGATCGTAAGTCATAACAATACGTATGTTTCAGATGGTTAGCCATGGCCCCCAAAATAAAAGAAATAACGGTAAc                                                                                  >  1:935257/1‑148 (MQ=255)
                                                              cGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAAc                                                                                  >  2:3286343/1‑148 (MQ=255)
                                                               gCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACa                                                                                 <  1:760071/148‑1 (MQ=255)
                                                                  gagGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaata                                                                              >  1:2342499/1‑148 (MQ=255)
                                                                    ggCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaataat                                                                            <  1:1087369/148‑1 (MQ=255)
                                                                       gTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGAc                                                                         <  1:1811751/148‑1 (MQ=255)
                                                                        tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCCTTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGAcc                                                                        >  2:3294390/1‑148 (MQ=255)
                                                                        tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGAcc                                                                        <  1:1060667/148‑1 (MQ=255)
                                                                        tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGAcc                                                                        <  2:2276648/148‑1 (MQ=255)
                                                                        tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGAcc                                                                        <  1:325945/148‑1 (MQ=255)
                                                                            tataAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTg                                                                    <  2:1927600/148‑1 (MQ=255)
                                                                            tataAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTg                                                                    >  2:1519730/1‑148 (MQ=255)
                                                                            tataAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTg                                                                    <  2:1901605/148‑1 (MQ=255)
                                                                            tataAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGATTTTCGATGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAAAAATTACCAGTg                                                                    >  2:112463/1‑148 (MQ=255)
                                                                               aaaaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGc                                                                 >  2:508981/1‑148 (MQ=255)
                                                                                aaaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGcc                                                                <  2:1863582/148‑1 (MQ=255)
                                                                                  aTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCcac                                                              >  1:1450148/1‑148 (MQ=255)
                                                                                    ttGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTGGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCcacac                                                            >  1:978172/1‑148 (MQ=255)
                                                                                    ttGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCcacac                                                            >  1:2890886/1‑148 (MQ=255)
                                                                                         ccAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCAt                                                       <  2:3171453/148‑1 (MQ=255)
                                                                                         ccAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCAt                                                       <  1:3297687/148‑1 (MQ=255)
                                                                                          cAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATc                                                      <  2:439992/148‑1 (MQ=255)
                                                                                           aGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCCAAATATCAGAAATCACGTTAACAAAAATTACCAGTGAGCCACACACCATCa                                                     >  2:1197090/1‑148 (MQ=255)
                                                                                            ggCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCAc                                                    >  2:427306/1‑148 (MQ=255)
                                                                                              cGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTc                                                  <  1:2623224/148‑1 (MQ=255)
                                                                                                gTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcg                                                <  1:2057873/148‑1 (MQ=255)
                                                                                                     cAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGa                                           >  1:3440457/1‑148 (MQ=255)
                                                                                                      aGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaa                                          <  1:855225/148‑1 (MQ=255)
                                                                                                      aGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaa                                          <  1:2793800/148‑1 (MQ=255)
                                                                                                       gCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaaa                                         <  1:2452626/148‑1 (MQ=255)
                                                                                                       gCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaaa                                         <  1:684062/148‑1 (MQ=255)
                                                                                                           gCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGc                                     <  2:3108179/148‑1 (MQ=255)
                                                                                                                  gggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTaa                              <  1:2751023/148‑1 (MQ=255)
                                                                                                                  gggcgCTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTaa                              <  1:997483/148‑1 (MQ=255)
                                                                                                                   ggcgCTTTTCGCTGTATACGCTATCGTAAGTCATAACATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCAACACTCCGGCGGAAATGGCTGCGTAAt                             >  1:487856/1‑148 (MQ=255)
                                                                                                                    gcgcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATc                            >  2:1389829/1‑148 (MQ=255)
                                                                                                                    gcgcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATc                            >  2:3292785/1‑148 (MQ=255)
                                                                                                                     cgcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATcc                           <  1:2248609/148‑1 (MQ=255)
                                                                                                                     cgcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATcc                           >  1:650012/1‑148 (MQ=255)
                                                                                                                     cgcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATcc                           <  1:2641370/148‑1 (MQ=255)
                                                                                                                     cgcTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATcc                           <  2:45586/148‑1 (MQ=255)
                                                                                                                       cTTTTCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctg                         >  2:109005/1‑148 (MQ=255)
                                                                                                                        ttttCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgc                        >  1:1582959/1‑148 (MQ=255)
                                                                                                                         tttCGCGGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgct                       <  2:2620167/148‑1 (MQ=255)
                                                                                                                          ttCGCTGTATACGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgctg                      <  1:1040378/148‑1 (MQ=255)
                                                                                                                                   taCGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGc             <  1:875669/148‑1 (MQ=255)
                                                                                                                                     cGCTATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCCCGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGt           >  2:1041651/1‑148 (MQ=255)
                                                                                                                                        tATCGTAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGttaa        <  2:2041934/148‑1 (MQ=255)
                                                                                                                                             tAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACGCCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAAc   >  1:1799045/1‑148 (MQ=255)
                                                                                                                                             tAAGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAAc   >  2:11756/1‑148 (MQ=255)
                                                                                                                                              aaGTCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAAcc  >  2:2033431/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACC  >  CP009273/1348627‑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: