Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,936,359 T→A 9.9% V259D (GTT→GAT) ‡ fucI → L‑fucose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,936,3590TA 9.9% 299.2 / 20.4 91V259D (GTT→GAT) ‡fucIL‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base T (41/41);  new base A (5/4);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.45e-01

CCGAACTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTG  >  NC_000913/2936215‑2936497
                                                                                                                                                |                                                                                                                                          
ccGAACTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTg                                                                                                                                        >  4:703651/1‑149 (MQ=255)
             gcagttCGGTCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTAc                                                                                                                           <  3:689545/144‑1 (MQ=255)
                tATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGa                                                                                                                        >  4:695524/1‑149 (MQ=255)
                         gaataTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGtt                                                                                                                              <  1:245449/134‑1 (MQ=255)
                         gaagaTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATTAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGtt                                                                                                                              >  2:245449/1‑134 (MQ=255)
                         gaagaTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAGCAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGAtgt                                                                                                               <  3:1064594/149‑1 (MQ=255)
                         gaagaTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCt                                                                                                                                         >  3:1188874/1‑123 (MQ=255)
                         gaagaTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCt                                                                                                                                         <  4:1188874/123‑1 (MQ=255)
                         gaagaTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGAtgt                                                                                                               >  2:16917/1‑149 (MQ=255)
                         gaagaTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGAtgt                                                                                                               >  1:392480/1‑149 (MQ=255)
                            gaTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTa                                                                                                            <  2:392480/149‑1 (MQ=255)
                             aTTTACGACGAAGCCGAATTGGCAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTAt                                                                                                           <  1:1255073/149‑1 (MQ=255)
                               ttACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATcc                                                                                                         <  4:146041/149‑1 (MQ=255)
                                   gacgaAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGa                                                                                                     >  4:528196/1‑149 (MQ=255)
                                   gacgaAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGa                                                                                                     >  1:508160/1‑149 (MQ=255)
                                    acgaAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGAc                                                                                                    >  2:666884/1‑149 (MQ=255)
                                         gCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACatgat                                                                                               <  2:508160/149‑1 (MQ=255)
                                           cGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGc                                                                                             <  3:695524/149‑1 (MQ=255)
                                                  gAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTAc                                                                                                                           <  1:601098/112‑1 (MQ=255)
                                                  gAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTAc                                                                                                                           >  2:601098/1‑112 (MQ=255)
                                                  gAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCaa                                                                                      <  3:703651/149‑1 (MQ=255)
                                                    aaTGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACa                                                                                    <  1:666884/149‑1 (MQ=255)
                                                      tGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGtt                                                                                                                              >  3:1012975/1‑105 (MQ=255)
                                                      tGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGtt                                                                                                                              <  4:1012975/105‑1 (MQ=255)
                                                          aCTGGCCTGGGCTGATAAAAACTTCCGCTATGACGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGtt                                                                                                                                           <  2:1453285/88‑1 (MQ=255)
                                                          aCTGGCCTGGGCTGATAAAAACTTCCGCTATGACGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGtt                                                                                                                                           >  1:1453285/1‑88 (MQ=255)
                                                              gCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGa                                                                                                                  <  1:344734/109‑1 (MQ=255)
                                                              gCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGa                                                                                                                  >  2:344734/1‑109 (MQ=255)
                                                                 tGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGttt                                                                                                                             >  3:1087033/1‑95 (MQ=255)
                                                                 tGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGttt                                                                                                                             <  4:1087033/95‑1 (MQ=255)
                                                                      tGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGAACTGcgcgg                                                                                                                                   <  1:1312856/84‑2 (MQ=255)
                                                                      tGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGAACTGcgcgg                                                                                                                                   >  2:1312856/1‑83 (MQ=255)
                                                                         tAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCaa                                                                                           <  3:287801/121‑1 (MQ=255)
                                                                         tAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCaa                                                                                           >  4:287801/1‑121 (MQ=255)
                                                                         tAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTc                                                               >  4:508887/1‑149 (MQ=255)
                                                                            aaaCTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTcgcg                                                            >  4:1423725/1‑149 (MQ=255)
                                                                               cTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTc                                                               >  3:676556/1‑143 (MQ=255)
                                                                               cTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTc                                                               <  4:676556/143‑1 (MQ=255)
                                                                                ttCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGga                                                        >  4:165641/1‑149 (MQ=255)
                                                                                ttCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGga                                                        >  1:1388666/1‑149 (MQ=255)
                                                                                ttCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGga                                                        >  1:507996/1‑149 (MQ=255)
                                                                                ttCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGga                                                        >  3:359752/1‑149 (MQ=255)
                                                                                 tCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaa                                                       <  3:528196/149‑1 (MQ=255)
                                                                                        tGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCAc                                                <  3:508887/149‑1 (MQ=255)
                                                                                          gCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACtt                                              <  2:507996/149‑1 (MQ=255)
                                                                                          gCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACtt                                              <  4:359752/149‑1 (MQ=255)
                                                                                          gCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACtt                                              <  2:1388666/149‑1 (MQ=255)
                                                                                          gCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACtt                                              <  1:16917/149‑1 (MQ=255)
                                                                                            gAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTgg                                            >  1:727260/1‑149 (MQ=255)
                                                                                             aaGATGAAAATAACGAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACatga                                                                                                >  3:101828/1‑96 (MQ=255)
                                                                                             aaGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACatga                                                                                                <  4:101828/96‑1 (MQ=255)
                                                                                                aTGAAAATAACAAACAGTATCAACGTAATGCCGAGGAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGCGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaa                                                       <  4:1004883/134‑1 (MQ=255)
                                                                                                aTGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATCGGTCGCGTGgaa                                                       >  3:1004883/1‑134 (MQ=255)
                                                                                                 tGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAg                                                                                                                                <  1:326750/60‑1 (MQ=255)
                                                                                                 tGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAg                                                                                                                                >  2:326750/1‑60 (MQ=255)
                                                                                                 tGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGcc                                                                         <  4:1465397/115‑1 (MQ=255)
                                                                                                 tGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGcc                                                                         >  3:1465397/1‑115 (MQ=255)
                                                                                                  gAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACaa                                      >  4:469586/1‑149 (MQ=255)
                                                                                                    aaaTAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACg                                    <  3:469586/149‑1 (MQ=255)
                                                                                                     aaTAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGc                                   >  3:1406406/1‑149 (MQ=255)
                                                                                                      aTAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGcc                                  <  3:1423725/149‑1 (MQ=255)
                                                                                                        aacaaacaGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCAt                                <  3:165641/149‑1 (MQ=255)
                                                                                                                 tATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTg                                                                            <  3:1046574/96‑1 (MQ=255)
                                                                                                                 tATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTg                                                                            >  4:1046574/1‑96 (MQ=255)
                                                                                                                   tCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACa                                                                                    >  3:582450/1‑86 (MQ=255)
                                                                                                                   tCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACa                                                                                    <  4:582450/86‑1 (MQ=255)
                                                                                                                   tCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCt                     >  2:508892/1‑149 (MQ=255)
                                                                                                                     aaCGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACtt                                              >  3:253473/1‑122 (MQ=255)
                                                                                                                     aaCGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACtt                                              <  4:253473/122‑1 (MQ=255)
                                                                                                                     aaCGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATAGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTc                   >  3:1210799/1‑149 (MQ=255)
                                                                                                                      aCGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGGAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTcc                  <  4:1210799/149‑1 (MQ=255)
                                                                                                                      aCGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATc                                                  <  1:165459/117‑1 (MQ=255)
                                                                                                                      aCGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATc                                                  >  2:165459/1‑117 (MQ=255)
                                                                                                                        gTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGc                                                                          >  4:744965/1‑91 (MQ=255)
                                                                                                                        gTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGc                                                                          <  3:744965/91‑1 (MQ=255)
                                                                                                                        gTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaa                                                       <  3:41313/110‑1 (MQ=255)
                                                                                                                        gTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaa                                                       >  4:41313/1‑110 (MQ=255)
                                                                                                                        gTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAg                >  3:1329392/1‑149 (MQ=255)
                                                                                                                              ccGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCaa          >  4:558398/1‑149 (MQ=255)
                                                                                                                                 aGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCgg                        >  2:484656/1‑132 (MQ=255)
                                                                                                                                 aGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCgg                        <  1:484656/132‑1 (MQ=255)
                                                                                                                                   cAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCa     <  4:1329392/149‑1 (MQ=255)
                                                                                                                                   aactttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCt                            >  1:340078/7‑126 (MQ=255)
                                                                                                                                   aactttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCt                            <  2:340078/120‑1 (MQ=255)
                                                                                                                                      tttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGc                                   >  3:580421/4‑116 (MQ=255)
                                                                                                                                      tttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGc                                   <  4:580421/113‑1 (MQ=255)
                                                                                                                                      tttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCt                            <  4:655004/120‑1 (MQ=255)
                                                                                                                                      tttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCt                            >  3:655004/4‑123 (MQ=255)
                                                                                                                                      tttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTg  >  1:225425/4‑149 (MQ=255)
                                                                                                                                          gcgcAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACt   >  4:312284/1‑144 (MQ=255)
                                                                                                                                          gcgcAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACt   <  3:312284/144‑1 (MQ=255)
                                                                                                                                           cgcAGTTCTGCGCGTAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGatat                                                                    <  4:1000358/78‑1 (MQ=255)
                                                                                                                                           cgcAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGatat                                                                    >  3:1000358/1‑78 (MQ=255)
                                                                                                                                                |                                                                                                                                          
CCGAACTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTG  >  NC_000913/2936215‑2936497

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: