Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,020,166 C→A 61.9% D445Y (GAC→TAC)  yedQ ← putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,020,1660CA61.9% 157.5 / 316.6 239D445Y (GAC→TAC) yedQputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (32/59);  new base A (74/74);  total (106/133)
Fisher's exact test for biased strand distribution p-value = 3.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 243 total aligned reads displayed.

AGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCA  >  CP009273/2020020‑2020313
                                                                                                                                                  |                                                                                                                                                   
aGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAAt                                                                                                                                                   <  1:2613953/149‑1 (MQ=255)
 gCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCAtt                                                                                                                                                  <  1:1794172/149‑1 (MQ=255)
  ccTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTa                                                                                                                                                 <  1:4053023/149‑1 (MQ=255)
    tCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGtaattaat                                                                                                                                               >  2:1899965/1‑149 (MQ=255)
    tCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAAt                                                                                                                                               <  2:3745430/149‑1 (MQ=255)
     ccGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATc                                                                                                                                              <  1:3155556/149‑1 (MQ=255)
      cGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCg                                                                                                                                             >  1:3573989/1‑149 (MQ=255)
        tCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCt                                                                                                                                           <  2:2055855/149‑1 (MQ=255)
           gACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTa                                                                                                                                        <  1:1276668/149‑1 (MQ=255)
           gACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTa                                                                                                                                        <  1:678209/149‑1 (MQ=255)
           gACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTa                                                                                                                                        <  1:2917414/149‑1 (MQ=255)
            aCTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTaa                                                                                                                                       >  2:2525413/1‑149 (MQ=255)
            aCTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTaa                                                                                                                                       <  2:3490332/149‑1 (MQ=255)
            aCTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTaa                                                                                                                                       >  2:3177079/1‑149 (MQ=255)
            aCTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTaa                                                                                                                                       <  2:4100201/149‑1 (MQ=255)
                  gcCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGACCTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGt                                                                                                                                 >  2:4590694/1‑149 (MQ=255)
                     tGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCaa                                                                                                                              >  1:4281613/1‑149 (MQ=255)
                     tGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCaa                                                                                                                              <  1:2406050/149‑1 (MQ=255)
                      ggCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAg                                                                                                                             >  2:4405120/1‑149 (MQ=255)
                      ggCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAg                                                                                                                             <  2:1093767/149‑1 (MQ=255)
                        cAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGAt                                                                                                                           >  1:1317821/1‑149 (MQ=255)
                                   aaaaCTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATg                                                                                                                >  2:4312933/1‑149 (MQ=255)
                                   aaaaCTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATg                                                                                                                >  2:4637186/1‑149 (MQ=255)
                                       ctcctcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAg                                                                                                            >  2:313365/1‑149 (MQ=255)
                                            caccacCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGaaaaa                                                                                                       <  2:1416386/149‑1 (MQ=255)
                                             accaccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAg                                                                                                      >  2:3849842/1‑149 (MQ=255)
                                             accaccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAg                                                                                                      <  1:4598883/149‑1 (MQ=255)
                                             accaccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAg                                                                                                      >  1:1881189/1‑149 (MQ=255)
                                              ccaccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAAGTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAgg                                                                                                     >  1:4120549/1‑149 (MQ=255)
                                               caccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAACAAggc                                                                                                    >  2:1691002/1‑148 (MQ=255)
                                                accGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAGGATg                                                                                                  >  1:2554175/1‑149 (MQ=255)
                                                accGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGAt                                                                                                   <  2:2722677/149‑1 (MQ=255)
                                                accGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGAt                                                                                                   <  1:2120163/149‑1 (MQ=255)
                                                accGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGAt                                                                                                   <  1:1303909/149‑1 (MQ=255)
                                                 ccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATg                                                                                                  >  1:1106405/1‑149 (MQ=255)
                                                 ccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATg                                                                                                  >  1:2488094/1‑149 (MQ=255)
                                                  cGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGt                                                                                                 <  1:1859041/149‑1 (MQ=255)
                                                    aCCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTg                                                                                               <  2:458019/149‑1 (MQ=255)
                                                    aCCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTg                                                                                               >  2:3545506/1‑149 (MQ=255)
                                                        gcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGc                                                                                           <  1:4190468/149‑1 (MQ=255)
                                                             gCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTg                                                                                      <  1:2482736/149‑1 (MQ=255)
                                                               aaCGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGTCAGCATGAGAAAGAACACTGTCGCCCGCATGATGACCAAATCGGTAATTAATCGATTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTATGac                                                                                    <  1:2288831/149‑1 (MQ=255)
                                                                   tCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACaa                                                                                <  2:1583538/149‑1 (MQ=255)
                                                                   tCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACaa                                                                                >  2:2258380/1‑149 (MQ=255)
                                                                     cTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAAtt                                                                              <  2:3689140/149‑1 (MQ=255)
                                                                           acgcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcga                                                                        <  2:3860241/149‑1 (MQ=255)
                                                                             gcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcgagc                                                                      <  1:3266098/149‑1 (MQ=255)
                                                                              caAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcgagcg                                                                     >  2:2305208/1‑149 (MQ=255)
                                                                               aaGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCgg                                                                    >  2:1509112/1‑149 (MQ=255)
                                                                                aGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGa                                                                   >  2:2234092/1‑149 (MQ=255)
                                                                                aGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGa                                                                   <  2:3010671/149‑1 (MQ=255)
                                                                                 ggAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGAc                                                                  <  1:1498053/149‑1 (MQ=255)
                                                                                 ggAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGAc                                                                  <  2:2540491/149‑1 (MQ=255)
                                                                                   aaCTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACgg                                                                <  2:4318870/149‑1 (MQ=255)
                                                                                     cTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGc                                                              >  1:2191604/1‑149 (MQ=255)
                                                                                      tGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCt                                                             <  2:1180894/149‑1 (MQ=255)
                                                                                       gCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCtt                                                            <  1:137899/149‑1 (MQ=255)
                                                                                         taattaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCtttt                                                          <  1:4618144/149‑1 (MQ=255)
                                                                                         taattaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCtttt                                                          <  1:3611541/149‑1 (MQ=255)
                                                                                          aattaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCttttt                                                         >  1:4638926/1‑149 (MQ=255)
                                                                                           attaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTc                                                        >  1:749755/1‑149 (MQ=255)
                                                                                              aatCCGGCAGCATGAGAAAGACCACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGaa                                                     <  2:330486/149‑1 (MQ=255)
                                                                                                tCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACCATTTAGCGAGCGGACGGGCTTTTTCGAAcc                                                   >  2:3066010/1‑148 (MQ=255)
                                                                                                tCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACa                                                   >  1:4491807/1‑149 (MQ=255)
                                                                                                 ccGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAg                                                  <  2:1894480/149‑1 (MQ=255)
                                                                                                  cGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGt                                                 >  1:2755854/1‑149 (MQ=255)
                                                                                                   ggCAGAATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTg                                                <  1:3979574/149‑1 (MQ=255)
                                                                                                    gcagcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTgc                                               >  1:2675637/1‑149 (MQ=255)
                                                                                                    gcagcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTgc                                               >  2:1158717/1‑149 (MQ=255)
                                                                                                      agcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTgcgc                                             >  2:4400621/1‑149 (MQ=255)
                                                                                                       gcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGcc                                            <  2:4669182/149‑1 (MQ=255)
                                                                                                       gcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTg                                                <  2:2636816/145‑1 (MQ=255)
                                                                                                       gcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGcc                                            <  1:1893317/149‑1 (MQ=255)
                                                                                                         aTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCAc                                          <  2:4249919/149‑1 (MQ=255)
                                                                                                          tGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACg                                         >  1:1949083/1‑149 (MQ=255)
                                                                                                            agaaagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGt                                       <  2:1596394/149‑1 (MQ=255)
                                                                                                            agaaagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGt                                       >  1:4461861/1‑149 (MQ=255)
                                                                                                            agaaagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGt                                       >  2:3614502/1‑149 (MQ=255)
                                                                                                             gaaagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGtt                                      >  1:2776720/1‑149 (MQ=255)
                                                                                                                agaaCCCGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtata                                   <  1:4397395/149‑1 (MQ=255)
                                                                                                                agaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtata                                   >  1:1950609/1‑149 (MQ=255)
                                                                                                                 gaaCACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtatat                                  >  2:1371103/1‑149 (MQ=255)
                                                                                                                  aaCACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtatata                                 <  1:2525413/149‑1 (MQ=255)
                                                                                                                   acacGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATaa                                >  1:3889130/1‑149 (MQ=255)
                                                                                                                    cacGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATaaa                               <  2:922339/149‑1 (MQ=255)
                                                                                                                    cacGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATaaa                               <  2:869778/149‑1 (MQ=255)
                                                                                                                         tCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGt                          <  1:1481671/149‑1 (MQ=255)
                                                                                                                         tCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGt                          <  2:321290/149‑1 (MQ=255)
                                                                                                                            cccgccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTaa                       <  2:1943579/149‑1 (MQ=255)
                                                                                                                            cccgccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTaa                       >  2:3907852/1‑149 (MQ=255)
                                                                                                                             ccgccTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAg                      <  2:2418834/149‑1 (MQ=255)
                                                                                                                             ccgccTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAg                      <  2:2418980/149‑1 (MQ=255)
                                                                                                                                 cTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgtgt                  <  1:2365842/149‑1 (MQ=255)
                                                                                                                                 cTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgtgt                  <  1:1036375/149‑1 (MQ=255)
                                                                                                                                      gaCCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGc             >  1:1757647/1‑149 (MQ=255)
                                                                                                                                       aCCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGcc            >  2:1548946/1‑149 (MQ=255)
                                                                                                                                         cAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCAc          >  2:2698157/1‑149 (MQ=255)
                                                                                                                                           aaGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGc        <  1:1579642/149‑1 (MQ=255)
                                                                                                                                            aGCTGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGcc       <  2:4464145/149‑1 (MQ=255)
                                                                                                                                            aGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGcc       >  1:2083108/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCA  >  CP009273/2020020‑2020313

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: