Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,452,566 C→A 9.5% 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,5660CA 9.5% 877.9 / 45.2 236T290K (ACG→AAG) gspEgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (110/103);  new base A (12/11);  total (122/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

TCCGCACGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTAC  >  CP009273/3452417‑3452715
                                                                                                                                                     |                                                                                                                                                     
tCCGCACGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAc                                                                                                                                                       <  1:3261125/150‑1 (MQ=255)
 ccGCACGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGATGGAAGGCATTGGGCAAACg                                                                                                                                                      <  1:1019831/150‑1 (MQ=255)
  cGCACGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGc                                                                                                                                                     <  1:846989/150‑1 (MQ=255)
    cACGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAg                                                                                                                                                   >  1:2388818/1‑150 (MQ=255)
     aCGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAgg                                                                                                                                                  <  2:702846/150‑1 (MQ=255)
      cGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGt                                                                                                                                                 <  2:2244081/150‑1 (MQ=255)
      cGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGt                                                                                                                                                 <  2:2022996/150‑1 (MQ=255)
       ggTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTg                                                                                                                                                >  2:3906445/1‑150 (MQ=255)
       ggTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTg                                                                                                                                                >  1:3301513/1‑150 (MQ=255)
        gTATTATCCTGGTGACAGGGCCGACAGGCGCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGa                                                                                                                                               >  2:359249/1‑150 (MQ=255)
          attatCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAAt                                                                                                                                             >  2:3504089/1‑150 (MQ=255)
          attatCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAAt                                                                                                                                             >  1:137690/1‑150 (MQ=255)
            tatCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATAc                                                                                                                                           <  2:1646810/150‑1 (MQ=255)
               ccTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCg                                                                                                                                        >  2:309552/1‑150 (MQ=255)
               ccTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCg                                                                                                                                        >  1:30249/1‑150 (MQ=255)
               ccTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCg                                                                                                                                        <  2:370079/150‑1 (MQ=255)
                   gTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGtgtg                                                                                                                                    <  1:1713273/150‑1 (MQ=255)
                    tGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTgg                                                                                                                                   >  2:1208977/1‑150 (MQ=255)
                    tGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTgg                                                                                                                                   >  2:1559437/1‑150 (MQ=255)
                     gACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGa                                                                                                                                  >  1:3407108/1‑150 (MQ=255)
                      aCAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGat                                                                                                                                 >  2:477604/1‑150 (MQ=255)
                       cAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGata                                                                                                                                >  1:686012/1‑150 (MQ=255)
                       cAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGata                                                                                                                                >  2:1883687/1‑150 (MQ=255)
                         gggCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATg                                                                                                                              <  1:3954767/150‑1 (MQ=255)
                            ccGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCt                                                                                                                           <  1:1192371/150‑1 (MQ=255)
                             cGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCtt                                                                                                                          <  1:2758702/150‑1 (MQ=255)
                              gACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCttt                                                                                                                         <  1:2805428/150‑1 (MQ=255)
                              gACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCttt                                                                                                                         >  2:600328/1‑150 (MQ=255)
                               aCAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTtc                                                                                                                        <  1:3158902/150‑1 (MQ=255)
                               aCAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTtc                                                                                                                        <  1:2724257/150‑1 (MQ=255)
                                cAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTtcg                                                                                                                       <  2:2593300/150‑1 (MQ=255)
                                cAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTtcg                                                                                                                       <  1:3465053/150‑1 (MQ=255)
                                  ggCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTtcgct                                                                                                                     <  1:1916941/150‑1 (MQ=255)
                                      ccGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTcgcg                                                                                                                 >  1:1307480/1‑150 (MQ=255)
                                       cGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCgg                                                                                                                >  2:3106754/1‑150 (MQ=255)
                                       cGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCgg                                                                                                                >  1:2950687/1‑150 (MQ=255)
                                       cGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTACGCTCGCgg                                                                                                                >  1:860735/1‑150 (MQ=255)
                                         gTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGcc                                                                                                              <  2:1634762/150‑1 (MQ=255)
                                         gTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGcc                                                                                                              <  1:1692870/150‑1 (MQ=255)
                                           aaaaGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTg                                                                                                            <  1:1080168/150‑1 (MQ=255)
                                               gCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTgcgcg                                                                                                        <  2:1834155/150‑1 (MQ=255)
                                                caccacGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTgcgcgc                                                                                                       >  1:1544950/1‑150 (MQ=255)
                                                 accacGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGcc                                                                                                      >  1:488215/1‑150 (MQ=255)
                                                    acGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATa                                                                                                   <  2:761615/150‑1 (MQ=255)
                                                    acGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATa                                                                                                   <  1:3130332/150‑1 (MQ=255)
                                                    acGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATa                                                                                                   >  2:3361393/1‑150 (MQ=255)
                                                     cGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATAc                                                                                                  >  2:3831313/1‑150 (MQ=255)
                                                      gCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACt                                                                                                 >  2:3916850/1‑150 (MQ=255)
                                                        tctACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTc                                                                                               <  2:1160850/150‑1 (MQ=255)
                                                        tctACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGAAGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCTTACTTc                                                                                               <  2:2379046/150‑1 (MQ=255)
                                                          tACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGc                                                                                             >  1:3916156/1‑150 (MQ=255)
                                                           aCGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGcc                                                                                            >  1:3764135/1‑150 (MQ=255)
                                                            cGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCa                                                                                           <  1:954909/150‑1 (MQ=255)
                                                             gCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAg                                                                                          >  2:1397770/1‑150 (MQ=255)
                                                              ccATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAgg                                                                                         <  2:1312542/150‑1 (MQ=255)
                                                              ccATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAgg                                                                                         <  1:2144557/150‑1 (MQ=255)
                                                                aTCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAGGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGAc                                                                                       <  2:1036486/150‑1 (MQ=255)
                                                                   cTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCg                                                                                    >  2:3954039/1‑150 (MQ=255)
                                                                      tCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTTTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGAt                                                                                 >  1:1622455/1‑150 (MQ=255)
                                                                       cGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATg                                                                                <  2:3335245/150‑1 (MQ=255)
                                                                         gcgcTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATgtc                                                                              <  2:2449569/150‑1 (MQ=255)
                                                                          cgcTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATgtcg                                                                             >  2:315346/1‑150 (MQ=255)
                                                                            cTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATgtcgtc                                                                           <  2:4075567/150‑1 (MQ=255)
                                                                            cTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATgtcgtc                                                                           >  1:3976665/1‑150 (MQ=255)
                                                                            cTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATgtcgtc                                                                           >  1:2607408/1‑150 (MQ=255)
                                                                              gAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAt                                                                         <  2:2500672/150‑1 (MQ=255)
                                                                               aaTACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtg                                                                        <  2:2456410/150‑1 (MQ=255)
                                                                               aaTACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtg                                                                        >  2:2269248/1‑150 (MQ=255)
                                                                                aTACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtgg                                                                       <  2:1575314/150‑1 (MQ=255)
                                                                                aTACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtgg                                                                       >  2:2108795/1‑150 (MQ=255)
                                                                                aTACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtgg                                                                       <  2:3800131/150‑1 (MQ=255)
                                                                                  aCTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtggtg                                                                     >  2:2940270/1‑150 (MQ=255)
                                                                                  aCTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtggtg                                                                     <  2:647289/150‑1 (MQ=255)
                                                                                        ggCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGaa                                                               >  1:765891/1‑150 (MQ=255)
                                                                                        ggCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGaa                                                               >  1:1130183/1‑150 (MQ=255)
                                                                                        ggCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGaa                                                               >  2:3259110/1‑150 (MQ=255)
                                                                                        ggCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGaa                                                               >  2:3143717/1‑150 (MQ=255)
                                                                                        ggCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGaa                                                               >  1:220711/1‑150 (MQ=255)
                                                                                               atatTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTg                                                        >  2:1993888/1‑150 (MQ=255)
                                                                                                     tGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAg                                                  <  1:1628932/150‑1 (MQ=255)
                                                                                                           tAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACcg                                            >  1:3549100/1‑150 (MQ=255)
                                                                                                             gaagaTCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACcgcg                                          >  1:1535597/1‑150 (MQ=255)
                                                                                                               agaTCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCa                                        >  2:2732199/1‑150 (MQ=255)
                                                                                                                 aTCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGa                                      >  1:272844/1‑150 (MQ=255)
                                                                                                                 aTCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGa                                      >  1:454907/1‑150 (MQ=255)
                                                                                                                   cccGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGAtt                                    <  2:3508385/150‑1 (MQ=255)
                                                                                                                   cccGTGGAATATGAGCTGGAAGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGAtt                                    >  1:3825928/1‑150 (MQ=255)
                                                                                                                        ggAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGt                               >  1:4131245/1‑150 (MQ=255)
                                                                                                                          aaTATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTc                             <  2:1884362/150‑1 (MQ=255)
                                                                                                                          aaTATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTc                             <  1:2742250/150‑1 (MQ=255)
                                                                                                                           atatGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCa                            <  2:2357631/150‑1 (MQ=255)
                                                                                                                           atatGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCa                            <  1:1775948/150‑1 (MQ=255)
                                                                                                                              tGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGc                         <  2:1242979/150‑1 (MQ=255)
                                                                                                                                      aaGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCa                 >  2:2603076/1‑150 (MQ=255)
                                                                                                                                      aaGGCATTGGGCAAAAGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCa                 <  2:879085/150‑1 (MQ=255)
                                                                                                                                        ggCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCAcc               <  1:2895970/150‑1 (MQ=255)
                                                                                                                                           aTTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGc            >  2:1757736/1‑150 (MQ=255)
                                                                                                                                             tGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCa          <  1:3234194/150‑1 (MQ=255)
                                                                                                                                               ggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATc        >  2:3818044/1‑150 (MQ=255)
                                                                                                                                                    aCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTa   <  2:1004866/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCGCACGGTATTATCCTGGTGACAGGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTAC  >  CP009273/3452417‑3452715

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: