Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,724 G→T 14.8% intergenic (‑41/+267) gatY ← / ← fbaB D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,170,7240GT14.8% 432.0 / 5.1 165intergenic (‑41/+267)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base G (51/80);  new base T (21/9);  total (73/92)
Fisher's exact test for biased strand distribution p-value = 2.34e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

TCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACA  >  CP009273/2170577‑2170870
                                                                                                                                                   |                                                                                                                                                  
  tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaa                                                                                                                                                  <  2:1627919/148‑1 (MQ=255)
  tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaa                                                                                                                                                  <  2:551868/148‑1 (MQ=255)
   accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                                 <  2:2680173/148‑1 (MQ=255)
   accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTCAAAATATAATGaaa                                                                                                                                                 >  2:774986/1‑148 (MQ=255)
    ccaccaCTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                                <  1:3509263/148‑1 (MQ=255)
    ccaccaCTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                                <  2:14376/148‑1 (MQ=255)
    ccaccaCTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                                <  2:143401/148‑1 (MQ=255)
    ccaccaCTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                                <  1:2659024/148‑1 (MQ=255)
     caccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAatt                                                                                                                                               <  2:2730147/148‑1 (MQ=255)
     caccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAatt                                                                                                                                               <  1:927674/148‑1 (MQ=255)
     caccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAatt                                                                                                                                               <  2:2609852/148‑1 (MQ=255)
     caccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAatt                                                                                                                                               <  2:637766/148‑1 (MQ=255)
      accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAatta                                                                                                                                              <  1:670425/148‑1 (MQ=255)
       ccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattat                                                                                                                                             <  2:2231789/148‑1 (MQ=255)
       ccacGTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCGGTCGTTTGTTTTCGATTTCAAAATATAATGAAattat                                                                                                                                             <  2:932476/148‑1 (MQ=255)
         acTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTAttt                                                                                                                                           <  1:1647930/148‑1 (MQ=255)
          cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTg                                                                                                                                          <  1:3283377/148‑1 (MQ=255)
          cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTg                                                                                                                                          <  1:3282297/148‑1 (MQ=255)
          cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATTAAATTAATTg                                                                                                                                          >  2:3499553/1‑148 (MQ=255)
          cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTCAAAATAAAATTAAAATAtttt                                                                                                                                          >  2:659600/1‑147 (MQ=255)
           ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGt                                                                                                                                         <  2:353091/148‑1 (MQ=255)
            tGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAAATAAAATAtttttt                                                                                                                                        >  2:3006374/1‑148 (MQ=255)
             gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttt                                                                                                                                       <  2:166265/148‑1 (MQ=255)
              cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGtttt                                                                                                                                      <  2:2005874/148‑1 (MQ=255)
              cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGtttt                                                                                                                                      <  2:3372983/148‑1 (MQ=255)
                tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTATATTTCAAAAAAAAATTAAAATATTTGtttttt                                                                                                                                    >  1:577156/1‑147 (MQ=255)
                 cttttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaa                                                                                                                                   <  1:3401535/146‑1 (MQ=255)
                 cGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaa                                                                                                                                   <  2:2287729/148‑1 (MQ=255)
                 cGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaa                                                                                                                                   <  1:1616954/148‑1 (MQ=255)
                   tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAat                                                                                                                                 <  1:1957838/148‑1 (MQ=255)
                   tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAat                                                                                                                                 <  1:419572/148‑1 (MQ=255)
                   tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAat                                                                                                                                 <  1:1957827/148‑1 (MQ=255)
                    ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAata                                                                                                                                <  1:1112263/148‑1 (MQ=255)
                     tCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                               <  2:1900231/148‑1 (MQ=255)
                       gagaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCg                                                                                                                             <  2:1349087/148‑1 (MQ=255)
                          aTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCAAGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgaga                                                                                                                          <  1:1019021/148‑1 (MQ=255)
                           ttGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGAt                                                                                                                         <  2:4554/148‑1 (MQ=255)
                            tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATa                                                                                                                        <  2:2881693/148‑1 (MQ=255)
                            tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTCGGGTTCAAAATATAATGAAAATATATTTTTTTTAatatatatata                                                                                                                        >  1:1391982/1‑141 (MQ=255)
                                aaTATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGa                                                                                                                    <  2:946287/148‑1 (MQ=255)
                                 atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGAt                                                                                                                   <  2:2463424/148‑1 (MQ=255)
                                 atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGAt                                                                                                                   <  1:2510218/148‑1 (MQ=255)
                                 atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGGTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGAt                                                                                                                   <  2:2014423/148‑1 (MQ=255)
                                  tatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATc                                                                                                                  <  1:1331167/148‑1 (MQ=255)
                                     tGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGTTTTAAAAATATAAATAAATTATTTTTTTTTTTATATAGTGATTACGATaaca                                                                                                               >  2:3248066/1‑148 (MQ=255)
                                         tGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                                           <  2:906695/148‑1 (MQ=255)
                                          gCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGTAATTATTTGTTTTTTAATATAGAGATAACGATCACAAAAAc                                                                                                          >  2:516810/1‑148 (MQ=255)
                                             ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCTATACCACTTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGAc                                                                                                       <  2:1261847/148‑1 (MQ=255)
                                             ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGAc                                                                                                       <  2:3244856/148‑1 (MQ=255)
                                                  cGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGTTTTGAAAATATAAATTAAATATTTTTTTTTTTTTATATGGATAACGATCACAAAAAACACaaaaa                                                                                                  >  1:1852310/1‑145 (MQ=255)
                                                    cATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTTAAAATATAATTAAAATATTTTTTTTTAAAAATAGAGATTAAGATCAAAAAAAACACaaaaatt                                                                                                >  2:1789857/1‑143 (MQ=255)
                                                    cATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTGAAAATATAAATAAAATAATTTTTTTTTAATATAGTGATAACGATAACAAAAAAGAAAATATGa                                                                                                >  1:1107647/1‑148 (MQ=255)
                                                      tAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTGAAAATATAAATAAATTATTTTTTTTTTAATATAGAGATAACGATAACAAAAAAAACaaaaataaa                                                                                              >  1:257170/1‑141 (MQ=255)
                                                         ccgccgGGGTGTGCGTTTTTCAGCATCTGATTTGACGATACCACGTACATTTTCATATCCTGTCTTTTGTTTTCTTTTTAAAAATATAATTAATTTTTTTGTTTTTAAATATCGATATAACGAACAAAAAAACGACAATATGAAAatt                                                                                           <  1:564228/148‑1 (MQ=255)
                                                         ccgccgCGCTGTGCGTTGTTCAGCATGTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATACTTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAatt                                                                                           <  2:1111849/148‑1 (MQ=255)
                                                         ccgccgCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTTTTTGTTTTTAAATTTCGAGATACCGATCACAAAAACGACAATATGAAAatt                                                                                           <  2:1019248/148‑1 (MQ=255)
                                                         ccgccgCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGAACATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAACACGACAATATGAAAatt                                                                                           <  2:213589/148‑1 (MQ=255)
                                                           gccgcGCTGTGCGTTTTTCAGCATCTGCTTTGTCTATACCACGTAAATTTTCATATCCTGTCTTTTGTTTTCTATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAattat                                                                                         <  2:800445/148‑1 (MQ=255)
                                                           gccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCTATACCACGTACATTTTCATATCCTGTCTTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAattat                                                                                         <  2:1263165/148‑1 (MQ=255)
                                                           gccgcGCTGTGCGTTGTTCAGCATCTGATTTTTCTATACCACGTACATTTTAATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAattat                                                                                         <  2:188698/148‑1 (MQ=255)
                                                            ccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATTAAATTATTTTTTTTTAAAATTAGAGATAAAGATAACAAAAACGACAAAATTAAAattatt                                                                                        >  2:1100563/1‑148 (MQ=255)
                                                                  tgtgCGTTGACCCCCATCTGCTTTGCCGAGCCCACGTCCATTTCCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATCAAATTATTTGTTTTTAAATAGCGAGAAAACGATCACAAAAACGACAATATGAAAATTATTCgagga                                                                                  <  1:1114912/148‑1 (MQ=255)
                                                                    tgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTGAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAAGATCACAAAAACCACAATATGAAAAATATTCGAGGAGt                                                                                >  2:1712865/1‑148 (MQ=255)
                                                                           ttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGATTTGAAAAAATAATGAAAATATTTTTTTTTAAATATAGAGATAAAGATAAAAAAAAACACAAAATTAAAAATATTATAGGAGTGAAAggg                                                                         >  2:3366304/1‑147 (MQ=255)
                                                                           ttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTGTTGGTTTTAAAAAAAAAAATAAATTAATTTTTTTTAAAAATAGAGATAAACAACACAAAAAACAAAAAATGAAAAATATTAGAGGGGGGAAAGGc                                                                         >  1:3125565/1‑148 (MQ=255)
                                                                               gCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATTAAATTATTTGTTTTTTAATATAGAGATAAAGATCAAAAAAAAGACAAAATGTAAATTATTAGAGGAGTGAAAGGCcaaa                                                                     >  2:2857653/1‑148 (MQ=255)
                                                                               gCATCTGCTTTGACGATACCACGTAATTTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATCAAATTATTTGTTTTTAAATATCGATATAACGATAACAAACACGACAATACGAAAATTATTCGAGGAGTGAAAGGCaaaa                                                                     <  1:1404837/148‑1 (MQ=255)
                                                                                              aTACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTGAAAATATAATGAAAATATTTGTTTTTTAATATCGAGATAAAGATCACCAAAACGACCAAATGAAAATTATTGGCGGAGGGGAAGGCCAAAAAAAGGCCTCCCgcg                                                      >  2:3260097/1‑146 (MQ=255)
                                                                                                cccccGTAAATTTTATTATCCTGTGGTTTTTTTTCTTTTTAAAAATATATTAAATTTTTTTTTTTTTAAATATCTATATAACTATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAg                                                    <  1:2513198/144‑1 (MQ=255)
                                                                                                 ccACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTCAAAATATAATGAAATTTTTTGTTTTTTAATTTGGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGTAAGGCAAAAAAACGGCCTCCCGCTAgg                                                   >  2:281113/1‑148 (MQ=255)
                                                                                                      tACATTTTCATATCCTGTCGTGTGTTTTGGGTTTCAAAATATAAATTAAATATTTTTTTTTTAATATCGTGATAAAGATCACAAAAAACACAATATGAAAATTTTTGGCGGGGTGGAAGGCAAAAAAACGGCCGCCCGCTAGGGCAGc                                              >  2:2173786/1‑148 (MQ=255)
                                                                                                           tttCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                         <  2:3143400/148‑1 (MQ=255)
                                                                                                            ttCATATCCTGTCGTTTGTTTTGGATTTAAAAATAAAATTAAATTATTTTTTTTTTAATATCGAGATAACGATAACAAAAAAGAACAAATTAAAATTATTAGAGGAGTGTAATGGAAAAAAAAGGCCTCCCGCTAGGGAAGCCGGAGc                                        >  2:1677030/1‑148 (MQ=255)
                                                                                                                      gTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATTAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                              >  2:1366138/1‑148 (MQ=255)
                                                                                                                          tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                          >  1:3117812/1‑148 (MQ=255)
                                                                                                                           ttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAACATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt                         >  2:2460549/1‑148 (MQ=255)
                                                                                                                           ttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt                         >  2:1784209/1‑148 (MQ=255)
                                                                                                                           ttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGgtt                         >  2:2460559/1‑148 (MQ=255)
                                                                                                                            tgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                        >  1:2837292/1‑148 (MQ=255)
                                                                                                                             gtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt                       >  2:362534/1‑148 (MQ=255)
                                                                                                                               tttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTGAt                     >  2:2248385/1‑148 (MQ=255)
                                                                                                                                 tCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                   >  1:3019851/1‑148 (MQ=255)
                                                                                                                                     tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa               >  2:745648/1‑148 (MQ=255)
                                                                                                                                     tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa               >  2:2489869/1‑148 (MQ=255)
                                                                                                                                     tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa               >  2:2488992/1‑148 (MQ=255)
                                                                                                                                     tttAAAAATATATTTAAATTATTTTTTTTTAAATATCGAGATAACGAAAACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa               <  2:964210/148‑1 (MQ=255)
                                                                                                                                      ttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt              >  1:2910852/1‑148 (MQ=255)
                                                                                                                                      ttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt              >  1:1476193/1‑148 (MQ=255)
                                                                                                                                        cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGc            >  1:91242/1‑148 (MQ=255)
                                                                                                                                        aaaaataaaatGAATTTATTTTTTTTTAAATATCGAGATAACGAACACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGc            <  1:226418/147‑1 (MQ=255)
                                                                                                                                          aaaTATAATGAAATTATTTTTTTTTAAATATCGAGATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg          <  1:3184476/148‑1 (MQ=255)
                                                                                                                                          aaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGTAAGGCAAAAAAACGGCCTCCCGAGAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg          >  1:3015162/1‑148 (MQ=255)
                                                                                                                                          aaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg          >  2:1324150/1‑148 (MQ=255)
                                                                                                                                            atatAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtg        >  2:506905/1‑148 (MQ=255)
                                                                                                                                             tataATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtgt       >  1:3508646/1‑148 (MQ=255)
                                                                                                                                               taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaa     >  2:2584508/1‑148 (MQ=255)
                                                                                                                                                aaTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa    >  2:1094927/1‑148 (MQ=255)
                                                                                                                                                aaTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa    >  2:13552/1‑148 (MQ=255)
                                                                                                                                                  tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACa  >  1:3215625/1‑148 (MQ=255)
                                                                                                                                                  tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACa  >  2:3068496/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACA  >  CP009273/2170577‑2170870

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