Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,139,450 G→T 17.8% E902* (GAG→TAG)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,139,4500GT17.8% 788.8 / 116.0 270E902* (GAG→TAG) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (117/105);  new base T (18/30);  total (135/135)
Fisher's exact test for biased strand distribution p-value = 7.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

Only 100 of 283 total aligned reads displayed.

GAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTT  >  CP009273/2139303‑2139599
                                                                                                                                                   |                                                                                                                                                     
gAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCtag                                                                                                                                                     >  1:3811233/1‑150 (MQ=255)
gAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCtag                                                                                                                                                     >  1:1032139/1‑150 (MQ=255)
gAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCgag                                                                                                                                                     >  2:2260392/1‑150 (MQ=255)
gAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCgag                                                                                                                                                     >  2:2652842/1‑150 (MQ=255)
 aaCGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCgaga                                                                                                                                                    <  1:3957108/150‑1 (MQ=255)
   cGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGatt                                                                                                                                                  <  2:2113317/150‑1 (MQ=255)
    gggCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGatta                                                                                                                                                 >  2:3958491/1‑150 (MQ=255)
    gggCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGatta                                                                                                                                                 <  2:3131767/150‑1 (MQ=255)
      gcggcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGattatt                                                                                                                                               <  1:2641214/150‑1 (MQ=255)
         gagATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGTTGATGCTAGCCCCAGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGTGATTTTtgat                                                                                                                                            <  1:2055535/148‑1 (MQ=255)
          cgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATtgatg                                                                                                                                           <  1:534311/150‑1 (MQ=255)
           gATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtgatga                                                                                                                                          >  1:4130508/1‑150 (MQ=255)
           gATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtgatga                                                                                                                                          <  2:2279126/150‑1 (MQ=255)
             tGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaac                                                                                                                                        <  1:307258/150‑1 (MQ=255)
              gTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaaca                                                                                                                                       <  1:3687082/150‑1 (MQ=255)
              gTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaaca                                                                                                                                       <  1:3343516/150‑1 (MQ=255)
               tCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaacaa                                                                                                                                      <  1:137915/150‑1 (MQ=255)
                 gCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaacaaac                                                                                                                                    <  1:468129/150‑1 (MQ=255)
                  cTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaacaaaca                                                                                                                                   >  2:3747304/1‑150 (MQ=255)
                    tgatgaACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAtt                                                                                                                                 >  1:886564/1‑150 (MQ=255)
                    tgatgaACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAtt                                                                                                                                 <  1:1706165/150‑1 (MQ=255)
                     gatgaACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAttt                                                                                                                                >  2:1527335/1‑150 (MQ=255)
                      attaaCAGTGGCGTATGATTAAAGAGATTCAGTTGATGATTCTCTCCCACGGTGTCGCTTCTCAACGGATCTCGTAATAGCGTAATTTGTGTCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATTTACAATCAtttg                                                                                                                               <  2:2447279/147‑2 (MQ=255)
                        gaACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTtcgt                                                                                                                             >  2:1790630/1‑150 (MQ=255)
                          aCAGTGGCGGTTGATTTAAGTGTATCAGTTGAGGATGCTCGCACAGGGTGTCGCTTCGCCACGGATCCCGGAAGCTAGTAATTTGTGTCTGATTTGACTTAAGCTCTGTAGTTGCGAAGGCGAGATTATTGATGTACTAACATTtcgtcg                                                                                                                           <  2:538067/150‑1 (MQ=255)
                            aGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTa                                                                                                                         <  2:1576390/150‑1 (MQ=255)
                             gTGGCGGAGGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGTAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAg                                                                                                                        <  2:3467445/150‑1 (MQ=255)
                               ggCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCt                                                                                                                      >  1:2230104/1‑150 (MQ=255)
                                gCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCtt                                                                                                                     <  2:1255481/150‑1 (MQ=255)
                                 cGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTc                                                                                                                    >  2:1325632/1‑150 (MQ=255)
                                 cGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTc                                                                                                                    >  2:234160/1‑150 (MQ=255)
                                    atgatTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGc                                                                                                                 <  2:3268889/150‑1 (MQ=255)
                                     tgatTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCg                                                                                                                <  2:1901020/150‑1 (MQ=255)
                                         tAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATcc                                                                                                            <  1:484880/150‑1 (MQ=255)
                                             gagaATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGgcg                                                                                                        <  1:3952336/150‑1 (MQ=255)
                                             gagaATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTACACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGgcg                                                                                                        >  1:1074636/1‑150 (MQ=255)
                                               gaATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCt                                                                                                      <  1:2975164/150‑1 (MQ=255)
                                                    aGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGcc                                                                                                 >  1:2897460/1‑150 (MQ=255)
                                                     gTTGATGATGCTCGCCCTCGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCTGCGATCCGGCGCTTAGCCa                                                                                                <  2:1671758/150‑1 (MQ=255)
                                                      ttGCTGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCAt                                                                                               <  2:3802763/150‑1 (MQ=255)
                                                       tgatgaTGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATg                                                                                              <  2:1095555/150‑1 (MQ=255)
                                                       tgatgaTGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATg                                                                                              <  1:1845683/150‑1 (MQ=255)
                                                        gatgatGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGc                                                                                             >  2:9225/1‑150 (MQ=255)
                                                            atgCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCtt                                                                                         >  2:2600626/1‑150 (MQ=255)
                                                             tgCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTg                                                                                        <  1:3412954/150‑1 (MQ=255)
                                                                   cccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGAccgcc                                                                                  >  2:4134066/1‑150 (MQ=255)
                                                                    ccACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGAccgccg                                                                                 >  2:364887/1‑150 (MQ=255)
                                                                      aCGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCggg                                                                               >  1:1291488/1‑150 (MQ=255)
                                                                        ggTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTa                                                                             <  1:3917072/150‑1 (MQ=255)
                                                                         gtgtCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTAt                                                                            <  1:3274883/150‑1 (MQ=255)
                                                                           gtCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTc                                                                          <  2:418531/150‑1 (MQ=255)
                                                                            tCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTcc                                                                         <  2:1036095/150‑1 (MQ=255)
                                                                             cGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCa                                                                        <  2:204570/150‑1 (MQ=255)
                                                                                ttCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGa                                                                     >  2:70922/1‑150 (MQ=255)
                                                                                        ggATCCCGTAAGCGCGTAATTTCTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATATAGACGAAGCTTCAGCGATCCGGCGCTTAGCCAACCTCTAGACCTCCGGCTCTTCAACGAATTTTAcc                                                             >  1:1752347/1‑150 (MQ=255)
                                                                                          aTCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcag                                                           <  1:861769/150‑1 (MQ=255)
                                                                                                  aGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGc                                                   >  1:2940615/1‑150 (MQ=255)
                                                                                                   gcgcGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCa                                                  >  2:1123184/1‑150 (MQ=255)
                                                                                                   gcgcGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCa                                                  >  2:281213/1‑150 (MQ=255)
                                                                                                       gTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAg                                              <  2:3430869/150‑1 (MQ=255)
                                                                                                        tAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGc                                             >  2:1386328/1‑150 (MQ=255)
                                                                                                        tAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGc                                             >  1:212322/1‑150 (MQ=255)
                                                                                                           tttGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGt                                          >  1:3941305/1‑150 (MQ=255)
                                                                                                              gtgGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGc                                       <  1:2674785/150‑1 (MQ=255)
                                                                                                               tgGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGcc                                      >  1:704358/1‑150 (MQ=255)
                                                                                                                ggCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCa                                     >  2:3091619/1‑150 (MQ=255)
                                                                                                                 gCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAg                                    >  1:1585480/1‑150 (MQ=255)
                                                                                                                  cTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGt                                   <  2:2538451/150‑1 (MQ=255)
                                                                                                                    gATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTaa                                 >  1:2283666/1‑150 (MQ=255)
                                                                                                                          aCTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCAt                           >  2:1937811/1‑150 (MQ=255)
                                                                                                                           cTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATa                          <  1:3569633/150‑1 (MQ=255)
                                                                                                                           cTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATa                          >  1:62601/1‑150 (MQ=255)
                                                                                                                           cTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATa                          >  2:2887307/1‑150 (MQ=255)
                                                                                                                              aaGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGc                       <  2:1734871/150‑1 (MQ=255)
                                                                                                                               aGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCa                      >  1:1202610/1‑150 (MQ=255)
                                                                                                                                gCTCTGGAGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAt                     <  1:2260626/150‑1 (MQ=255)
                                                                                                                                  tctGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAtct                   <  1:1992512/150‑1 (MQ=255)
                                                                                                                                  tctGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAtct                   >  1:3454272/1‑150 (MQ=255)
                                                                                                                                  tctGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAtct                   >  1:1959441/1‑150 (MQ=255)
                                                                                                                                   ctGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAtctc                  >  1:1979841/1‑150 (MQ=255)
                                                                                                                                       aGTTGCGATGGCGTGATTATTGATGAACAAACATTGAGTCGGAGCTTCAGCGATCCGGCGCTTAGACATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCTCAAAAGCGTTTGCAAGTAAAGGCATAAGCATatccctc              <  2:662138/150‑1 (MQ=255)
                                                                                                                                       aGTTGCGAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATctccctc              <  2:864670/150‑1 (MQ=255)
                                                                                                                                         ttGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCtccctccc            <  1:2799551/150‑1 (MQ=255)
                                                                                                                                         ttGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGCTGCCAGTAAAGGCATAAGCATCtccctccc            <  1:3036104/150‑1 (MQ=255)
                                                                                                                                          tGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTcccc           >  2:2938618/1‑150 (MQ=255)
                                                                                                                                          tGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTcccc           >  1:586663/1‑150 (MQ=255)
                                                                                                                                           gCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTccccc          <  1:1926152/150‑1 (MQ=255)
                                                                                                                                            cGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCt         >  2:3183029/1‑150 (MQ=255)
                                                                                                                                            cGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCt         >  1:2444242/1‑150 (MQ=255)
                                                                                                                                            cGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCt         >  2:2706919/1‑150 (MQ=255)
                                                                                                                                            cGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCt         >  2:1764658/1‑150 (MQ=255)
                                                                                                                                            cGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCt         >  2:1945734/1‑150 (MQ=255)
                                                                                                                                             gAAGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCtt        <  1:1475822/150‑1 (MQ=255)
                                                                                                                                             gAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCtt        >  2:1439123/1‑150 (MQ=255)
                                                                                                                                             gAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCtt        >  1:4059255/1‑150 (MQ=255)
                                                                                                                                               aGGCTAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTc      >  2:4165730/1‑150 (MQ=255)
                                                                                                                                                ggCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTcc     >  2:1195438/1‑150 (MQ=255)
                                                                                                                                                ggCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTcc     >  2:4042616/1‑150 (MQ=255)
                                                                                                                                                  cGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGt   >  2:720159/1‑150 (MQ=255)
                                                                                                                                                   gagaTTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGtt  <  2:3347120/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCCCTCCCCCTTTCCGTT  >  CP009273/2139303‑2139599

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: