Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,407,969 C→A 94.8% N249K (AAC→AAA acrE → cytoplasmic membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,407,9690CA94.8% 659.5 / 26.6 233N249K (AAC→AAAacrEcytoplasmic membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (5/7);  new base A (129/92);  total (134/99)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

Only 100 of 233 total aligned reads displayed.

AACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATTGA  >  CP009273/3407824‑3408112
                                                                                                                                                 |                                                                                                                                               
aaCTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAgg                                                                                                                                               >  1:2784944/1‑148 (MQ=255)
aaCTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAgg                                                                                                                                               >  1:2412295/1‑148 (MQ=255)
 aCTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCCATCCAGCAACGATTTTATGAGTCTGAAGCAATCCGTATAGCAAGGAAATTTGCATAAGGACAACGGCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGt                                                                                                                                              <  2:258718/148‑1 (MQ=255)
    ggCGCCTGTCCAGCAGCTCGATCCTATCTACTTTGATGTGACCAAATCCAGCAACTATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGTAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCaa                                                                                                                                           <  1:1561365/148‑1 (MQ=255)
    ggCGAGTGTCCAGCAGCTCTATCCTATCTACTTTTATGTCACCCAATCCAGAAACTATTTTATGTGGCTTAAGCAATCCTTATAGAAATTAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCaa                                                                                                                                           <  1:1504068/148‑1 (MQ=255)
     gCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCACCGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCaaa                                                                                                                                          <  1:2252984/148‑1 (MQ=255)
     gCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCaaa                                                                                                                                          <  2:576389/148‑1 (MQ=255)
      cGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTCACCCAAACCCGCAACGATTTTATGCGGCTGAAGCAATCCGTAGAGCAAGCAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAAc                                                                                                                                         <  1:638135/148‑1 (MQ=255)
         cTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTa                                                                                                                                      <  1:1926920/148‑1 (MQ=255)
          tGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTAt                                                                                                                                     <  2:2465396/148‑1 (MQ=255)
           gTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATc                                                                                                                                    >  1:2055548/1‑148 (MQ=255)
                   gcTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGaaa                                                                                                                            <  2:1061372/148‑1 (MQ=255)
                    cTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAg                                                                                                                           >  1:79072/1‑148 (MQ=255)
                     tCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAgg                                                                                                                          >  2:2174724/1‑148 (MQ=255)
                     tCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAgg                                                                                                                          >  1:2237958/1‑148 (MQ=255)
                      cGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGt                                                                                                                         <  1:2338962/148‑1 (MQ=255)
                        aTCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAAATTGCATAAGGAAAACGCCACCAGCAACGTAGAGTAGGTCATGGAAAACGGTCAAACCAATCCCCAGAAAGGTAc                                                                                                                       >  1:2324777/1‑148 (MQ=255)
                         tCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACg                                                                                                                      >  1:2583335/1‑148 (MQ=255)
                          ccTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGACAAAGGTCAAACCTCTCCCCTGAAAGGTACGc                                                                                                                     >  2:3484704/1‑148 (MQ=255)
                              tCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCa                                                                                                                 >  1:2194369/1‑148 (MQ=255)
                               cTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCaa                                                                                                                >  2:3054161/1‑148 (MQ=255)
                                  cGTTGATGTGACCCAATCCACCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAACACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATtc                                                                                                             <  1:1647964/148‑1 (MQ=255)
                                    ttGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATtctc                                                                                                           >  1:183979/1‑148 (MQ=255)
                                    ttGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATtctc                                                                                                           <  2:2371072/148‑1 (MQ=255)
                                    ttGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATtctc                                                                                                           <  2:2662071/148‑1 (MQ=255)
                                         gtgACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtg                                                                                                      <  2:324228/148‑1 (MQ=255)
                                             cccACTCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGAc                                                                                                  <  2:2924875/148‑1 (MQ=255)
                                             cccAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGAc                                                                                                  <  1:2433698/148‑1 (MQ=255)
                                                  tCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTtg                                                                                             >  2:2470050/1‑148 (MQ=255)
                                                     aGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTtgatg                                                                                          <  1:68726/148‑1 (MQ=255)
                                                     aGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTtgatg                                                                                          <  2:897261/148‑1 (MQ=255)
                                                       cAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGaa                                                                                        >  1:1819516/1‑148 (MQ=255)
                                                        aaCGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGaaa                                                                                       >  1:2550931/1‑148 (MQ=255)
                                                         aCGATTTTATGGGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAg                                                                                      >  1:2020373/1‑148 (MQ=255)
                                                         aCGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAg                                                                                      >  1:3060784/1‑148 (MQ=255)
                                                          cGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGc                                                                                     >  1:1638890/1‑148 (MQ=255)
                                                          cGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGc                                                                                     >  2:1177517/1‑148 (MQ=255)
                                                           gATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCa                                                                                    <  1:3432899/148‑1 (MQ=255)
                                                            aTTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCAc                                                                                   <  2:1105398/148‑1 (MQ=255)
                                                                  tGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTc                                                                             >  1:1001441/1‑148 (MQ=255)
                                                                   gAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTcc                                                                            >  1:3143421/1‑148 (MQ=255)
                                                                    aGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCa                                                                           >  1:421109/1‑148 (MQ=255)
                                                                     ggCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCAt                                                                          >  2:3428860/1‑148 (MQ=255)
                                                                      gCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATa                                                                         >  2:1740324/1‑148 (MQ=255)
                                                                          aaGCAATCCTTAGAGCAAGGAAATTTGAATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAAccc                                                                     <  1:1547716/148‑1 (MQ=255)
                                                                          aaGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAAccc                                                                     >  2:866437/1‑148 (MQ=255)
                                                                           aGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCt                                                                    >  2:1377036/1‑148 (MQ=255)
                                                                              aaTCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACg                                                                 <  2:672205/148‑1 (MQ=255)
                                                                               aTCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGt                                                                >  2:760264/1‑148 (MQ=255)
                                                                                tCCGTAGAGCAAGGAAATTTGAATAAGCAACACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTg                                                               <  1:2428577/148‑1 (MQ=255)
                                                                                 ccGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGc                                                              <  2:1464619/148‑1 (MQ=255)
                                                                                 ccGTAGAGAAAGTAAATTTGCATAAGTAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGc                                                              <  1:807180/148‑1 (MQ=255)
                                                                                    tAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGt                                                           >  2:1658140/1‑148 (MQ=255)
                                                                                     agagCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTc                                                          >  1:754112/1‑148 (MQ=255)
                                                                                      gagCAAGGAAATTTGCATAAGTAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGCTCAAACCTATCCCCTGAAAGGTACGCTGCAAATCTCCGATGTGACCGTAGATGAACGCACCGGCACAATAACCCTCCCTGctctct                                                         >  1:3019897/1‑148 (MQ=255)
                                                                                       agCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCtt                                                        <  2:1879988/148‑1 (MQ=255)
                                                                                        gaaaGGAAATTTGCATAAGCAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTc                                                       <  1:1274219/146‑1 (MQ=255)
                                                                                        gCACGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTc                                                       <  1:1409488/148‑1 (MQ=255)
                                                                                         cAAGGAAATTTGCATAATGAAAACGCCACCAGCATAGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTTAAAGGTACGATGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTcc                                                      <  1:870824/148‑1 (MQ=255)
                                                                                         cAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTcc                                                      <  2:709211/148‑1 (MQ=255)
                                                                                         cAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCAGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTca                                                      >  2:2122020/1‑147 (MQ=255)
                                                                                          aaGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTccc                                                     <  2:2412295/148‑1 (MQ=255)
                                                                                           aGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCt                                                    >  1:2841949/1‑148 (MQ=255)
                                                                                           aGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCt                                                    >  2:1044076/1‑148 (MQ=255)
                                                                                           aGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCt                                                    >  1:2219397/1‑148 (MQ=255)
                                                                                            ggAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTa                                                   >  2:1828520/1‑148 (MQ=255)
                                                                                                aTTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCTGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAAccc                                               >  1:1067446/1‑148 (MQ=255)
                                                                                                     cATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAAc                                          <  1:284193/148‑1 (MQ=255)
                                                                                                     cATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAAc                                          >  1:2742162/1‑148 (MQ=255)
                                                                                                      aTAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACa                                         >  1:447218/1‑148 (MQ=255)
                                                                                                      aTAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACa                                         >  2:1911688/1‑148 (MQ=255)
                                                                                                      aTAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACa                                         <  2:2053429/148‑1 (MQ=255)
                                                                                                       tAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACAt                                        <  1:563639/148‑1 (MQ=255)
                                                                                                        aaGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATa                                       <  2:2841949/148‑1 (MQ=255)
                                                                                                        aaGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATa                                       >  2:1463796/1‑148 (MQ=255)
                                                                                                        aaGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATa                                       <  2:2564069/148‑1 (MQ=255)
                                                                                                          ggAAAACGCCACCAGCAACGTAGAGTTGGTCATGAAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACg                                     >  1:2954126/1‑148 (MQ=255)
                                                                                                             aaaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCtt                                  <  2:2788606/148‑1 (MQ=255)
                                                                                                              aaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCttt                                 >  2:1866662/1‑148 (MQ=255)
                                                                                                               aCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCtttt                                <  2:609233/148‑1 (MQ=255)
                                                                                                                  ccaccaGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGcc                             >  1:1007026/1‑148 (MQ=255)
                                                                                                                  ccaccaGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGcc                             <  1:1123713/148‑1 (MQ=255)
                                                                                                                     ccaGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCggg                          >  2:3328019/1‑148 (MQ=255)
                                                                                                                       aGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTa                        >  2:834444/1‑148 (MQ=255)
                                                                                                                        gCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTAt                       >  1:3501279/1‑148 (MQ=255)
                                                                                                                        gCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCAGAAAGGTAAGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCATAACACGCAACATACGCTTTTGCCGGGtct                       >  1:3500330/1‑146 (MQ=255)
                                                                                                                         cAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATg                      >  1:2125154/1‑148 (MQ=255)
                                                                                                                         cAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATg                      <  1:2174724/148‑1 (MQ=255)
                                                                                                                           aCGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGtt                    >  1:912392/1‑148 (MQ=255)
                                                                                                                           aCGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGtt                    >  1:2689369/1‑148 (MQ=255)
                                                                                                                             gTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTtg                  <  2:936464/148‑1 (MQ=255)
                                                                                                                              tAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTtgt                 >  1:2326686/1‑148 (MQ=255)
                                                                                                                              tAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTtgt                 >  2:2279100/1‑148 (MQ=255)
                                                                                                                                    tGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTgtgcgtgc           >  2:1497235/1‑148 (MQ=255)
                                                                                                                                        cATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACgg       <  1:123597/148‑1 (MQ=255)
                                                                                                                                        cATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACgg       <  1:1175347/148‑1 (MQ=255)
                                                                                                                                         aTGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGa      >  1:1035704/1‑148 (MQ=255)
                                                                                                                                          tGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGAt     >  1:2759901/1‑148 (MQ=255)
                                                                                                                                           ggAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGAtt    >  2:2961863/1‑148 (MQ=255)
                                                                                                                                             aaaaAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCACCATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATtga  >  2:2347258/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATTGA  >  CP009273/3407824‑3408112

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: