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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,865,8630TA85.0% 925.7 / 136.8 326T166S (ACC→TCC) dgoDgalactonate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (25/24);  new base A (119/158);  total (144/182)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 330 total aligned reads displayed.

CGCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGCCAACCCAACTG  >  CP009273/3865717‑3866008
                                                                                                                                                  |                                                                                                                                                 
cGCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGgagt                                                                                                                                                 >  2:4241953/1‑149 (MQ=255)
  ccAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTa                                                                                                                                               >  2:4069333/1‑149 (MQ=255)
   cAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTaa                                                                                                                                              >  1:2543905/1‑149 (MQ=255)
   cAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTaa                                                                                                                                              <  2:4154588/149‑1 (MQ=255)
    aGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAc                                                                                                                                             <  1:4711115/149‑1 (MQ=255)
    aGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAc                                                                                                                                             >  1:1148465/1‑149 (MQ=255)
        ccGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAccgcc                                                                                                                                         <  1:593445/149‑1 (MQ=255)
        ccGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAccgcc                                                                                                                                         >  2:2265980/1‑149 (MQ=255)
              cctcAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCt                                                                                                                                   <  2:3219530/149‑1 (MQ=255)
              cctcAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCt                                                                                                                                   <  1:1658963/149‑1 (MQ=255)
              cctcAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCt                                                                                                                                   <  2:3721863/149‑1 (MQ=255)
              cctcAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCGGATTGCCAAAAGCCTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCt                                                                                                                                   <  2:3411605/149‑1 (MQ=255)
               ctcAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTa                                                                                                                                  >  1:3576891/1‑149 (MQ=255)
                 cAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTAcc                                                                                                                                <  2:2035182/149‑1 (MQ=255)
                     aaaCAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACcgcgc                                                                                                                            <  1:630919/149‑1 (MQ=255)
                        cAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCGAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCggg                                                                                                                         >  1:3638209/1‑149 (MQ=255)
                        cAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCGAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCggg                                                                                                                         >  1:1248029/1‑149 (MQ=255)
                           cGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGt                                                                                                                      <  2:256446/149‑1 (MQ=255)
                             ggCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTa                                                                                                                    <  2:731536/149‑1 (MQ=255)
                             ggCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTa                                                                                                                    <  2:1740443/149‑1 (MQ=255)
                             ggCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTa                                                                                                                    <  2:1423693/149‑1 (MQ=255)
                              gCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTAt                                                                                                                   <  1:1789777/149‑1 (MQ=255)
                               cGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTATCCGCCGCATCTACCGCGCGGGAGTTATc                                                                                                                  <  1:2974575/149‑1 (MQ=255)
                                 aTACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCaa                                                                                                                >  1:29402/1‑149 (MQ=255)
                                 aTACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCaa                                                                                                                >  1:1306283/1‑149 (MQ=255)
                                 aTACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCaa                                                                                                                >  1:3695266/1‑149 (MQ=255)
                                   aCGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAAtt                                                                                                              <  1:1918047/149‑1 (MQ=255)
                                    cGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTa                                                                                                             <  1:979517/149‑1 (MQ=255)
                                     ggCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTCCCGCGCGGGAGTTATCAATTAg                                                                                                            >  1:429660/1‑149 (MQ=255)
                                      tcTCCAGTTCTTTAATCAGCACTTTCGCCATGGGCGCGCTGACGCGACCGTGGAAATCAAGTCCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGc                                                                                                           <  2:2404917/148‑1 (MQ=255)
                                        tCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGccc                                                                                                         >  1:2233249/1‑149 (MQ=255)
                                        tCCAGGTCTTTACTCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGccc                                                                                                         <  2:2970098/149‑1 (MQ=255)
                                          cAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCa                                                                                                       <  1:292298/149‑1 (MQ=255)
                                          cAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCa                                                                                                       >  2:310843/1‑149 (MQ=255)
                                            gTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGt                                                                                                     >  1:3457814/1‑149 (MQ=255)
                                            gTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGt                                                                                                     >  2:1964109/1‑149 (MQ=255)
                                               cTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttct                                                                                                  >  1:158571/1‑149 (MQ=255)
                                               cTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttct                                                                                                  <  2:2744958/149‑1 (MQ=255)
                                               cTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttct                                                                                                  >  1:379974/1‑149 (MQ=255)
                                               cTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttct                                                                                                  <  2:579541/149‑1 (MQ=255)
                                                tttAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttctt                                                                                                 <  2:4484739/149‑1 (MQ=255)
                                                 ttAATCAGCAGTGTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCCAGACCAACCTCAATCTGATTGCCTACAGCTTCACGAATTTGTGCCACGGTGTTAAGCGCCTCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttcttc                                                                                                <  1:4191043/149‑1 (MQ=255)
                                                 ttAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCTCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttcttc                                                                                                >  1:1984920/1‑149 (MQ=255)
                                                    aTCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTcaca                                                                                             >  1:1302540/1‑149 (MQ=255)
                                                       aGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAAcc                                                                                          <  1:824110/149‑1 (MQ=255)
                                                        gCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCg                                                                                         <  2:293172/149‑1 (MQ=255)
                                                         cACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGt                                                                                        <  1:3028980/149‑1 (MQ=255)
                                                             ttCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAg                                                                                    <  1:2797299/149‑1 (MQ=255)
                                                             ttCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAg                                                                                    <  2:2366859/149‑1 (MQ=255)
                                                             ttCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAg                                                                                    <  2:4708801/149‑1 (MQ=255)
                                                                gCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGttt                                                                                 >  1:3154205/1‑149 (MQ=255)
                                                                 ccATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTg                                                                                <  2:3046658/149‑1 (MQ=255)
                                                                 ccATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTg                                                                                >  2:2106901/1‑149 (MQ=255)
                                                                  cATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGa                                                                               >  1:804064/1‑149 (MQ=255)
                                                                    tCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAg                                                                             >  1:1414048/1‑149 (MQ=255)
                                                                    tCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAg                                                                             <  2:3147224/149‑1 (MQ=255)
                                                                        cgcgcTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTAt                                                                         <  1:2185597/149‑1 (MQ=255)
                                                                            cTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGaa                                                                     <  2:4593832/149‑1 (MQ=255)
                                                                            cTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGaa                                                                     <  1:1129173/149‑1 (MQ=255)
                                                                                   gACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGAtt                                                              >  1:747654/1‑149 (MQ=255)
                                                                                    aCCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGAttt                                                             <  1:2601679/149‑1 (MQ=255)
                                                                                     ccGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTc                                                            <  1:4773062/149‑1 (MQ=255)
                                                                                        tGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTcgcg                                                         >  1:4415094/1‑149 (MQ=255)
                                                                                        tGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTcgcg                                                         >  2:1312000/1‑149 (MQ=255)
                                                                                         ggAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGt                                                        <  1:1848977/149‑1 (MQ=255)
                                                                                         ggAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGt                                                        <  2:2373858/149‑1 (MQ=255)
                                                                                             aTCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCg                                                    >  2:2697439/1‑149 (MQ=255)
                                                                                               cAAGACCTAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGtt                                                  <  1:310843/149‑1 (MQ=255)
                                                                                                aaGACCACACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGttt                                                 <  1:1885268/149‑1 (MQ=255)
                                                                                                aaGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGttt                                                 >  2:24839/1‑149 (MQ=255)
                                                                                                 aGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGtttt                                                <  2:2615820/149‑1 (MQ=255)
                                                                                                   acaaaaCTCAATCTGATTGCCAAAAGCTTCAAGAATATGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTACACAACCGATCAGTTTGAAGGAATCGAAGCCGATTTCGCGTAGCGTATTaa                                              >  2:4779553/4‑149 (MQ=255)
                                                                                                    ccAAACTCACTCTGATTGACAAAAGCTTCACGAATTTGTGCCTCGGAGTTAACCGCCGCATCTACCGCGAGGGAGTTATCAATTAGCCCCAGTTCTTCTCATCCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAAt                                             <  2:2666592/149‑1 (MQ=255)
                                                                                                    ccAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAAt                                             <  1:3090616/149‑1 (MQ=255)
                                                                                                     cAAACTCACTCTGATTGCCAACAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATg                                            <  2:3628594/149‑1 (MQ=255)
                                                                                                     cAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCTCAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATg                                            >  2:2943515/1‑149 (MQ=255)
                                                                                                       aaCTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGcc                                          >  2:3344281/1‑149 (MQ=255)
                                                                                                        aCTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCg                                         <  2:848187/149‑1 (MQ=255)
                                                                                                                  agtGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACAt                               <  1:4483977/147‑1 (MQ=255)
                                                                                                                        aaaaGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATccgccg                         >  1:4008196/1‑149 (MQ=255)
                                                                                                                          aaGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGa                       >  1:349592/1‑149 (MQ=255)
                                                                                                                          aaGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGa                       >  2:3846148/1‑149 (MQ=255)
                                                                                                                             cTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGa                    <  2:3290668/149‑1 (MQ=255)
                                                                                                                             cTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGa                    <  1:2288847/149‑1 (MQ=255)
                                                                                                                              ttCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCTGACGAt                   <  2:4421386/149‑1 (MQ=255)
                                                                                                                              ttCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGAt                   <  1:2691285/149‑1 (MQ=255)
                                                                                                                               tCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATc                  <  1:343541/149‑1 (MQ=255)
                                                                                                                                cACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcg                 >  1:861182/1‑149 (MQ=255)
                                                                                                                                cACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcg                 >  1:1016954/1‑149 (MQ=255)
                                                                                                                                 aCGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgc                <  2:1192458/149‑1 (MQ=255)
                                                                                                                                  cGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgcc               <  1:1068206/149‑1 (MQ=255)
                                                                                                                                    aaTTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgccgc             >  1:1439531/1‑149 (MQ=255)
                                                                                                                                    aaTTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgccgc             >  2:398790/1‑149 (MQ=255)
                                                                                                                                     aTTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCgccgcc            <  2:1248016/149‑1 (MQ=255)
                                                                                                                                        tgtgCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGGTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGccaac         <  2:1461749/149‑1 (MQ=255)
                                                                                                                                          tgCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGccaaccc       <  1:3921412/149‑1 (MQ=255)
                                                                                                                                            ccACGGTGTTTACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGccaacccaa     >  2:2889821/1‑149 (MQ=255)
                                                                                                                                            ccACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGccaacccaa     <  2:2245045/149‑1 (MQ=255)
                                                                                                                                             cACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGccaacccaac    <  2:1608407/149‑1 (MQ=255)
                                                                                                                                              aCGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGCCAACCCAACt   <  2:1283035/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGCCAACCCAACTG  >  CP009273/3865717‑3866008

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: