Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,020,166 C→A 94.4% 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,1660CA94.4% 722.2 / 30.0 232D445Y (GAC→TAC) yedQputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (7/6);  new base A (97/122);  total (104/128)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 233 total aligned reads displayed.

CAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTG  >  CP009273/2020019‑2020310
                                                                                                                                                   |                                                                                                                                                
cAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAAt                                                                                                                                                <  2:142603/150‑1 (MQ=255)
 aGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAAtt                                                                                                                                               <  1:3093708/150‑1 (MQ=255)
 aGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAAtt                                                                                                                                               <  1:3093234/150‑1 (MQ=255)
 aGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAAtt                                                                                                                                               <  2:3560092/150‑1 (MQ=255)
    cTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGtaattaat                                                                                                                                            >  2:78683/1‑150 (MQ=255)
     tCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATc                                                                                                                                           <  1:2230090/150‑1 (MQ=255)
      ccGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCg                                                                                                                                          <  2:320155/150‑1 (MQ=255)
       cGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGc                                                                                                                                         >  1:4286690/1‑150 (MQ=255)
       cGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGc                                                                                                                                         >  1:387656/1‑150 (MQ=255)
        gTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCt                                                                                                                                        <  2:661873/150‑1 (MQ=255)
           aGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTa                                                                                                                                     >  1:855743/1‑150 (MQ=255)
               tCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAAt                                                                                                                                 <  1:2875380/150‑1 (MQ=255)
                 gcgcCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATgg                                                                                                                               <  1:1674135/150‑1 (MQ=255)
                  cgcCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGt                                                                                                                              >  1:3430705/1‑150 (MQ=255)
                   gcCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTc                                                                                                                             >  1:2096889/1‑150 (MQ=255)
                    ccTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCa                                                                                                                            <  2:3215439/150‑1 (MQ=255)
                    ccTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCa                                                                                                                            >  2:2660267/1‑150 (MQ=255)
                    ccTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCa                                                                                                                            >  2:486729/1‑150 (MQ=255)
                     cTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCaa                                                                                                                           >  1:3419305/1‑150 (MQ=255)
                        gctgAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGAt                                                                                                                        <  2:4014595/147‑1 (MQ=255)
                        gCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGAt                                                                                                                        <  2:1824400/150‑1 (MQ=255)
                        gCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGAt                                                                                                                        <  2:1149852/150‑1 (MQ=255)
                        gCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGAt                                                                                                                        <  2:4318983/150‑1 (MQ=255)
                         cAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATc                                                                                                                       >  2:3648979/1‑150 (MQ=255)
                            aaTCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGAc                                                                                                                    <  2:619478/150‑1 (MQ=255)
                            aaTCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGAc                                                                                                                    >  2:59479/1‑150 (MQ=255)
                             aTCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGAcc                                                                                                                   >  1:2750427/1‑150 (MQ=255)
                             aTCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGAcc                                                                                                                   <  1:3233253/150‑1 (MQ=255)
                              tCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCt                                                                                                                  <  2:805941/150‑1 (MQ=255)
                              tCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCt                                                                                                                  <  2:3076382/150‑1 (MQ=255)
                                  acaAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAAt                                                                                                              >  1:934412/1‑150 (MQ=255)
                                   caAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAACGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATg                                                                                                             <  2:2522723/150‑1 (MQ=255)
                                    aaaaCTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGa                                                                                                            >  1:1357676/1‑150 (MQ=255)
                                      aaCTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACa                                                                                                          <  2:1962605/150‑1 (MQ=255)
                                        ctactcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTAAAGATCGACCTGAATGACAGa                                                                                                        >  2:747213/4‑150 (MQ=255)
                                             caccacCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAg                                                                                                   >  1:3610417/1‑150 (MQ=255)
                                                  ccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGt                                                                                              <  2:3450985/150‑1 (MQ=255)
                                                  ccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGt                                                                                              <  1:613068/150‑1 (MQ=255)
                                                  ccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGATAATTACTCCGGCAGCATGAGAAAGAACACGGTCGCCCTCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGt                                                                                              <  1:3409405/150‑1 (MQ=255)
                                                   cGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGtt                                                                                             <  2:89250/150‑1 (MQ=255)
                                                    gACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTg                                                                                            <  1:3590946/150‑1 (MQ=255)
                                                    gACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTg                                                                                            >  1:165659/1‑150 (MQ=255)
                                                    gACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTg                                                                                            >  2:2676828/1‑150 (MQ=255)
                                                      cccgcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGt                                                                                          >  2:3746223/1‑150 (MQ=255)
                                                       ccgcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTg                                                                                         <  2:2862416/150‑1 (MQ=255)
                                                        cgcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGc                                                                                        <  1:1163498/150‑1 (MQ=255)
                                                         gcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCg                                                                                       >  1:1689903/1‑150 (MQ=255)
                                                            cgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCt                                                                                    <  1:2059581/150‑1 (MQ=255)
                                                                aaCGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGaca                                                                                >  2:336275/1‑150 (MQ=255)
                                                                  cGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGAcaca                                                                              <  1:3331395/150‑1 (MQ=255)
                                                                      cTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAGGGATGTTGGTGCGTCTGACACAAttt                                                                          >  1:869745/1‑150 (MQ=255)
                                                                      cTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAAttt                                                                          <  1:729512/150‑1 (MQ=255)
                                                                      cTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAAttt                                                                          <  1:2818874/150‑1 (MQ=255)
                                                                        gcgcACGCAGGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTag                                                                        <  2:2407897/150‑1 (MQ=255)
                                                                          gcacgcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcg                                                                      <  2:855743/150‑1 (MQ=255)
                                                                             cgcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcgagc                                                                   >  2:3938178/1‑150 (MQ=255)
                                                                              gcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcgagcg                                                                  <  2:1688686/150‑1 (MQ=255)
                                                                                   gAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACgg                                                             >  1:1198409/1‑150 (MQ=255)
                                                                                   gAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACgg                                                             <  1:267373/150‑1 (MQ=255)
                                                                                    aaCTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACggg                                                            >  2:3320164/1‑150 (MQ=255)
                                                                                     aCTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGc                                                           >  1:496279/1‑150 (MQ=255)
                                                                                      cTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCt                                                          >  1:3602917/1‑150 (MQ=255)
                                                                                        gCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCttt                                                        <  2:2308396/150‑1 (MQ=255)
                                                                                        gCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCttt                                                        >  2:426483/1‑150 (MQ=255)
                                                                                        gCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCttt                                                        >  2:3189794/1‑150 (MQ=255)
                                                                                           aattaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTc                                                     >  1:1913717/1‑150 (MQ=255)
                                                                                               aatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAAc                                                 >  1:1875677/1‑150 (MQ=255)
                                                                                               aatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAAc                                                 >  1:3362948/1‑150 (MQ=255)
                                                                                                atCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACa                                                <  1:4170382/150‑1 (MQ=255)
                                                                                                  ccGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGt                                              <  1:1018165/150‑1 (MQ=255)
                                                                                                   cGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTg                                             >  1:3868658/1‑150 (MQ=255)
                                                                                                   cGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTg                                             <  1:4248646/150‑1 (MQ=255)
                                                                                                    ggCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTgc                                            >  2:408042/1‑150 (MQ=255)
                                                                                                     gcagcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTgcg                                           <  1:2336755/150‑1 (MQ=255)
                                                                                                         caagagaAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCAc                                       <  2:4182287/147‑1 (MQ=255)
                                                                                                          aTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACg                                      <  2:1292192/150‑1 (MQ=255)
                                                                                                          aTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACg                                      >  2:2035904/1‑150 (MQ=255)
                                                                                                            gagaAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGt                                    <  1:2442584/150‑1 (MQ=255)
                                                                                                                aagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtata                                >  2:3607334/1‑150 (MQ=255)
                                                                                                                 agaaCACGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtatat                               <  2:3617325/150‑1 (MQ=255)
                                                                                                                  gaaCACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtatata                              <  2:3892170/150‑1 (MQ=255)
                                                                                                                     cacGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAAc                           <  1:284218/150‑1 (MQ=255)
                                                                                                                     cacGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAAc                           >  1:4339157/1‑150 (MQ=255)
                                                                                                                      acGGTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAAcg                          <  2:1200143/150‑1 (MQ=255)
                                                                                                                        ggTCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAAcgcg                        <  1:1029483/150‑1 (MQ=255)
                                                                                                                          tCGCCCGCCTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGtt                      <  1:273906/150‑1 (MQ=255)
                                                                                                                            gcccgccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTaa                    <  1:2983597/150‑1 (MQ=255)
                                                                                                                             cccgccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAg                   <  2:1046897/150‑1 (MQ=255)
                                                                                                                              ccgccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAgg                  >  2:2451978/1‑150 (MQ=255)
                                                                                                                               cgccTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgt                 <  1:1031166/150‑1 (MQ=255)
                                                                                                                                gccTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgtg                >  2:3812895/1‑150 (MQ=255)
                                                                                                                                  cTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTc              <  2:3726509/150‑1 (MQ=255)
                                                                                                                                  cTGATGACCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTc              <  2:84708/150‑1 (MQ=255)
                                                                                                                                       gaCCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGcc         >  2:2621279/1‑150 (MQ=255)
                                                                                                                                        aCCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGTGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCa        >  1:1116482/1‑150 (MQ=255)
                                                                                                                                        aCCAAAGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCa        >  1:1864745/1‑150 (MQ=255)
                                                                                                                                           aaaGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGc     >  1:3898905/1‑150 (MQ=255)
                                                                                                                                           aaaGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGc     >  1:2001405/1‑150 (MQ=255)
                                                                                                                                             aGCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCt   <  2:3043741/150‑1 (MQ=255)
                                                                                                                                              gCGGTAATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTg  <  2:1481863/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTG  >  CP009273/2020019‑2020310

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: