Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,936,360 T→A 7.3% V259V (GTT→GTA) ‡ fucI → L‑fucose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,936,3600TA 7.3% 281.6 / 6.5 93V259V (GTT→GTA) ‡fucIL‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base T (43/43);  new base A (4/3);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

GACCGAACTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTGGACCGA  >  NC_000913/2936213‑2936503
                                                                                                                                                   |                                                                                                                                               
gACCGAACTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCt                                                                                                                                               >  2:1386082/1‑150 (MQ=255)
   cGAACTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTgcg                                                                                                                                            <  2:628431/150‑1 (MQ=255)
     aaCTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGACACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGcgcg                                                                                                                                          >  2:979125/1‑150 (MQ=255)
       cTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGAACTGcgcggc                                                                                                                                        <  1:987429/150‑3 (MQ=255)
        tGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGAACTGcgcggct                                                                                                                                       <  2:1385278/150‑4 (MQ=255)
         gCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAg                                                                                                                                      <  1:1040259/150‑1 (MQ=255)
         gCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAg                                                                                                                                      <  1:459212/150‑1 (MQ=255)
            gcGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGttt                                                                                                                                   <  1:613859/149‑1 (MQ=255)
                cGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTg                                                                                                                               >  1:563502/1‑150 (MQ=255)
                   aTCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATg                                                                                                                            <  2:1431042/150‑1 (MQ=255)
                     cGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGc                                                                                                                          <  1:376953/150‑1 (MQ=255)
                     cGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGc                                                                                                                          <  1:1289625/150‑1 (MQ=255)
                         cAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAg                                                                                                                                      <  1:1111926/134‑1 (MQ=255)
                         cAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAg                                                                                                                                      >  2:1111926/1‑134 (MQ=255)
                         cAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGAtg                                                                                                                      <  1:406671/150‑1 (MQ=255)
                          agaagaTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGAtgt                                                                                                                     <  1:1299775/150‑1 (MQ=255)
                                  tACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAACATTACTGATGGCGATGTCTATCcgc                                                                                                             >  1:324833/1‑150 (MQ=255)
                                   acgacgAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAATATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCcgcg                                                                                                            <  1:816096/150‑1 (MQ=255)
                                   acgacgAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCcgcg                                                                                                            <  1:1120981/150‑1 (MQ=255)
                                    cgccgaAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGCTGGCGATGTGTATCCGCGa                                                                                                           >  2:1257719/4‑150 (MQ=255)
                                     gacgaAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGAc                                                                                                          <  1:250210/150‑1 (MQ=255)
                                     gacgaAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGAc                                                                                                          >  2:763416/1‑150 (MQ=255)
                                       cgaAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACat                                                                                                        <  1:1514847/150‑1 (MQ=255)
                                        gaAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACatg                                                                                                       <  1:581776/150‑1 (MQ=255)
                                        gaAGCCGAATTGGAAATGACACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACatg                                                                                                       <  1:979125/150‑1 (MQ=255)
                                         aaGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACatga                                                                                                      <  1:844547/150‑1 (MQ=255)
                                           gCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACatgatg                                                                                                    >  2:464611/1‑150 (MQ=255)
                                            ccGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGc                                                                                                   <  1:464611/150‑1 (MQ=255)
                                                   ggAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCaa                                                                                            >  1:1533788/1‑150 (MQ=255)
                                                   ggAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAACCCGCGCAGTTCTGCGCCAAAGTTTATTGATGGCGATGTGTATCCGCGACATGATGCAAGGCaa                                                                                            >  2:1524719/1‑150 (MQ=255)
                                                    gAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAAc                                                                                           <  2:563502/150‑1 (MQ=255)
                                                    gAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAAc                                                                                           <  2:174158/150‑1 (MQ=255)
                                                     aaaTGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACa                                                                                          <  2:67920/150‑1 (MQ=255)
                                                       aTGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGc                                                                                        >  1:775248/1‑150 (MQ=255)
                                                        tGGCACTGGCCTGGGCTGTTAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCa                                                                                       >  1:1089146/1‑150 (MQ=255)
                                                           cACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAAc                                                                                    <  2:453469/150‑1 (MQ=255)
                                                           cACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAAc                                                                                    <  1:1367698/150‑1 (MQ=255)
                                                             cTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTg                                                                                  >  1:195422/1‑150 (MQ=255)
                                                               ggCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGc                                                                                >  1:1187663/1‑150 (MQ=255)
                                                                gCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGcc                                                                               <  1:1042687/150‑1 (MQ=255)
                                                                gCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGcc                                                                               <  2:533152/150‑1 (MQ=255)
                                                                   tGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGat                                                                            >  1:1428094/1‑150 (MQ=255)
                                                                    gggCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGata                                                                           >  2:334046/1‑150 (MQ=255)
                                                                    gggCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGata                                                                           >  1:1038487/1‑150 (MQ=255)
                                                                      gCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATAtt                                                                         >  2:678752/1‑150 (MQ=255)
                                                                      gCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATAtt                                                                         >  1:1484557/1‑150 (MQ=255)
                                                                          aTAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTc                                                                     <  2:775248/150‑1 (MQ=255)
                                                                            aaaaaCTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTcgc                                                                   >  1:1127963/1‑150 (MQ=255)
                                                                             aaaaCTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGATCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGa                                                                             <  2:1533932/139‑1 (MQ=255)
                                                                             aaaaCTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGa                                                                             >  1:1533932/1‑139 (MQ=255)
                                                                              aaaCTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGt                                                                 >  2:1262805/1‑150 (MQ=255)
                                                                               aaCTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTg                                                                >  2:1431989/1‑150 (MQ=255)
                                                                                 cTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGTGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGga                                                              >  2:154018/1‑150 (MQ=255)
                                                                                  ttCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaa                                                             <  1:373156/150‑1 (MQ=255)
                                                                                   tCCTCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCATAGCCGCGCAGTTCTGCGCGAAAGTTTACTGTTGGCGATGTGTATCCGAGACATGATGCAATGCAACAGCAAACTGGCCGATTTTGGTCGCGTGgaag                                                            <  2:1340899/150‑1 (MQ=255)
                                                                                   tCCGCTATGGCGAAGATGAAAATAACAAAAAGTATAAACGTAATGACGAGAAAAGCCGCGCAGTTATGCGCGAAAGTTTACTGATGGCGATGTGTATACGAGACATGATGCACGGCCACAGCAAACTGGCCGATATGGGTCGCGTGgaag                                                            >  2:730505/1‑150 (MQ=255)
                                                                                            gCGAAGATGAAAATAACAACCAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTg                                                   >  1:460849/1‑150 (MQ=255)
                                                                                             cGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTgg                                                  <  1:1262805/150‑1 (MQ=255)
                                                                                              gAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGc                                                 <  2:1187663/150‑1 (MQ=255)
                                                                                                 gATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTAc                                              <  1:1386082/150‑1 (MQ=255)
                                                                                                     aaaaTAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGTTATTGGTCGCGTGGAAGAATCACTTGGCTACAACg                                          >  2:1398084/1‑150 (MQ=255)
                                                                                                         tAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGACCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCAt                                      <  1:964925/150‑1 (MQ=255)
                                                                                                         tAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGAGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGATACAAAGCCAt                                      >  1:271868/1‑150 (MQ=255)
                                                                                                           acaaacaGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGc                                                 <  2:1326625/137‑1 (MQ=255)
                                                                                                           acaaacaGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGc                                                 >  1:1326625/1‑137 (MQ=255)
                                                                                                               acaGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGc                                >  2:1200582/1‑150 (MQ=255)
                                                                                                               acaGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGc                                >  2:603061/1‑150 (MQ=255)
                                                                                                                caGTATCAACGTAATGCCGAGCAACGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGAGATTGGGCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCg                               <  1:730505/150‑1 (MQ=255)
                                                                                                                     tCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGACGAATCAGTTGGCTACAACGCCATCGCTGCGGGCtt                          >  2:1298675/1‑150 (MQ=255)
                                                                                                                       aaCGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTcc                        >  2:620141/1‑150 (MQ=255)
                                                                                                                       aaCGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTcc                        >  1:1042658/1‑150 (MQ=255)
                                                                                                                       aaCGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTcc                        <  1:763416/150‑1 (MQ=255)
                                                                                                                           tAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAggg                    <  1:1053920/150‑1 (MQ=255)
                                                                                                                           tAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAggg                    <  1:1431989/150‑1 (MQ=255)
                                                                                                                           tAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAggg                    <  1:620141/150‑1 (MQ=255)
                                                                                                                              tGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCa                 >  2:496290/1‑150 (MQ=255)
                                                                                                                                ccGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAAc               >  1:1266933/1‑150 (MQ=255)
                                                                                                                                 cGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGAAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGATTCCAGGGGCAACg              >  2:571972/1‑150 (MQ=255)
                                                                                                                                   aGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGTGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTc            <  1:154018/150‑1 (MQ=255)
                                                                                                                                   aGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTc            <  2:1484557/150‑1 (MQ=255)
                                                                                                                                    gCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTcc           >  2:755062/1‑149 (MQ=255)
                                                                                                                                    gCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCa           >  2:380237/1‑150 (MQ=255)
                                                                                                                                    gCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCa           >  2:736101/1‑150 (MQ=255)
                                                                                                                                      actttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACt         >  2:218192/6‑150 (MQ=255)
                                                                                                                                      aaaGCCGCGCAGTTCTGCGCGAAAGTTTTCTGATGGCGTTGTGTATCCGCGACAAGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACACCTTCGCTGCGGGCTTCCAGGGGCAACGTCACt         <  2:1449271/150‑1 (MQ=255)
                                                                                                                                       ctttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTg        >  1:1481115/5‑150 (MQ=255)
                                                                                                                                       aaGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTg        >  1:677184/1‑150 (MQ=255)
                                                                                                                                        tttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTgg       >  2:911312/4‑150 (MQ=255)
                                                                                                                                        aGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTgg       <  2:94969/150‑1 (MQ=255)
                                                                                                                                         ttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGACGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTGGa      >  1:649282/3‑150 (MQ=255)
                                                                                                                                           cgcgCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTGGAcc    >  2:211344/1‑150 (MQ=255)
                                                                                                                                            gcgcAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGACGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTGGACCg   <  2:649282/150‑1 (MQ=255)
                                                                                                                                             cgcAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTGGACCGa  >  1:353294/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GACCGAACTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTGGACCGA  >  NC_000913/2936213‑2936503

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑