Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 24.3% intergenic (‑195/+113) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC24.3% 517.8 / 168.6 259intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (88/108);  new base C (60/3);  total (148/111)
Fisher's exact test for biased strand distribution p-value = 3.94e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-02

Only 100 of 264 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084
                                                                                                  |                                                                                                              
 ctgagggATGGCGCTGCACGCAGAGGCAGCGCGTAAACGGGTGTACAGATTGTTGTCAAAGTGCGATTTGTGTATGGCGTAGCCCGCGTTAACACCttc                                                                                                               <  2:362567/94‑2 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                          <  2:2211313/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                        <  1:3776470/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCGCATCCGGc                                                                                                         <  2:1624355/99‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                        <  1:4138406/100‑2 (MQ=37)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGCGCGCTTCGTTTATGCCCGATGCGGCGCTAACGCCTTCCCCGGCCTa                                                                                                      >  2:1857371/1‑100 (MQ=25)
          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                                     >  2:348834/1‑100 (MQ=255)
            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGGTAACGCCTTATCCGGCCTAC‑aa                                                                                                  >  1:105899/1‑100 (MQ=255)
               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                               <  1:3178730/100‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                             <  1:3737482/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                           >  1:2927163/1‑100 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                          <  2:3907159/100‑1 (MQ=255)
                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                         <  1:3900806/100‑1 (MQ=255)
                         ggCCCAGCCCCACCGAGTATACATATCGTTGACAAAGTGCGCTTTGTTTGTGCCGGAGGCGGCGTTACCGACTTATCCGGCCTAC‑ACAACCGTGCa                                                                                         <  2:3158118/96‑1 (MQ=17)
                         agCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                                     >  2:1754591/2‑100 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                                    <  2:1387920/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                                 <  2:3087705/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                                 <  1:144595/100‑1 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                                                <  1:3314644/100‑1 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCAGCCTAC‑AAAATCGTGCAAATTCAATa                                                                                <  1:2713647/100‑1 (MQ=255)
                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                               >  1:2310467/1‑100 (MQ=255)
                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                               <  1:3291938/100‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCGAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑g                                                                      <  2:1675371/100‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGAATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                           <  2:2903065/97‑1 (MQ=255)
                                       gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑Ga                                                                     <  1:2768258/100‑1 (MQ=255)
                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTCC‑AAAATCGTGCAAATTCAATACATTGc                                                                          >  2:2491604/1‑96 (MQ=255)
                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑Gaa                                                                    >  1:919460/1‑100 (MQ=255)
                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAt‑‑t‑‑‑t                                                            <  1:2488960/100‑1 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGTATTCGGCGTTAACTCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAt‑‑t‑‑‑t‑t                                                          >  1:2268689/1‑100 (MQ=255)
                                              cgatctgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑g                                                      <  2:2288436/95‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTa                                                    >  1:2568812/1‑100 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTa                                                    <  2:4151717/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTa                                                    >  1:3397320/1‑100 (MQ=255)
                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGc                                                 >  2:2981095/1‑100 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGcc                                                <  1:2124849/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGcc                                                >  1:4052532/1‑100 (MQ=255)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTg                                              >  2:2378451/1‑100 (MQ=255)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTg                                              >  2:857927/1‑100 (MQ=255)
                                                      cgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTg                                              >  2:281776/2‑100 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATa                                           <  1:871671/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATa                                           <  1:4047399/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATa                                           <  2:2646751/100‑1 (MQ=255)
                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATaa                                          >  2:3743007/1‑100 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCg                                       >  2:2916819/1‑100 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCgg                                      >  2:75980/1‑100 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGa                                     >  1:233793/1‑99 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAg                                    <  2:3729640/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  1:281779/100‑1 (MQ=255)
                                                                   attttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCTAATTCAATATATTGC‑A‑GAAA‑‑C‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAgcgc                                 >  1:1406379/10‑100 (MQ=17)
                                                                    ttGTTTATGCCGGATGCGCCGTTAACGCCTTATCCGGCCTAC‑AAACTCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCa                                <  1:1157764/100‑1 (MQ=255)
                                                                    atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GAATC‑C‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCAt                               >  1:2671933/9‑100 (MQ=18)
                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCAt                               <  2:1924994/100‑1 (MQ=255)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑A‑TGAG‑‑C‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCa                                >  1:1294540/7‑98 (MQ=17)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                             <  2:2395782/100‑1 (MQ=255)
                                                                        gtaggCCGGATAAGGCGCTTACGCCGCATCCGGCCTAC‑AAGAACGTGCAAATTCAACAAATTGTGATGCAT‑‑C‑‑‑C‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCa                             >  1:1023285/5‑100 (MQ=18)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGC‑A‑‑ATA‑‑T‑‑‑G‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:3794943/4‑100 (MQ=18)
                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTCC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGACCTGATAAGAGGAGCGCATCAgg                           <  2:2567147/100‑1 (MQ=255)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑A‑GAAA‑‑C‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:3633099/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑‑AAT‑‑TCAAC‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                           >  1:3687641/3‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGC‑ATGGGC‑‑C‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:3332077/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGT‑G‑GGCT‑‑C‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:1516650/2‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑A‑GGAG‑‑C‑‑‑T‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:3830939/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑A‑GGAG‑‑C‑‑‑T‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2057983/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑A‑ATCT‑‑C‑‑‑C‑GT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:3660336/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                          >  1:1366761/2‑99 (MQ=18)
                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  2:933234/1‑100 (MQ=255)
                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  1:1754180/1‑100 (MQ=255)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGC‑A‑GAAA‑‑T‑‑‑C‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  1:1641078/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GAATC‑C‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  1:743393/1‑100 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑T‑CGAT‑‑C‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  1:3517283/1‑99 (MQ=21)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑AAT‑‑C‑‑‑T‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  1:3561992/9‑98 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑AGGGAC‑‑T‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagc                       >  1:2655868/9‑98 (MQ=18)
                                                                             ccGGATGCGGAGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGACCTGATAAGCTGAGCGCATAACGCAtt                       >  2:2547629/1‑100 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GAAT‑‑C‑‑‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                    >  1:4005093/7‑98 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑A‑‑ATTCAA‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                      >  1:2392307/7‑97 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑A‑AGAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                     >  1:2524187/7‑96 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑A‑GGAG‑‑C‑‑‑T‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:535835/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑AGGAAC‑‑G‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                     >  1:3390970/7‑96 (MQ=25)
                                                                               ggacaaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑A‑AGAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                     >  2:1197906/7‑96 (MQ=17)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑A‑GGAC‑‑C‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                    >  2:2828621/6‑100 (MQ=18)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:3438729/100‑1 (MQ=255)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:2947696/100‑1 (MQ=255)
                                                                                    aggcgTTAACGCCTTTTTTGGCCTAC‑TAAATCATGCAAATTCAATTCATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTAATAAGCGGAGCGCATCAGGCATTtttgc                  >  2:3554594/2‑98 (MQ=18)
                                                                                     ngcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATCtttgcctt               >  2:1041620/2‑100 (MQ=255)
                                                                                     ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               <  2:11452/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  2:320447/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  2:481015/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGACCAAC‑‑AAATTATGCAACTCCACAATATTGC‑A‑GAAA‑‑T‑‑‑T‑AA‑‑GTAGGACTTATAAGCATAACGCATAAAACAtttattcgttt              >  2:329578/1‑90 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TT‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  2:883385/100‑1 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  1:3215199/100‑1 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  2:1096380/100‑1 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt         >  1:3365096/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:1718252/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAACGCATCAGGCATTTTTGCCTTTGCCATCa       >  2:3967615/1‑100 (MQ=255)
                                                                                              gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGGCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      <  1:3960608/100‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCan      <  2:535133/100‑2 (MQ=255)
                                                                                               ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  1:3105986/1‑100 (MQ=255)
                                                                                               ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  1:3015510/1‑100 (MQ=255)
                                                                                                 ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  1:441288/100‑1 (MQ=255)
                                                                                                  tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  <  1:1890908/100‑1 (MQ=255)
                                                                                                  |                                                                                                              
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084

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