Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,865,863 T→A 100% T166S (ACC→TCC)  dgoD ← galactonate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,865,8630TA96.8% 907.9 / 14.1 279T166S (ACC→TCC) dgoDgalactonate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (8/1);  new base A (137/133);  total (145/134)
Fisher's exact test for biased strand distribution p-value = 3.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 280 total aligned reads displayed.

GCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGCCAACCCAACTGTA  >  CP009273/3865718‑3866010
                                                                                                                                                 |                                                                                                                                                   
gCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGtt                                                                                                                                                  <  1:3842182/149‑1 (MQ=255)
   aGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAc                                                                                                                                               <  2:988326/149‑1 (MQ=255)
   aGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAc                                                                                                                                               >  2:2258053/1‑149 (MQ=255)
   aGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAc                                                                                                                                               <  2:3044031/149‑1 (MQ=255)
   aGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAc                                                                                                                                               >  1:183303/1‑149 (MQ=255)
      aCCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAccgc                                                                                                                                            >  1:1420024/1‑149 (MQ=255)
        cGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAccgccg                                                                                                                                          >  1:1016866/1‑149 (MQ=255)
        cGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAccgccg                                                                                                                                          >  1:1670607/1‑149 (MQ=255)
         ggCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACcgccgc                                                                                                                                         >  2:831996/1‑149 (MQ=255)
         ggCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACcgccgc                                                                                                                                         >  2:366920/1‑149 (MQ=255)
                cAATAACCAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTAcc                                                                                                                                  <  2:2940084/149‑1 (MQ=255)
                cAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTAcc                                                                                                                                  <  1:319802/149‑1 (MQ=255)
                   tAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACcgcg                                                                                                                               >  1:2650814/1‑149 (MQ=255)
                   tAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACcgcg                                                                                                                               >  1:673498/1‑149 (MQ=255)
                   tAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCGACcgcg                                                                                                                               >  1:2735767/1‑149 (MQ=255)
                      aCAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCCTCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGAGTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCgg                                                                                                                            <  2:3471452/149‑1 (MQ=255)
                      aCAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCgg                                                                                                                            >  2:798533/1‑149 (MQ=255)
                            ggCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTa                                                                                                                      <  1:2872017/149‑1 (MQ=255)
                             gCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTAt                                                                                                                     <  1:841468/149‑1 (MQ=255)
                              cGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATc                                                                                                                    <  2:3026589/149‑1 (MQ=255)
                                aTACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCaa                                                                                                                  <  1:1405232/149‑1 (MQ=255)
                                aTACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCaa                                                                                                                  <  1:3893396/149‑1 (MQ=255)
                                    ggCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGGTATCAATTAg                                                                                                              >  2:2329920/1‑149 (MQ=255)
                                      cTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGcc                                                                                                            >  2:3218611/1‑149 (MQ=255)
                                            ttCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGtt                                                                                                      >  2:2238232/1‑149 (MQ=255)
                                            ttCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGtt                                                                                                      <  1:2306090/149‑1 (MQ=255)
                                            ttCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTAAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGtt                                                                                                      >  1:1698674/1‑149 (MQ=255)
                                             tCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttc                                                                                                     >  2:519107/1‑149 (MQ=255)
                                             tCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttc                                                                                                     <  2:3525110/149‑1 (MQ=255)
                                             tCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttc                                                                                                     >  1:743476/1‑149 (MQ=255)
                                              cTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttct                                                                                                    <  1:2129007/149‑1 (MQ=255)
                                               tttAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttctt                                                                                                   >  2:3576052/1‑149 (MQ=255)
                                                 tAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTca                                                                                                 <  1:526936/149‑1 (MQ=255)
                                                  aaTCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTCCCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTcac                                                                                                >  1:2559891/1‑149 (MQ=255)
                                                   aTCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTcaca                                                                                               >  2:3820214/1‑149 (MQ=255)
                                                    tCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACaa                                                                                              <  2:50223/149‑1 (MQ=255)
                                                    tCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACaa                                                                                              >  1:3407341/1‑149 (MQ=255)
                                                    tCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACaa                                                                                              >  1:444613/1‑149 (MQ=255)
                                                    tCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACaa                                                                                              <  1:431823/149‑1 (MQ=255)
                                                     cAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAAc                                                                                             <  2:3520391/149‑1 (MQ=255)
                                                        cACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGt                                                                                          <  2:2563178/149‑1 (MQ=255)
                                                         aCTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGtt                                                                                         >  1:296622/1‑149 (MQ=255)
                                                          cTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTc                                                                                        <  2:2634420/149‑1 (MQ=255)
                                                           tttCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCa                                                                                       >  1:949486/1‑149 (MQ=255)
                                                             tCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGt                                                                                     <  1:1894198/149‑1 (MQ=255)
                                                             tCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGt                                                                                     <  2:119214/149‑1 (MQ=255)
                                                             tCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGt                                                                                     <  1:3051884/149‑1 (MQ=255)
                                                               gCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGttt                                                                                   >  1:2690474/1‑149 (MQ=255)
                                                                  aTCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGaa                                                                                >  1:2008229/1‑149 (MQ=255)
                                                                  aTCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGaa                                                                                <  2:1854766/149‑1 (MQ=255)
                                                                     ggCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGt                                                                             >  1:4055674/1‑149 (MQ=255)
                                                                      gcgcgcTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTa                                                                            <  2:1016866/149‑1 (MQ=255)
                                                                      gcgcgcTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTa                                                                            <  1:1095548/149‑1 (MQ=255)
                                                                        gcgcTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCAAGTTCTTCACAACCGTTCAGTTTGAAGGTATc                                                                          <  2:2771671/149‑1 (MQ=255)
                                                                          gcTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGa                                                                        <  1:2020311/149‑1 (MQ=255)
                                                                           cTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGaa                                                                       >  1:554207/1‑149 (MQ=255)
                                                                           cTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGaa                                                                       <  2:3337534/149‑1 (MQ=255)
                                                                           cTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGaa                                                                       <  1:2892378/149‑1 (MQ=255)
                                                                           cTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGaa                                                                       >  1:3728170/1‑149 (MQ=255)
                                                                           cTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGaa                                                                       <  1:341996/149‑1 (MQ=255)
                                                                             gACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGc                                                                     >  1:1544824/1‑149 (MQ=255)
                                                                             gACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGc                                                                     <  1:3457895/149‑1 (MQ=255)
                                                                                 cgACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGAt                                                                 <  2:2720797/149‑1 (MQ=255)
                                                                                  gACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGAtt                                                                <  1:3132506/149‑1 (MQ=255)
                                                                                  gACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGAtt                                                                <  1:3132118/149‑1 (MQ=255)
                                                                                  gACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGAtt                                                                >  2:2172936/1‑149 (MQ=255)
                                                                                   aCCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGAttt                                                               >  1:1753888/1‑149 (MQ=255)
                                                                                    ccGTGGAAATCACGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTc                                                              <  2:869845/149‑1 (MQ=255)
                                                                                     cGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTcg                                                             >  2:1762605/1‑149 (MQ=255)
                                                                                         gAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTCCCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTa                                                         >  2:73351/1‑149 (MQ=255)
                                                                                         gAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTa                                                         >  1:3847926/1‑149 (MQ=255)
                                                                                             tCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGt                                                     <  2:927930/149‑1 (MQ=255)
                                                                                              cAAGACCAAACTCACTCTGATTGCCAACAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGGTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGtt                                                    <  2:3952555/149‑1 (MQ=255)
                                                                                               aaGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGttt                                                   >  2:639476/1‑149 (MQ=255)
                                                                                                aGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGtttt                                                  >  2:3781866/1‑149 (MQ=255)
                                                                                                 gACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTa                                                 >  1:1155636/1‑149 (MQ=255)
                                                                                                   ccAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAAt                                               >  1:2345842/1‑149 (MQ=255)
                                                                                                   ccAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAAt                                               <  2:1698320/149‑1 (MQ=255)
                                                                                                   ccAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAAt                                               <  2:3642350/149‑1 (MQ=255)
                                                                                                    cAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATg                                              >  1:2077792/1‑149 (MQ=255)
                                                                                                     aaaaTCAATCTTATTGCCAAAAGCTTCACGATTTTGTGCCACGGAGTTAACCGCCGGATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGc                                             <  1:3784768/149‑1 (MQ=255)
                                                                                                        cTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGt                                          <  2:1607130/149‑1 (MQ=255)
                                                                                                         tCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTc                                         >  2:1490248/1‑149 (MQ=255)
                                                                                                                gATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACa                                  <  2:1793939/149‑1 (MQ=255)
                                                                                                                 aTTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGAGTTCGCGTAGCGTTTTAATGCCGTCGATAACAt                                 >  2:4083134/1‑149 (MQ=255)
                                                                                                                     ccAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATccgc                             <  2:4007718/149‑1 (MQ=255)
                                                                                                                      cAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATccgcc                            >  2:1121202/1‑149 (MQ=255)
                                                                                                                          aGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGAc                        >  1:3979407/1‑149 (MQ=255)
                                                                                                                           gCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACg                       >  2:1335314/1‑149 (MQ=255)
                                                                                                                           gCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACg                       >  1:2262198/1‑149 (MQ=255)
                                                                                                                              tCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATc                    <  2:2897634/149‑1 (MQ=255)
                                                                                                                                   aaTTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgccgc               >  1:3294719/1‑149 (MQ=255)
                                                                                                                                   aaTTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgccgc               >  1:1352/1‑149 (MQ=255)
                                                                                                                                    aTTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCgccgcc              <  1:3375030/149‑1 (MQ=255)
                                                                                                                                    aTTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCgccgcc              >  1:1133696/1‑149 (MQ=255)
                                                                                                                                      ttGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGccaa            <  2:3792124/149‑1 (MQ=255)
                                                                                                                                       tgtgCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGccaac           <  1:1252489/149‑1 (MQ=255)
                                                                                                                                         tgCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGccaaccc         <  1:1457312/149‑1 (MQ=255)
                                                                                                                                           ccACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGAAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCAccaacccaa       >  2:2426557/1‑149 (MQ=255)
                                                                                                                                                gtgtTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGCCAACCCAACTGTa  <  1:798533/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGCCAACCCAACTGTA  >  CP009273/3865718‑3866010

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: