Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,621,733 T→A 13.2% intergenic (‑96/‑41) dcp ← / → ydfG dipeptidyl carboxypeptidase II/NADP‑dependent 3‑hydroxy acid dehydrogenase; malonic semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,621,7330TA13.2% 730.1 / 84.5 257intergenic (‑96/‑41)dcp/ydfGdipeptidyl carboxypeptidase II/NADP‑dependent 3‑hydroxy acid dehydrogenase; malonic semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base T (105/118);  new base A (26/8);  total (131/126)
Fisher's exact test for biased strand distribution p-value = 1.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 262 total aligned reads displayed.

GGAGCCAGATACGGCAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGG  >  CP009273/1621585‑1621881
                                                                                                                                                    |                                                                                                                                                    
ggAGCCAGATACGGCAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGt                                                                                                                                                      >  2:1699376/1‑149 (MQ=255)
 gAGCCAGATACGGCAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCCTCCCGTg                                                                                                                                                     >  2:2536662/1‑149 (MQ=255)
  aGCCAGATACGGCAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGc                                                                                                                                                    <  1:391206/149‑1 (MQ=255)
  aGCCAGATACGGCAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGc                                                                                                                                                    <  1:3704420/149‑1 (MQ=255)
   gCCAGATACGGCAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCt                                                                                                                                                   <  2:1172866/149‑1 (MQ=255)
   gCCAGATACGGCAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCt                                                                                                                                                   <  2:390436/149‑1 (MQ=255)
     cAGATACGGCAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGAGCTAt                                                                                                                                                 >  2:308569/1‑149 (MQ=255)
       gATACGGCAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGt                                                                                                                                               >  2:3145141/1‑149 (MQ=255)
        aTACGGCAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGtt                                                                                                                                              >  1:510800/1‑149 (MQ=255)
        aTACGGCAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGtt                                                                                                                                              <  1:3775103/149‑1 (MQ=255)
              cAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGac                                                                                                                                        >  1:4000418/1‑149 (MQ=255)
                  gtgCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGacacac                                                                                                                                    >  2:2577468/1‑149 (MQ=255)
                   tgCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGAcacaca                                                                                                                                   >  2:556802/1‑149 (MQ=255)
                      tttGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACaaaa                                                                                                                                >  2:865654/1‑149 (MQ=255)
                      tttGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACaaaa                                                                                                                                >  2:659673/1‑149 (MQ=255)
                      tttGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACaaaa                                                                                                                                >  2:1888074/1‑149 (MQ=255)
                          cacaAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGt                                                                                                                            <  2:2219743/149‑1 (MQ=255)
                           acaAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGtt                                                                                                                           <  2:1224826/149‑1 (MQ=255)
                             aaggaaaggaTTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGa                                                                                                                         <  1:177949/149‑1 (MQ=255)
                                 aaaggaTTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGaa                                                                                                                     >  1:868864/1‑149 (MQ=255)
                                 aaaggaTTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGaa                                                                                                                     <  1:963034/149‑1 (MQ=255)
                                  aaggaTTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAAc                                                                                                                    <  2:2368514/149‑1 (MQ=255)
                                  aaggaTTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGAGCTATGTTATTGACACACAAAAGCGTTGAGGAAc                                                                                                                    >  1:2723287/1‑149 (MQ=255)
                                    ggaTTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAg                                                                                                                  >  2:4030350/1‑149 (MQ=255)
                                      aTTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTg                                                                                                                <  2:3248720/149‑1 (MQ=255)
                                       ttCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTTACACACAAAAGCGTTGAGTAACAGTga                                                                                                               >  2:1571308/1‑149 (MQ=255)
                                       ttCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTga                                                                                                               <  1:2772854/149‑1 (MQ=255)
                                       ttCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGAGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTga                                                                                                               <  2:1549792/149‑1 (MQ=255)
                                       ttCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGAGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTga                                                                                                               >  2:3882622/1‑149 (MQ=255)
                                         cATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTgaga                                                                                                             <  2:3358475/149‑1 (MQ=255)
                                           ttgttgTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAATCGTTGAGGAACAGTGAgatg                                                                                                           >  1:1885409/1‑149 (MQ=255)
                                           ttgttgTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAgatg                                                                                                           <  1:2292114/149‑1 (MQ=255)
                                            tgttgtCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAgatga                                                                                                          >  2:2001874/1‑149 (MQ=255)
                                               tgtCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGAGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCg                                                                                                       <  2:208027/149‑1 (MQ=255)
                                                  cATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGttt                                                                                                    <  1:304084/149‑1 (MQ=255)
                                                  cATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGAGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGttt                                                                                                    >  1:2026574/1‑149 (MQ=255)
                                                         ttctCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAAc                                                                                             <  1:3217410/149‑1 (MQ=255)
                                                         ttctCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAAc                                                                                             >  2:2996437/1‑149 (MQ=255)
                                                          tctcCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACt                                                                                            <  2:3991212/149‑1 (MQ=255)
                                                               tGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGc                                                                                       <  2:2555240/149‑1 (MQ=255)
                                                               tGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGc                                                                                       >  2:699053/1‑149 (MQ=255)
                                                                gAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCa                                                                                      <  2:741610/149‑1 (MQ=255)
                                                                  aCGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAAc                                                                                    >  2:852464/1‑149 (MQ=255)
                                                                  aCGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAAc                                                                                    <  2:1774881/149‑1 (MQ=255)
                                                                      gAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGAGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCa                                                                                >  1:70727/1‑149 (MQ=255)
                                                                       aGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCNCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAg                                                                               >  1:2087697/1‑149 (MQ=255)
                                                                       aGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAg                                                                               >  2:689849/1‑149 (MQ=255)
                                                                       aGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGGTGAGGAACAGGGAGATGATCGTTGTAGTACGTGGAGCTACGGccg                                                                               >  2:3481071/1‑147 (MQ=255)
                                                                          tgtTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGtt                                                                            <  1:1037793/149‑1 (MQ=255)
                                                                            ttCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGtttt                                                                          <  2:2614082/149‑1 (MQ=255)
                                                                               cATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGt                                                                       >  2:775380/1‑149 (MQ=255)
                                                                               cATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGATCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGt                                                                       <  1:2631080/149‑1 (MQ=255)
                                                                                  aGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGaa                                                                    <  1:1863794/149‑1 (MQ=255)
                                                                                      tAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCa                                                                >  1:2424650/1‑149 (MQ=255)
                                                                                        ggCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCAtt                                                              <  2:2468087/149‑1 (MQ=255)
                                                                                        ggCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCAtt                                                              >  1:1513372/1‑149 (MQ=255)
                                                                                          cTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTAc                                                            <  2:491128/149‑1 (MQ=255)
                                                                                          cTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTAc                                                            <  1:1567471/149‑1 (MQ=255)
                                                                                              cTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACtcgt                                                        >  1:1686886/1‑149 (MQ=255)
                                                                                                 aTAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGAGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTcgtcgt                                                     >  2:2307882/1‑149 (MQ=255)
                                                                                                    aCGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGtttt                                                  <  2:1471746/149‑1 (MQ=255)
                                                                                                    aCGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGtttt                                                  <  2:2770108/149‑1 (MQ=255)
                                                                                                    aCGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGATATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGtttt                                                  <  2:1978092/149‑1 (MQ=255)
                                                                                                     cGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTa                                                 <  2:2144922/149‑1 (MQ=255)
                                                                                                       gAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTAtt                                               >  1:2342134/1‑149 (MQ=255)
                                                                                                        aGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTc                                              <  2:3333829/149‑1 (MQ=255)
                                                                                                              aTCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAg                                        >  1:1170580/1‑149 (MQ=255)
                                                                                                                cTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAggg                                      >  2:3130859/1‑149 (MQ=255)
                                                                                                                 ttGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGc                                     <  1:616069/149‑1 (MQ=255)
                                                                                                                 ttGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGc                                     <  1:2894625/149‑1 (MQ=255)
                                                                                                                 ttGCAATCCAGTATTACCCGCTCTTAAGCATCCCGAGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGc                                     >  2:402707/1‑149 (MQ=255)
                                                                                                                 ttGCAATCCAGTAATACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTTTTAAAAAAGGGc                                     >  1:49789/1‑149 (MQ=255)
                                                                                                                  tGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCa                                    <  1:184484/149‑1 (MQ=255)
                                                                                                                    cAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATa                                  <  1:98155/149‑1 (MQ=255)
                                                                                                                      aTCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATaaa                                >  1:1621226/1‑149 (MQ=255)
                                                                                                                       tCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAg                               <  1:2999151/149‑1 (MQ=255)
                                                                                                                       tCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAg                               <  1:3420477/149‑1 (MQ=255)
                                                                                                                          aGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTa                            <  1:1111519/149‑1 (MQ=255)
                                                                                                                          aGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTa                            >  2:2428097/1‑149 (MQ=255)
                                                                                                                          aGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTa                            <  1:1332163/149‑1 (MQ=255)
                                                                                                                           gTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTAt                           >  1:355399/1‑149 (MQ=255)
                                                                                                                            tATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATc                          <  1:2492194/149‑1 (MQ=255)
                                                                                                                             aTTACCCGCTCTTAAGCATCCCGAGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAATGGCATAAAGTTAtag                         >  1:777565/1‑147 (MQ=255)
                                                                                                                                aCCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCa                      >  2:3082322/1‑149 (MQ=255)
                                                                                                                                 cccGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCAc                     <  2:3728459/149‑1 (MQ=255)
                                                                                                                                 cccGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCAc                     >  2:2190965/1‑149 (MQ=255)
                                                                                                                                  ccGCTCTTAAGCATCCCGAGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACt                    >  2:1862860/1‑149 (MQ=255)
                                                                                                                                   cGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTg                   >  1:3069070/1‑149 (MQ=255)
                                                                                                                                      tctTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTGTTATTGAAAAAGGGCATAAAGTTATCGCCACTGGcc                >  1:30335/1‑149 (MQ=255)
                                                                                                                                        ttAAGCATCCCGAGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTAATCAACAAGGGCATAAAGTTATCGCCACTGGccct              >  2:1463417/1‑147 (MQ=255)
                                                                                                                                          aaGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCg            >  1:2064014/1‑149 (MQ=255)
                                                                                                                                          aaGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCg            >  1:3238247/1‑149 (MQ=255)
                                                                                                                                           aGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGc           >  1:1711891/1‑149 (MQ=255)
                                                                                                                                            gCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGcc          <  1:556802/149‑1 (MQ=255)
                                                                                                                                              aTCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAg        >  1:3355490/1‑149 (MQ=255)
                                                                                                                                              aTCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAg        <  1:3296188/149‑1 (MQ=255)
                                                                                                                                               tCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGACCAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAgg       >  2:964885/1‑149 (MQ=255)
                                                                                                                                                cccGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGa      <  1:4030350/149‑1 (MQ=255)
                                                                                                                                                  cGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAAc    >  1:745601/1‑149 (MQ=255)
                                                                                                                                                   gTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACg   >  2:2396318/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGAGCCAGATACGGCAGTGTGCTTTGCACAAGGAAAGGATTCATTGTTGTCATTTTCTTCTCCTGAACGCGAGGTGTTCCATAGCGTAGGCTTACTGATAACGGAGTGCAATCTTGCAATCCAGTATTACCCGCTCTTAAGCATCCCGTGCTATGTTATTGACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGG  >  CP009273/1621585‑1621881

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: