Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 868,473 A→C 17.9% I13I (ATA→ATC yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273868,4730AC17.9% 633.8 / 165.4 291I13I (ATA→ATCyliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (121/118);  new base C (26/26);  total (147/144)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 295 total aligned reads displayed.

TAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACG  >  CP009273/868324‑868622
                                                                                                                                                     |                                                                                                                                                     
tAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATa                                                                                                                                                       >  2:1043784/1‑150 (MQ=255)
    tttCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAattt                                                                                                                                                   >  2:1359402/1‑150 (MQ=255)
      tCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCatttta                                                                                                                                                 <  1:4597448/150‑1 (MQ=255)
      tCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAatttta                                                                                                                                                 <  1:4860318/150‑1 (MQ=255)
        tgctgATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAattttatt                                                                                                                                               >  1:4721095/1‑150 (MQ=255)
        tgctgATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAattttatt                                                                                                                                               <  2:4102809/150‑1 (MQ=255)
            gATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAAttttattttta                                                                                                                                           >  1:514147/1‑150 (MQ=255)
            gATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAAttttattttta                                                                                                                                           >  1:1633576/1‑150 (MQ=255)
              ttttGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATtttatttttatt                                                                                                                                         >  1:669841/1‑150 (MQ=255)
                ttGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTgg                                                                                                                                       <  1:1937145/150‑1 (MQ=255)
                ttGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTgg                                                                                                                                       <  2:116889/150‑1 (MQ=255)
                ttGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTgg                                                                                                                                       <  2:4840744/150‑1 (MQ=255)
                 tGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGgc                                                                                                                                      >  1:2405272/1‑150 (MQ=255)
                    aGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGgcagc                                                                                                                                   <  1:1120484/150‑1 (MQ=255)
                       aTCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACt                                                                                                                                <  1:1008149/150‑1 (MQ=255)
                       aTCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACt                                                                                                                                <  1:1008321/150‑1 (MQ=255)
                        tCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTg                                                                                                                               >  1:4932148/1‑150 (MQ=255)
                         cGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGt                                                                                                                              >  1:9525/1‑150 (MQ=255)
                         cGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGt                                                                                                                              >  1:452932/1‑150 (MQ=255)
                          gggTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTc                                                                                                                             >  2:4039414/1‑150 (MQ=255)
                          gggTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTc                                                                                                                             >  1:1882867/1‑150 (MQ=255)
                           ggTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCa                                                                                                                            >  1:3489851/1‑150 (MQ=255)
                              aaTTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCAttt                                                                                                                         <  1:2258002/150‑1 (MQ=255)
                                tntGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATGTAt                                                                                                                       >  1:3555861/1‑150 (MQ=255)
                                 ttGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATTTAtt                                                                                                                      >  1:1434649/1‑150 (MQ=255)
                                     ttAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTc                                                                                                                  <  2:3777913/150‑1 (MQ=255)
                                     ttAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTc                                                                                                                  >  1:3077685/1‑150 (MQ=255)
                                      tAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCt                                                                                                                 >  1:5018715/1‑150 (MQ=255)
                                      tAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCt                                                                                                                 <  1:3367789/150‑1 (MQ=255)
                                         aaaGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCtttt                                                                                                              <  1:570926/150‑1 (MQ=255)
                                         aaaGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCtttt                                                                                                              <  1:1484549/150‑1 (MQ=255)
                                            gCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCttttttt                                                                                                           <  2:4630195/150‑1 (MQ=255)
                                              cGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTtca                                                                                                         <  1:4374609/150‑1 (MQ=255)
                                                   caAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCAtt                                                                                                    <  1:2835729/150‑1 (MQ=255)
                                                    aaGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTa                                                                                                   <  2:3611144/150‑1 (MQ=255)
                                                    aaGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATATTTATTGGAAGAACTGTCATTTATTGTTCTTTTTTTCATCATTa                                                                                                   >  1:3444827/1‑150 (MQ=255)
                                                        gACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTa                                                                                               <  2:1446979/150‑1 (MQ=255)
                                                         aCAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTAt                                                                                              <  2:3706283/150‑1 (MQ=255)
                                                            aTTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCaa                                                                                           >  2:3252182/1‑150 (MQ=255)
                                                              ttttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACt                                                                                         >  2:1685660/1‑150 (MQ=255)
                                                               tttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTg                                                                                        <  1:2583162/150‑1 (MQ=255)
                                                                ttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTgg                                                                                       <  2:3357166/150‑1 (MQ=255)
                                                                   cAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTAt                                                                                    <  2:2565613/150‑1 (MQ=255)
                                                                    aGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATc                                                                                   <  1:4280240/150‑1 (MQ=255)
                                                                     gCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCg                                                                                  <  2:4730421/150‑1 (MQ=255)
                                                                           caGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAgcgc                                                                            >  2:341454/1‑150 (MQ=255)
                                                                             gCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAgcgcgc                                                                          <  2:669841/150‑1 (MQ=255)
                                                                                   aGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTaaa                                                                    >  1:3857995/1‑150 (MQ=255)
                                                                                   aGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTaaa                                                                    <  2:4916815/150‑1 (MQ=255)
                                                                                    gATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTaaaa                                                                   <  2:4292882/150‑1 (MQ=255)
                                                                                     aTAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAAc                                                                  >  2:2737669/1‑150 (MQ=255)
                                                                                       aGTTATTTTGTTAAATTTATTAACATGCTGAGTTTATACTAAAATATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGTAGCGCGCGGTAAAACGt                                                                <  2:3192588/150‑1 (MQ=255)
                                                                                       aGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGt                                                                >  1:1058040/1‑150 (MQ=255)
                                                                                        gTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTg                                                               <  2:153904/150‑1 (MQ=255)
                                                                                        gTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTg                                                               >  2:2381150/1‑150 (MQ=255)
                                                                                          tATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGcc                                                             <  2:1995728/150‑1 (MQ=255)
                                                                                             tttGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGAt                                                          <  1:397298/150‑1 (MQ=255)
                                                                                              ttGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATa                                                         <  2:1927547/150‑1 (MQ=255)
                                                                                              ttGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATAAAATGGTATCGGAGCGCGAGGTAAAACGTGCCGATa                                                         >  2:3996193/1‑150 (MQ=255)
                                                                                               tGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTATCTATCAACTGGTATCGGAGCGGGCGGTAATACGTGCCGATAg                                                        >  2:3609113/1‑150 (MQ=255)
                                                                                                gTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGc                                                       <  1:2025154/150‑1 (MQ=255)
                                                                                                    aaTGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCt                                                   >  2:4587552/1‑150 (MQ=255)
                                                                                                      tGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCttt                                                 <  2:8413/150‑1 (MQ=255)
                                                                                                       gTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTg                                                <  1:2439851/150‑1 (MQ=255)
                                                                                                           tAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACt                                            >  2:2140255/1‑150 (MQ=255)
                                                                                                           tAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTATTTTTTTCATCATTAACTATCAACTGGTATCGGAGAGCGCGGTAAAACTTGCCGATAGCCGCTTTGAACt                                            >  1:4358719/1‑150 (MQ=255)
                                                                                                             aCATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTa                                          >  1:1374932/1‑150 (MQ=255)
                                                                                                             aCATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTa                                          >  2:2604408/1‑150 (MQ=255)
                                                                                                             aCATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTa                                          >  2:4459307/1‑150 (MQ=255)
                                                                                                              cATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTAt                                         >  1:3131994/1‑150 (MQ=255)
                                                                                                               aTGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTAtt                                        >  1:4147406/1‑150 (MQ=255)
                                                                                                               aTGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTAGTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTAtt                                        >  1:2231038/1‑150 (MQ=255)
                                                                                                                 gCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCa                                      <  2:3857995/150‑1 (MQ=255)
                                                                                                                     aGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGaaa                                  <  1:3337532/150‑1 (MQ=255)
                                                                                                                     aGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGaaa                                  <  2:3663251/150‑1 (MQ=255)
                                                                                                                        ttATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAAc                               <  1:1748135/150‑1 (MQ=255)
                                                                                                                          ataCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGt                             >  1:3707580/1‑150 (MQ=255)
                                                                                                                          ataCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATAAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGt                             >  1:4083299/1‑150 (MQ=255)
                                                                                                                           taCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGtt                            <  1:821697/150‑1 (MQ=255)
                                                                                                                            aCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTg                           <  1:1715132/150‑1 (MQ=255)
                                                                                                                            aCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTg                           <  1:1715081/150‑1 (MQ=255)
                                                                                                                            aCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTg                           <  2:1230867/150‑1 (MQ=255)
                                                                                                                             cGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTgg                          <  1:4007964/150‑1 (MQ=255)
                                                                                                                                 aagataaagataAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTATTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAt                      <  1:2568024/150‑1 (MQ=255)
                                                                                                                                    ataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTc                   <  2:728622/150‑1 (MQ=255)
                                                                                                                                    ataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTc                   <  2:4932157/150‑1 (MQ=255)
                                                                                                                                     taaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCt                  >  2:3782835/1‑150 (MQ=255)
                                                                                                                                       aagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCttt                <  2:3303717/150‑1 (MQ=255)
                                                                                                                                        agataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTTTTTa               <  2:514147/150‑1 (MQ=255)
                                                                                                                                        agataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTa               >  2:3562557/1‑150 (MQ=255)
                                                                                                                                        agataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTa               <  2:1882867/150‑1 (MQ=255)
                                                                                                                                        agataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACATATTCAGAAAAACGTTGGCTATTTCTTTa               >  1:1962043/1‑150 (MQ=255)
                                                                                                                                            aaGGCTGATCATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGa           >  2:4133118/1‑150 (MQ=255)
                                                                                                                                            aaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGa           <  2:2051371/150‑1 (MQ=255)
                                                                                                                                            aaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGa           <  2:2482640/150‑1 (MQ=255)
                                                                                                                                                cTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATAtt       <  2:4958751/150‑1 (MQ=255)
                                                                                                                                                 tGATCATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTg      <  1:4410115/150‑1 (MQ=255)
                                                                                                                                                 tGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTg      >  1:2710513/1‑150 (MQ=255)
                                                                                                                                                 tGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTg      >  2:1620757/1‑150 (MQ=255)
                                                                                                                                                  gATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGa     >  1:1271662/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAACTTTCTGCTGATTTTGTAGAATCGGGTAATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACG  >  CP009273/868324‑868622

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: