Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,209,162 C→A 71.3% I9I (ATC→ATA yohO → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,209,1620CA71.3% 374.5 / 282.9 289I9I (ATC→ATAyohOputative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (42/41);  new base A (119/87);  total (161/128)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 289 total aligned reads displayed.

CAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAAC  >  CP009273/2209017‑2209305
                                                                                                                                                 |                                                                                                                                               
cagcagCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGccc                                                                                                                                             >  1:4282515/1‑150 (MQ=255)
   cagcaACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCcctgt                                                                                                                                          >  2:3426012/1‑150 (MQ=255)
    agcaACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAATTTTGGGTCATAGCcctgtt                                                                                                                                         >  1:412329/1‑150 (MQ=255)
    agcaACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCcctgtt                                                                                                                                         >  1:116600/1‑150 (MQ=255)
       aaCTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCcctgttcct                                                                                                                                      >  2:2202335/1‑150 (MQ=255)
          tCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCcctgttcctgtt                                                                                                                                   >  2:681963/1‑150 (MQ=255)
           cccccggtGGCAAGAACACGGCCAGGACATCTACCCACTCGGCATTTACTGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGttt                                                                                                                                  <  1:4295291/143‑1 (MQ=255)
           cACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGttt                                                                                                                                  <  1:3462148/150‑1 (MQ=255)
             ccAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTAt                                                                                                                                <  2:2101635/150‑1 (MQ=255)
               aGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATgg                                                                                                                              <  2:4018716/150‑1 (MQ=255)
                gTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGc                                                                                                                             <  1:107240/150‑1 (MQ=255)
                   ggCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGtt                                                                                                                          <  1:4471134/150‑1 (MQ=255)
                   ggCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGtt                                                                                                                          >  1:3905118/1‑150 (MQ=255)
                      aaGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAgg                                                                                                                       >  1:1131964/1‑150 (MQ=255)
                      aaGAACAAGGCCATGATATTTTCCCACTCGGCATTAAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAgg                                                                                                                       <  2:2977277/150‑1 (MQ=255)
                         aaCAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCgg                                                                                                                    >  2:2513043/1‑150 (MQ=255)
                          aCAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGa                                                                                                                   <  1:3399265/150‑1 (MQ=255)
                          aCAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGa                                                                                                                   <  1:3258569/150‑1 (MQ=255)
                            aaGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAAt                                                                                                                 >  1:3256682/1‑150 (MQ=255)
                                  aTGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtggtgg                                                                                                           <  1:760423/150‑1 (MQ=255)
                                   tGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGc                                                                                                          <  2:407765/150‑1 (MQ=255)
                                   tGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGc                                                                                                          <  2:1294410/150‑1 (MQ=255)
                                   tGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGc                                                                                                          <  2:407841/150‑1 (MQ=255)
                                    gATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCg                                                                                                         <  2:1959493/150‑1 (MQ=255)
                                     atatTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGt                                                                                                        >  2:1148975/1‑150 (MQ=255)
                                      tatTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTc                                                                                                       >  1:1396101/1‑150 (MQ=255)
                                      tatTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTc                                                                                                       >  1:65727/1‑150 (MQ=255)
                                      tatTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGACATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTc                                                                                                       >  1:3664403/1‑150 (MQ=255)
                                       atTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCa                                                                                                      <  2:31207/150‑1 (MQ=255)
                                       atTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCa                                                                                                      >  1:3990556/1‑150 (MQ=255)
                                         tttCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATg                                                                                                    <  1:356205/150‑1 (MQ=255)
                                         tttCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATg                                                                                                    >  1:2307237/1‑150 (MQ=255)
                                          ttCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGc                                                                                                   >  1:1888010/1‑150 (MQ=255)
                                          ttCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGc                                                                                                   >  2:1535191/1‑150 (MQ=255)
                                           tCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCt                                                                                                  >  1:2572463/1‑150 (MQ=255)
                                           tCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCt                                                                                                  <  2:94591/150‑1 (MQ=255)
                                           tCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCt                                                                                                  <  2:431820/150‑1 (MQ=255)
                                            cccACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTc                                                                                                 <  2:1631313/150‑1 (MQ=255)
                                            cccACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTc                                                                                                 >  2:1057773/1‑150 (MQ=255)
                                             ccACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCg                                                                                                <  2:3732410/150‑1 (MQ=255)
                                             ccACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCg                                                                                                <  2:555287/150‑1 (MQ=255)
                                                cTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAg                                                                                             <  2:3275399/150‑1 (MQ=255)
                                                cTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAg                                                                                             >  2:3726802/1‑150 (MQ=255)
                                                cTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCAGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAg                                                                                             >  1:3720910/1‑150 (MQ=255)
                                                 tCGGCATTAAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAgg                                                                                            <  2:89809/150‑1 (MQ=255)
                                                  cGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGt                                                                                           >  1:3590487/1‑150 (MQ=255)
                                                      aTTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTtata                                                                                       >  1:2136043/1‑150 (MQ=255)
                                                         tAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCt                                                                                    >  1:4388766/1‑150 (MQ=255)
                                                          aaTGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGCGGCGTCATGCTCGCAGGTTATACCtt                                                                                   <  1:3944720/150‑1 (MQ=255)
                                                                 cGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTg                                                                            >  2:2682946/1‑150 (MQ=255)
                                                                  gTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTATCATAAGGAGGGCATTAGGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGGCATGCTCGCAGGGTATACCTTTATTTTGc                                                                           >  2:705412/1‑150 (MQ=255)
                                                                  gTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGc                                                                           >  2:3886061/1‑150 (MQ=255)
                                                                   ttAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCg                                                                          <  1:3988672/150‑1 (MQ=255)
                                                                   ttAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCGGGTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCGGGTTATACCTTTATTTTGCg                                                                          >  2:3117309/1‑150 (MQ=255)
                                                                    tAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGt                                                                         >  2:3679563/1‑150 (MQ=255)
                                                                       cAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCt                                                                      <  1:146010/150‑1 (MQ=255)
                                                                          aTGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGc                                                                   >  2:1217525/1‑150 (MQ=255)
                                                                           tGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCt                                                                  >  1:864543/1‑150 (MQ=255)
                                                                           tGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCt                                                                  >  1:3142629/1‑150 (MQ=255)
                                                                            gcgcTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTa                                                                 <  1:2595210/150‑1 (MQ=255)
                                                                             cgcTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTaa                                                                >  2:1611070/1‑150 (MQ=255)
                                                                               cTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAgc                                                              >  1:1689504/1‑150 (MQ=255)
                                                                                 ttACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAgcgc                                                            >  1:465091/1‑150 (MQ=255)
                                                                                  tACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGcc                                                           <  1:4136172/150‑1 (MQ=255)
                                                                                    cTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTg                                                         >  2:4444512/1‑150 (MQ=255)
                                                                                    cTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTg                                                         >  1:643320/1‑150 (MQ=255)
                                                                                        aCAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCAcc                                                     <  1:998043/150‑1 (MQ=255)
                                                                                             ccTTTCGTTAACATAAGGAGTGCATTAAGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCt                                                <  2:3039892/150‑1 (MQ=255)
                                                                                                ttCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTc                                             >  2:3061476/1‑150 (MQ=255)
                                                                                                  cGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTCTTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCaa                                           >  2:1035143/1‑150 (MQ=255)
                                                                                                   gTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCaaa                                          <  1:3588014/150‑1 (MQ=255)
                                                                                                         aTAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCg                                    >  2:455158/1‑150 (MQ=255)
                                                                                                         aTAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCg                                    >  2:455076/1‑150 (MQ=255)
                                                                                                          tAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCgg                                   >  2:846820/1‑150 (MQ=255)
                                                                                                           aaGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGt                                  >  2:271281/1‑150 (MQ=255)
                                                                                                                gTGCATTATGCGCATAGCTACAATTTGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCg                             <  1:1779965/150‑1 (MQ=255)
                                                                                                                gTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCg                             <  1:277381/150‑1 (MQ=255)
                                                                                                                gTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCg                             >  2:2400706/1‑150 (MQ=255)
                                                                                                                 tGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCga                            <  2:4387456/150‑1 (MQ=255)
                                                                                                                    aTTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCgatga                         >  2:1418261/1‑150 (MQ=255)
                                                                                                                    aTTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCgatga                         >  2:1806723/1‑150 (MQ=255)
                                                                                                                    aTTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGACGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCgatga                         >  1:1024825/1‑150 (MQ=255)
                                                                                                                      tATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATc                       >  1:322305/1‑150 (MQ=255)
                                                                                                                       aTGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcg                      >  2:1105844/1‑150 (MQ=255)
                                                                                                                       aTGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcg                      >  2:2781430/1‑150 (MQ=255)
                                                                                                                        tGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGACTTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcgc                     >  1:4066750/1‑150 (MQ=255)
                                                                                                                         gcgcATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcgcc                    <  1:2052489/150‑1 (MQ=255)
                                                                                                                          cgcATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcgccg                   <  1:489337/150‑1 (MQ=255)
                                                                                                                          cgcATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcgccg                   >  1:2411439/1‑150 (MQ=255)
                                                                                                                            cATAGCTATAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGTATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTTAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGTTcgccgcg                 <  2:4374843/150‑2 (MQ=255)
                                                                                                                            cATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCgccgcc                 <  2:116600/150‑1 (MQ=255)
                                                                                                                             aTAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCa                <  1:2202335/150‑1 (MQ=255)
                                                                                                                             aTAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCa                <  1:2513043/150‑1 (MQ=255)
                                                                                                                               aGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGt              >  2:3837332/1‑150 (MQ=255)
                                                                                                                               aGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGACTTGGTGGCGTCATGCTCGCAGGTTATACCTGTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGt              >  1:3310133/1‑150 (MQ=255)
                                                                                                                               aGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGt              >  1:3310675/1‑150 (MQ=255)
                                                                                                                                   aacaTTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGcca          <  2:1034215/147‑1 (MQ=255)
                                                                                                                                       ttGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAg      <  2:350684/150‑1 (MQ=255)
                                                                                                                                         ggggTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTa    >  2:2451142/1‑150 (MQ=255)
                                                                                                                                          gggTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTaa   <  2:3340891/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAAC  >  CP009273/2209017‑2209305

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: