Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,452,566 C→A 44.3% T290K (ACG→AAG)  gspE → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,452,5660CA44.3% 216.9 / 301.0 238T290K (ACG→AAG) gspEgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (66/66);  new base A (59/46);  total (126/112)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

TCCGCACGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTA  >  CP009273/3452417‑3452714
                                                                                                                                                     |                                                                                                                                                    
  cGCACGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCACTATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGc                                                                                                                                                    <  1:1503804/150‑1 (MQ=255)
   gCACGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCa                                                                                                                                                   <  2:6894/150‑1 (MQ=255)
    cACGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAg                                                                                                                                                  <  1:32748/150‑1 (MQ=255)
         tattatCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGaa                                                                                                                                             <  2:4001458/150‑1 (MQ=255)
         tattatCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGaa                                                                                                                                             <  1:3868064/150‑1 (MQ=255)
         tattatCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATAGTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGaa                                                                                                                                             <  2:3651703/150‑1 (MQ=255)
          attatCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAAt                                                                                                                                            <  1:892635/150‑1 (MQ=255)
           ttatCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATa                                                                                                                                           <  1:3829084/150‑1 (MQ=255)
            tatCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATAc                                                                                                                                          <  2:1726702/150‑1 (MQ=255)
            tatCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATAc                                                                                                                                          >  2:2009611/1‑150 (MQ=255)
               ccTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCg                                                                                                                                       >  2:4563139/1‑150 (MQ=255)
                     gACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGa                                                                                                                                 >  2:1974679/1‑150 (MQ=255)
                     gACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGa                                                                                                                                 >  2:3858597/1‑150 (MQ=255)
                     gACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGa                                                                                                                                 >  2:4064502/1‑150 (MQ=255)
                        aGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGatat                                                                                                                              <  2:1152184/150‑1 (MQ=255)
                         gggCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATg                                                                                                                             >  2:2581741/1‑150 (MQ=255)
                          ggCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGt                                                                                                                            <  1:500078/150‑1 (MQ=255)
                             cGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGACTATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCtt                                                                                                                         >  2:2948608/1‑150 (MQ=255)
                             cGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCtt                                                                                                                         >  2:1532971/1‑150 (MQ=255)
                             cGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCtt                                                                                                                         >  1:4326921/1‑150 (MQ=255)
                               aCAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTtc                                                                                                                       >  2:847944/1‑150 (MQ=255)
                                 aGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTtcgc                                                                                                                     >  1:3475469/1‑150 (MQ=255)
                                     tCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTtcgctcgc                                                                                                                 >  1:2443953/1‑150 (MQ=255)
                                      ccGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTcgcg                                                                                                                <  1:1058425/150‑1 (MQ=255)
                                       cGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCgg                                                                                                               <  1:1168737/150‑1 (MQ=255)
                                        ggTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGc                                                                                                              >  1:4163488/1‑150 (MQ=255)
                                         gTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGcc                                                                                                             >  1:3154542/1‑150 (MQ=255)
                                          tAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCt                                                                                                            >  1:2828355/1‑150 (MQ=255)
                                           aaaaGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTg                                                                                                           <  1:4563447/150‑1 (MQ=255)
                                           aaaaGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTg                                                                                                           >  1:2720478/1‑150 (MQ=255)
                                              aGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTgcgc                                                                                                        >  1:2589456/1‑150 (MQ=255)
                                                caccacGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTgcgcgc                                                                                                      >  2:3986900/1‑150 (MQ=255)
                                                caccacGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTgcgcgc                                                                                                      >  2:1917523/1‑150 (MQ=255)
                                                caccacGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTgcgcgc                                                                                                      >  2:1694932/1‑150 (MQ=255)
                                                caccacGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTgcgcgc                                                                                                      >  2:1080460/1‑150 (MQ=255)
                                                 accacGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGcc                                                                                                     >  2:4212801/1‑150 (MQ=255)
                                                   cacGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCAt                                                                                                   <  1:1676041/150‑1 (MQ=255)
                                                      gCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACt                                                                                                >  2:1662488/1‑150 (MQ=255)
                                                         ctACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCg                                                                                             <  1:2089761/150‑1 (MQ=255)
                                                           aCGCCATCCTTTCGGCGCTGAATACTCCCGGCCACAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGcc                                                                                           <  1:2275605/150‑1 (MQ=255)
                                                            cGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCa                                                                                          >  2:3333170/1‑150 (MQ=255)
                                                            cGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCa                                                                                          <  1:972182/150‑1 (MQ=255)
                                                            cGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCa                                                                                          <  1:2990454/150‑1 (MQ=255)
                                                              ccATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAgg                                                                                        >  2:3312094/1‑150 (MQ=255)
                                                                aTCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGAc                                                                                      >  2:665198/1‑150 (MQ=255)
                                                                aTCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGAc                                                                                      <  1:2289005/150‑1 (MQ=255)
                                                                  ccTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGAccc                                                                                    >  1:1699814/1‑150 (MQ=255)
                                                                  ccTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGAccc                                                                                    >  2:2535031/1‑150 (MQ=255)
                                                                   cTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCg                                                                                   <  1:3042645/150‑1 (MQ=255)
                                                                   cTTTCGGCGCTGAATACTCCCGGCCGCAAGATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCg                                                                                   <  2:3096119/150‑1 (MQ=255)
                                                                    tttCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCgg                                                                                  <  1:706703/150‑1 (MQ=255)
                                                                      tCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGAt                                                                                <  2:1096510/150‑1 (MQ=255)
                                                                        ggCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATgt                                                                              >  2:4396137/1‑150 (MQ=255)
                                                                           gcTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGGGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATggcgt                                                                           >  2:1481882/1‑150 (MQ=255)
                                                                              gAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAt                                                                        <  1:2225718/150‑1 (MQ=255)
                                                                              gAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAt                                                                        <  1:3377479/150‑1 (MQ=255)
                                                                                  aCTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtggtg                                                                    >  2:4088876/1‑150 (MQ=255)
                                                                                    tCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggg                                                                  <  1:1650483/150‑1 (MQ=255)
                                                                                     cccGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggg                                                                 >  1:4236079/1‑150 (MQ=255)
                                                                                     cccGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggg                                                                 >  2:3701975/1‑150 (MQ=255)
                                                                                     cccGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggg                                                                 <  2:517872/150‑1 (MQ=255)
                                                                                          ccGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGCCCCGGATGTCGTCATGGTGGGGGAAAt                                                            >  2:2487887/1‑150 (MQ=255)
                                                                                             cAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAGGGTGGGGGAAATTCg                                                         >  2:249968/1‑150 (MQ=255)
                                                                                               atatTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTg                                                       >  1:4589544/1‑150 (MQ=255)
                                                                                                tatTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGa                                                      >  2:2315761/1‑150 (MQ=255)
                                                                                                 atTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGAt                                                     >  1:2417582/1‑150 (MQ=255)
                                                                                                 atTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGAt                                                     >  1:311634/1‑150 (MQ=255)
                                                                                                   tCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATAc                                                   <  2:2237067/150‑1 (MQ=255)
                                                                                                    cTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACa                                                  <  1:2560690/150‑1 (MQ=255)
                                                                                                        cGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGaaa                                              <  2:1738144/150‑1 (MQ=255)
                                                                                                           tAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACcg                                           <  1:21652/150‑1 (MQ=255)
                                                                                                            agaagaTCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACcgc                                          >  2:1349756/1‑150 (MQ=255)
                                                                                                            agaagaTCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACcgc                                          <  1:3231745/150‑1 (MQ=255)
                                                                                                             gaagaTCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACcgcg                                         >  1:658992/1‑150 (MQ=255)
                                                                                                              aagaTCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCgcgc                                        >  2:2383321/1‑150 (MQ=255)
                                                                                                              aagaTCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCgcgc                                        >  2:1770077/1‑150 (MQ=255)
                                                                                                                  tCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGAt                                    <  2:2082972/150‑1 (MQ=255)
                                                                                                                     cGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGc                                 >  2:1872792/1‑150 (MQ=255)
                                                                                                                         gAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGtt                             >  1:1365988/1‑150 (MQ=255)
                                                                                                                         gAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGtt                             <  2:2327776/150‑1 (MQ=255)
                                                                                                                           atatGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCa                           <  1:1672021/150‑1 (MQ=255)
                                                                                                                               gAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCATGACACGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGcc                       <  2:4466736/150‑1 (MQ=255)
                                                                                                                               gAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGcc                       <  2:4370272/150‑1 (MQ=255)
                                                                                                                                aGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACGTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCt                      <  2:1190201/150‑1 (MQ=255)
                                                                                                                                 gCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTc                     <  1:2652783/150‑1 (MQ=255)
                                                                                                                                 gCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTc                     <  1:328187/150‑1 (MQ=255)
                                                                                                                                     gAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTTGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTc                 <  1:1902887/150‑1 (MQ=255)
                                                                                                                                     gAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTc                 >  2:553377/1‑150 (MQ=255)
                                                                                                                                      aaGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCa                <  1:3536248/150‑1 (MQ=255)
                                                                                                                                      aaGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCa                <  2:2690128/150‑1 (MQ=255)
                                                                                                                                      aaGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCTGGCCTCGCTCa                >  2:3515966/1‑150 (MQ=255)
                                                                                                                                           aTTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGc           >  1:913621/1‑150 (MQ=255)
                                                                                                                                           aTTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGc           >  1:715532/1‑150 (MQ=255)
                                                                                                                                             tGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCa         <  2:4505518/150‑1 (MQ=255)
                                                                                                                                               ggCCAACGCAGGCGAATACCCGTGTGGATATGTCTGTCGCTCGCGGCCAGCGCGCCCTACGCCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGACACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCCCCGGCCATc       <  2:329415/150‑1 (MQ=255)
                                                                                                                                               ggCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATc       <  2:2175614/150‑1 (MQ=255)
                                                                                                                                                gCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCt      <  1:790430/150‑1 (MQ=255)
                                                                                                                                                 cAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTg     >  1:827946/1‑150 (MQ=255)
                                                                                                                                                    aCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTa  >  2:211783/1‑150 (MQ=255)
                                                                                                                                                    aCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTa  >  2:3812157/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCCGCACGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTA  >  CP009273/3452417‑3452714

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: