Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 961,339 C→A 10.2% L522I (CTC→ATC)  ycaI → ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273961,3390CA10.2% 796.8 / 46.4 225L522I (CTC→ATC) ycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (97/105);  new base A (13/10);  total (110/115)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 242 total aligned reads displayed.

GGCGAACGTTGCCAGCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGATTCAATATTGCAT  >  CP009273/961190‑961488
                                                                                                                                                     |                                                                                                                                                     
ggCGAACGTTGCCCGCAATGTGGGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAATAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATtc                                                                                                                                                       >  2:4272907/1‑150 (MQ=255)
ggCGAACGTTGCCAGCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATtc                                                                                                                                                       >  1:3795821/1‑150 (MQ=255)
  cGAACGTTGCCAGCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATtctc                                                                                                                                                     <  1:2160236/150‑1 (MQ=255)
  cGAACGTTGCCAGCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATtatc                                                                                                                                                     >  1:817575/1‑150 (MQ=255)
       gTTGCCAGCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGa                                                                                                                                                >  1:1200613/1‑150 (MQ=255)
       gTTGCCAGCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGa                                                                                                                                                >  1:1200685/1‑150 (MQ=255)
        ttGCCAGCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGac                                                                                                                                               <  2:4035494/150‑1 (MQ=255)
           ccAGCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGAcaca                                                                                                                                            <  2:944413/150‑1 (MQ=255)
           ccAGCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGAcaca                                                                                                                                            <  1:758875/150‑1 (MQ=255)
            cAGCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAg                                                                                                                                           >  1:3842300/1‑150 (MQ=255)
             aGCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAgg                                                                                                                                          >  1:3970628/1‑150 (MQ=255)
              gCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGa                                                                                                                                         <  2:3062725/150‑1 (MQ=255)
              gCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGa                                                                                                                                         <  1:2227283/150‑1 (MQ=255)
               cAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGAc                                                                                                                                        <  2:4393533/150‑1 (MQ=255)
                aaTGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGAct                                                                                                                                       <  1:3766962/150‑1 (MQ=255)
                 aTGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGActg                                                                                                                                      >  2:2360765/1‑150 (MQ=255)
                 aTGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGActg                                                                                                                                      <  2:2863417/150‑1 (MQ=255)
                   gtgtgtgGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGActggc                                                                                                                                    <  2:1914881/150‑1 (MQ=255)
                   gtgtgtgGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGActggc                                                                                                                                    >  1:2296865/1‑150 (MQ=255)
                     gtgtgGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGActggcct                                                                                                                                  <  1:4384591/150‑1 (MQ=255)
                     gtgtgGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGActggcct                                                                                                                                  <  1:2199901/150‑1 (MQ=255)
                     gtgtgGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGActggcct                                                                                                                                  <  2:3656985/150‑1 (MQ=255)
                       gtgGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGActggcctgg                                                                                                                                <  1:1168064/150‑1 (MQ=255)
                        tgGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGActggcctggc                                                                                                                               >  1:1077156/1‑150 (MQ=255)
                        tgGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGActggcctggc                                                                                                                               >  1:1693487/1‑150 (MQ=255)
                        tgGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTATCTATGACACAGGActggcctggc                                                                                                                               >  2:2040710/1‑150 (MQ=255)
                            tGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGa                                                                                                                           <  2:2290344/150‑1 (MQ=255)
                            tGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGa                                                                                                                           <  1:4404715/150‑1 (MQ=255)
                              gAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAg                                                                                                                         <  1:316679/150‑1 (MQ=255)
                               aGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAgg                                                                                                                        <  2:2800935/150‑1 (MQ=255)
                               aGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAgg                                                                                                                        <  2:3167015/150‑1 (MQ=255)
                                 gCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAgggg                                                                                                                      >  1:541640/1‑150 (MQ=255)
                                       tGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTg                                                                                                                >  1:235408/1‑150 (MQ=255)
                                       tGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTg                                                                                                                >  2:1079706/1‑150 (MQ=255)
                                       tGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTg                                                                                                                <  1:1923755/150‑1 (MQ=255)
                                       tGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTATCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTg                                                                                                                <  2:4240287/150‑1 (MQ=255)
                                           gtgCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGca                                                                                                            <  1:2095766/150‑1 (MQ=255)
                                              cTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGcaaca                                                                                                         >  2:322377/1‑150 (MQ=255)
                                              cTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGcaaca                                                                                                         <  2:1737183/150‑1 (MQ=255)
                                               tGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTATCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGcaacaa                                                                                                        >  2:3348554/1‑150 (MQ=255)
                                                 gCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACt                                                                                                      <  1:2790400/150‑1 (MQ=255)
                                                 gCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACt                                                                                                      <  2:3715413/150‑1 (MQ=255)
                                                   cGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGa                                                                                                    >  2:4031082/1‑150 (MQ=255)
                                                            aaaaaCCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCt                                                                                           <  2:685572/150‑1 (MQ=255)
                                                            aaaaaCCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCt                                                                                           >  1:4502802/1‑150 (MQ=255)
                                                            aaaaaCCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTATCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCt                                                                                           >  1:4125038/1‑150 (MQ=255)
                                                            aaaaaCCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTATCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCt                                                                                           >  1:4125093/1‑150 (MQ=255)
                                                             aaaaCCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTg                                                                                          >  1:714880/1‑150 (MQ=255)
                                                             aaaaCCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTAGCCCCTg                                                                                          >  1:3626717/1‑150 (MQ=255)
                                                              aaaCCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTgg                                                                                         >  2:4445078/1‑150 (MQ=255)
                                                                  cTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTcc                                                                                     <  1:4262943/150‑1 (MQ=255)
                                                                  cTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTcc                                                                                     <  1:2160569/150‑1 (MQ=255)
                                                                   tCGACCTGACGAGTGGCAGCTGTACATGCTTTATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCa                                                                                    >  2:3385741/1‑150 (MQ=255)
                                                                   tCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCa                                                                                    >  2:3329875/1‑150 (MQ=255)
                                                                   tCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTATCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCa                                                                                    >  1:2534045/1‑150 (MQ=255)
                                                                   tCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTATCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCa                                                                                    >  1:234322/1‑150 (MQ=255)
                                                                          gACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCAt                                                                             >  2:4134455/1‑150 (MQ=255)
                                                                           aCGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATa                                                                            <  1:3969491/150‑1 (MQ=255)
                                                                            cGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGATACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATaa                                                                           <  1:2971749/150‑1 (MQ=255)
                                                                            cGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATaa                                                                           >  1:2866647/1‑150 (MQ=255)
                                                                            agAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGTCAACAACTGATTATCCCCTGGCTCCACTGGCATaa                                                                           <  1:2911293/149‑1 (MQ=255)
                                                                              aGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATc                                                                         <  1:3384604/150‑1 (MQ=255)
                                                                              aGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATc                                                                         >  1:2842577/1‑150 (MQ=255)
                                                                              aGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATc                                                                         >  1:2393760/1‑150 (MQ=255)
                                                                              aGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTATCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATc                                                                         >  1:1747780/1‑150 (MQ=255)
                                                                                tGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAATCGGCAAAGCGTTTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCtt                                                                       <  2:1367699/150‑1 (MQ=255)
                                                                                tGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCtt                                                                       <  2:622314/150‑1 (MQ=255)
                                                                                 ggCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTg                                                                      >  1:1000120/1‑150 (MQ=255)
                                                                                 ggCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTg                                                                      <  1:1439001/150‑1 (MQ=255)
                                                                                    aGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAAc                                                                   >  1:864771/1‑150 (MQ=255)
                                                                                     gCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTATCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAAcc                                                                  <  1:1827083/150‑1 (MQ=255)
                                                                                       tGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCgg                                                                >  2:2996328/1‑150 (MQ=255)
                                                                                         tACATGCTTGATGTCGGGCAAGGGCTGGCAATGGAGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGaa                                                              >  2:4119872/1‑150 (MQ=255)
                                                                                          aCATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAg                                                             <  1:3588537/150‑1 (MQ=255)
                                                                                                 tGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTAt                                                      >  2:508787/1‑150 (MQ=255)
                                                                                                  gATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTAtt                                                     <  1:3114509/150‑1 (MQ=255)
                                                                                                       cGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAg                                                >  2:1194660/1‑150 (MQ=255)
                                                                                                          gCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCa                                             <  1:2819001/150‑1 (MQ=255)
                                                                                                          gCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCa                                             <  2:3842300/150‑1 (MQ=255)
                                                                                                            aaGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATg                                           >  2:1304514/1‑150 (MQ=255)
                                                                                                            aaGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATg                                           >  1:3243282/1‑150 (MQ=255)
                                                                                                            aaGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATg                                           >  1:3357695/1‑150 (MQ=255)
                                                                                                             aGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGa                                          <  1:2182528/150‑1 (MQ=255)
                                                                                                                  tGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATc                                     <  2:2113176/150‑1 (MQ=255)
                                                                                                                  tGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATc                                     <  2:2112969/150‑1 (MQ=255)
                                                                                                                  tGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATc                                     >  1:352743/1‑150 (MQ=255)
                                                                                                                    gCAATGGTGATAGCCAGAAACGGCAATGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTg                                   >  2:1440600/1‑150 (MQ=255)
                                                                                                                    gCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTg                                   <  1:1516878/150‑1 (MQ=255)
                                                                                                                    gCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTg                                   <  2:4042610/150‑1 (MQ=255)
                                                                                                                       aTGGTGATAGCCAGAAACGGCAAAGCGATTATCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGAt                                <  1:2624644/150‑1 (MQ=255)
                                                                                                                               aGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCgggg                        >  2:140432/1‑150 (MQ=255)
                                                                                                                               aGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCACCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAAAATCTGGATCAACgggg                        >  2:1158930/1‑150 (MQ=255)
                                                                                                                                    gAAACGGAAAAGCTATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATTATCTTGAACCGGAAGGCGTTTTTCTGAGCCATGAACATCTGGATCACCGGGGAGGGc                   <  2:2644805/150‑1 (MQ=255)
                                                                                                                                      aaCGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTg                 >  1:4045506/1‑150 (MQ=255)
                                                                                                                                       aCGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTgg                <  1:2023354/150‑1 (MQ=255)
                                                                                                                                       aCGGCAAAGCGATTATCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTgg                <  1:3852979/150‑1 (MQ=255)
                                                                                                                                         ggCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGAt              >  1:999342/1‑150 (MQ=255)
                                                                                                                                                cGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGATTCAatat       >  2:3377104/1‑150 (MQ=255)
                                                                                                                                                 gATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGATTCAATAtt      <  1:1333471/150‑1 (MQ=255)
                                                                                                                                                     ctctATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGACCATCTGGATCACCGGGGAGGGCTGGATTCAATATTGCat  <  2:4480824/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGAACGTTGCCAGCAATGTGTGTGGCTGGAGGCTTGCTGATGTGCTGGCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGATTCAATATTGCAT  >  CP009273/961190‑961488

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