Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,501,734 G→T 26.0% D233Y (GAT→TAT)  ydcP → putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,501,7340GT26.0% 511.7 / 166.5 230D233Y (GAT→TAT) ydcPputative peptidase
Reads supporting (aligned to +/- strand):  ref base G (90/80);  new base T (27/33);  total (117/113)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

AGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACAT  >  CP009273/1501586‑1501883
                                                                                                                                                    |                                                                                                                                                     
aGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCtgattga                                                                                                                                                      >  1:585293/1‑150 (MQ=255)
  ggCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATg                                                                                                                                                    <  2:2899498/150‑1 (MQ=255)
  ggCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATg                                                                                                                                                    >  2:4082430/1‑150 (MQ=255)
      tAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTgg                                                                                                                                                >  1:4122136/1‑150 (MQ=255)
      tAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTgg                                                                                                                                                >  1:1613346/1‑150 (MQ=255)
          gcCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTa                                                                                                                                            <  1:1283532/150‑1 (MQ=255)
          gcCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTa                                                                                                                                            <  1:3147029/150‑1 (MQ=255)
          gcCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTa                                                                                                                                            <  1:3147875/150‑1 (MQ=255)
           ccAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTAc                                                                                                                                           <  1:4314344/150‑1 (MQ=255)
            cAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACg                                                                                                                                          >  1:4251333/1‑150 (MQ=255)
                cGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTcc                                                                                                                                      >  1:4782855/1‑150 (MQ=255)
                 gTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGACAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCt                                                                                                                                     <  2:234085/150‑1 (MQ=255)
                 gTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCt                                                                                                                                     >  1:78199/1‑150 (MQ=255)
                  tGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCtt                                                                                                                                    <  1:3488783/150‑1 (MQ=255)
                        ttGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGAt                                                                                                                              >  2:3340562/1‑150 (MQ=255)
                            tCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGaa                                                                                                                          <  1:3381125/150‑1 (MQ=255)
                              gCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAgg                                                                                                                        >  1:4040358/1‑150 (MQ=255)
                                aGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGc                                                                                                                      <  1:368112/150‑1 (MQ=255)
                                    gTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTa                                                                                                                  >  1:2948144/1‑150 (MQ=255)
                                         gTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAg                                                                                                             >  1:4462146/1‑150 (MQ=255)
                                           ttGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGat                                                                                                           <  1:2605639/150‑1 (MQ=255)
                                            tGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGata                                                                                                          <  1:2796559/150‑1 (MQ=255)
                                             gCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGatat                                                                                                         <  1:1937516/150‑1 (MQ=255)
                                                aTACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                      >  1:4798762/1‑150 (MQ=255)
                                                 tACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGc                                                                                                     <  1:2250997/150‑1 (MQ=255)
                                                 tACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGc                                                                                                     <  1:173660/150‑1 (MQ=255)
                                                  acacATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCt                                                                                                    >  2:175091/1‑150 (MQ=255)
                                                   cacaTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTa                                                                                                   <  1:3614458/150‑1 (MQ=255)
                                                        tGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTga                                                                                              <  2:3235272/150‑1 (MQ=255)
                                                         gAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                             >  1:4473718/1‑150 (MQ=255)
                                                           aaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaaga                                                                                           <  2:758488/150‑1 (MQ=255)
                                                           aaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaaga                                                                                           <  2:3974945/150‑1 (MQ=255)
                                                             gACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAat                                                                                         <  2:2027721/150‑1 (MQ=255)
                                                             gACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAat                                                                                         <  2:3946201/150‑1 (MQ=255)
                                                             gACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAat                                                                                         <  2:735668/150‑1 (MQ=255)
                                                                  tCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTTAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATAAc                                                                                    >  1:4868160/1‑150 (MQ=255)
                                                                  tCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCAc                                                                                    >  2:2512339/1‑150 (MQ=255)
                                                                   cAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCAcc                                                                                   >  2:670234/1‑150 (MQ=255)
                                                                   cAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCAcc                                                                                   <  1:4957994/150‑1 (MQ=255)
                                                                   cAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCAcc                                                                                   >  1:1411260/1‑150 (MQ=255)
                                                                    aGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCAGCg                                                                                  >  2:4249152/1‑150 (MQ=255)
                                                                    aGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCg                                                                                  >  2:573615/1‑150 (MQ=255)
                                                                    aGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCg                                                                                  >  1:1930913/1‑150 (MQ=255)
                                                                      gggCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGcc                                                                                >  2:3397596/1‑150 (MQ=255)
                                                                         cGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCAt                                                                             >  1:4884914/1‑150 (MQ=255)
                                                                         cGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCAt                                                                             >  2:4875157/1‑150 (MQ=255)
                                                                         cGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCAt                                                                             >  2:1584493/1‑150 (MQ=255)
                                                                          gggTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCAtt                                                                            <  1:2659271/150‑1 (MQ=255)
                                                                             tggtTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATc                                                                         <  1:1074581/150‑1 (MQ=255)
                                                                             tggtTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATc                                                                         >  1:83665/1‑150 (MQ=255)
                                                                             tggtTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATc                                                                         >  1:1618932/1‑150 (MQ=255)
                                                                                 ttCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCa                                                                     >  2:239216/1‑150 (MQ=255)
                                                                                    cTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGAt                                                                  <  1:1465994/150‑1 (MQ=255)
                                                                                    cTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGAt                                                                  <  2:2989189/150‑1 (MQ=255)
                                                                                    cTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGAt                                                                  <  1:1085958/150‑1 (MQ=255)
                                                                                     tATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATg                                                                 >  2:1464980/1‑150 (MQ=255)
                                                                                           aaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGAt                                                           <  2:419019/150‑1 (MQ=255)
                                                                                           aaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGAt                                                           <  2:1100700/150‑1 (MQ=255)
                                                                                            aaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATg                                                          >  1:2983432/1‑150 (MQ=255)
                                                                                            aaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATg                                                          <  1:917556/150‑1 (MQ=255)
                                                                                              cATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGcc                                                        <  1:586887/150‑1 (MQ=255)
                                                                                                 ctgctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                     <  1:3587024/150‑1 (MQ=255)
                                                                                                 ctgctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                     <  1:3671005/150‑1 (MQ=255)
                                                                                                  tgctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatta                                                    >  1:1834435/1‑150 (MQ=255)
                                                                                                   gctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAATGATATGAGCTACGAGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                   <  2:2782609/150‑1 (MQ=255)
                                                                                                    ctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattatt                                                  >  2:3136945/1‑150 (MQ=255)
                                                                                                    ctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattatt                                                  >  2:1499785/1‑150 (MQ=255)
                                                                                                     tgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTg                                                 >  1:470461/1‑150 (MQ=255)
                                                                                                      gTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCAGTATTga                                                >  2:2018915/1‑150 (MQ=255)
                                                                                                       tCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaa                                               >  2:452505/1‑150 (MQ=255)
                                                                                                       tCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaa                                               >  1:2409239/1‑150 (MQ=255)
                                                                                                         gATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaaga                                             >  1:4006340/1‑150 (MQ=255)
                                                                                                          aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGATCTACGTGAAGAATATCACCTCCCATTATCGCCAGATGCTTGATGCCATTATTgaagaa                                            <  2:3180385/150‑1 (MQ=255)
                                                                                                          aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaagaa                                            >  1:1260000/1‑150 (MQ=255)
                                                                                                           tGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAAc                                           >  1:4637200/1‑150 (MQ=255)
                                                                                                            gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACg                                          <  2:3722694/150‑1 (MQ=255)
                                                                                                            gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCGACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACg                                          <  2:2011510/150‑1 (MQ=255)
                                                                                                             aaaGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGt                                         >  1:2473966/1‑150 (MQ=255)
                                                                                                               aGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTgg                                       <  1:3097525/150‑1 (MQ=255)
                                                                                                                 aTAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCg                                     <  2:2621051/150‑1 (MQ=255)
                                                                                                                  tAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGa                                    <  1:2025263/150‑1 (MQ=255)
                                                                                                                     cGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCt                                 <  2:1069182/150‑1 (MQ=255)
                                                                                                                          aGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGgcgc                            >  1:2731085/1‑150 (MQ=255)
                                                                                                                          aGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGgcgc                            >  1:5064606/1‑150 (MQ=255)
                                                                                                                             cTGCCAACCTCGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGgcgcgcg                         >  1:3046042/1‑150 (MQ=255)
                                                                                                                             cTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGACGAACGTGGCGATCTGgcgcgcg                         >  2:13427/1‑150 (MQ=255)
                                                                                                                               gCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCt                       >  1:410605/1‑150 (MQ=255)
                                                                                                                                ccAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGTGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCtt                      >  2:548496/1‑150 (MQ=255)
                                                                                                                                  aaCCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTtca                    >  1:669714/1‑150 (MQ=255)
                                                                                                                                     cTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTtcatca                 >  1:3276245/1‑150 (MQ=255)
                                                                                                                                       cGGCGCGCTGATTTATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAgg               <  1:2298682/150‑1 (MQ=255)
                                                                                                                                           gcgcTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGt           <  2:2247597/150‑1 (MQ=255)
                                                                                                                                              cTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACt        >  2:2192783/1‑150 (MQ=255)
                                                                                                                                               tgatttatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTg       <  1:3870830/150‑1 (MQ=255)
                                                                                                                                               tgatttatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTg       >  2:1240593/1‑150 (MQ=255)
                                                                                                                                               tgattgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTg       <  1:4082419/150‑1 (MQ=255)
                                                                                                                                                  ttgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAAc    >  1:5125628/1‑150 (MQ=255)
                                                                                                                                                  ttgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAAc    <  1:2003751/150‑1 (MQ=255)
                                                                                                                                                    gatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACAt  <  1:4189573/150‑1 (MQ=255)
                                                                                                                                                    gatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACAt  <  2:4789625/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGGGCGTAGCGCCAACCGTGGCGATTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACAT  >  CP009273/1501586‑1501883

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: ‑