Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 42.0% 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,9850TC42.0% 287.0 / inf 317intergenic (‑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 (105/79);  new base C (126/7);  total (231/86)
Fisher's exact test for biased strand distribution p-value = 3.95e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-04

Only 100 of 323 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084
                                                                                                   |                                                                                                       
tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                         >  2:75385/1‑100 (MQ=255)
tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                         <  1:2929442/100‑1 (MQ=255)
  ctgcCGGATGGCGCTGCACGCGGAGGCACAGCGCCACGGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                       >  1:2363313/1‑100 (MQ=255)
  ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                       >  1:3288729/1‑100 (MQ=255)
     ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                    >  2:1244450/1‑100 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                 <  1:1244446/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                 <  1:2237018/100‑2 (MQ=37)
         aTGGCGCTGCACGCAGTGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                <  2:2171244/100‑3 (MQ=25)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                <  1:2585685/100‑1 (MQ=255)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                <  2:2704199/100‑1 (MQ=255)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGTATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                            >  1:3638596/1‑100 (MQ=255)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                            >  1:3143893/1‑100 (MQ=255)
             cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                           <  2:846729/100‑1 (MQ=255)
               atGCACGCAGAGGCACAGCGCCCCCGGGTATAGAGGTTGTTGACAAAGGGCGCTTTGTTTATGACGGATTCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                                         >  1:3377083/2‑100 (MQ=21)
                 gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                       >  1:1343124/1‑100 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                      >  1:152041/1‑100 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                      >  2:3239664/1‑100 (MQ=255)
                    cgatctggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                    <  1:556371/94‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                    >  2:3185546/1‑100 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCATg                                                                                    >  1:3846442/1‑100 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                   <  2:1161768/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                   <  2:1562290/100‑1 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                  <  2:1373877/100‑1 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCCCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                  <  1:3082814/100‑1 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                            >  1:2267071/1‑100 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                          <  2:1412288/100‑1 (MQ=255)
                                  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                      >  2:2606558/1‑100 (MQ=255)
                                             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAt‑ttt                                                        >  1:1074708/1‑100 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCg                                                      <  1:3501177/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCg                                                      <  2:2995149/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAggac                                                >  1:2389369/1‑98 (MQ=255)
                                                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGa                                             >  1:3744466/1‑100 (MQ=255)
                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGc                                        <  1:1976622/100‑1 (MQ=255)
                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTA‑CA‑GAAT‑TTTCGTAGGCCTGATAAGc                                        <  1:2470434/100‑1 (MQ=255)
                                                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCgg                                      <  1:2205020/100‑1 (MQ=255)
                                                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCa                                >  1:26943/1‑100 (MQ=255)
                                                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCAt                               >  2:3772392/1‑100 (MQ=255)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                              >  1:141880/1‑100 (MQ=255)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                              <  2:1791456/100‑1 (MQ=255)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                              <  1:1690462/100‑1 (MQ=255)
                                                                        tgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                              >  2:2264423/6‑99 (MQ=17)
                                                                         gtaggCCGGATAAGGCGTTCACGCTGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CCAT‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAg                            >  2:645012/5‑100 (MQ=18)
                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAg                            >  2:1889616/5‑100 (MQ=18)
                                                                         gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATGCAATACATTG‑CA‑GAAG‑TTTCGTAGACCTGATAAGCGGAGCGCATCAg                            >  1:2980673/5‑100 (MQ=18)
                                                                         gtaggCCGGATAAGGCGCTTACGCCGCATCCGGCCTAC‑AAGAACGTGCAAATTCAACAAATTGTGATGCAT‑CCT‑GTAGGCCTGATAAGCGTAGCGCATCa                             >  2:515598/5‑100 (MQ=18)
                                                                          taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAgg                           >  1:1483394/4‑100 (MQ=17)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTG‑CA‑GAGA‑TTATGTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:2553222/4‑100 (MQ=17)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAgg                           >  1:3266774/1‑100 (MQ=255)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAgg                           <  1:3801484/100‑1 (MQ=255)
                                                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAAGCGTGCAAATTCAATACATTG‑CA‑TGGG‑CCATGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:4033404/3‑100 (MQ=17)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTG‑CA‑GGAG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:1713708/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTG‑CAGGGAC‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAGGc                          >  2:3792127/3‑100 (MQ=18)
                                                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                          <  2:1966884/100‑1 (MQ=255)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                         >  2:2898305/1‑100 (MQ=255)
                                                                            ggctggataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTG‑CA‑CAGC‑TAGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2329378/11‑100 (MQ=25)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                        >  2:4283658/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                        >  2:693336/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCACTATATTG‑CA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:732366/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:966751/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTG‑CA‑ATCT‑TCTGGTAGGCCTGATAAGCATCGCGCATCAGGCa                         >  1:1574412/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCAATATATTG‑CA‑GAGA‑TTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:4098772/2‑100 (MQ=11)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATGCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTG‑AT‑CGAT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2679109/2‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑TG‑GGCT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2075891/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTG‑CAGGGAC‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                         >  2:2470596/2‑100 (MQ=18)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        <  2:3083076/100‑1 (MQ=255)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  2:2747060/1‑100 (MQ=18)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTG‑CA‑TGAG‑CCACGTAGTCCTGATAAGCGTAGCGCATCAGGCag                        >  1:1008699/1‑99 (MQ=17)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTG‑CA‑TGAG‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  1:3151441/1‑99 (MQ=17)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CCAT‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:2803502/1‑99 (MQ=21)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTG‑CA‑GGAG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:2227271/1‑99 (MQ=18)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                       >  1:1177795/9‑98 (MQ=18)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                       >  1:1376719/9‑98 (MQ=18)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GAGT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                       >  2:457522/9‑98 (MQ=17)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                       >  2:4384513/9‑98 (MQ=18)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTG‑CA‑GTTT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaat                       >  2:3873554/9‑98 (MQ=25)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTG‑CA‑AGAT‑TTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacg                       >  2:4129596/9‑98 (MQ=17)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTG‑CA‑GAGA‑TCATGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  2:4052065/9‑98 (MQ=14)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTG‑CA‑GGAG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                      >  2:1286656/9‑100 (MQ=18)
                                                                              ccGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTG‑CA‑GGAG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                      >  1:1279166/1‑100 (MQ=21)
                                                                              ccGGATAAGGCGGTCACTCCGCATCCGACCTAC‑‑AAATCATGCAAATCCAAAATATTG‑CA‑GGAG‑GTATGTATGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  1:4302543/1‑99 (MQ=18)
                                                                               cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTG‑CAGGGAC‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAGGCagct                      >  2:4234326/8‑97 (MQ=18)
                                                                               cggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                      >  2:2418941/8‑100 (MQ=18)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                    >  2:3501717/7‑98 (MQ=21)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GAGT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:3355097/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:2287956/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:3937798/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:1778834/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑TG‑GGCT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:3320088/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTG‑CAGGAAC‑GAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                     >  2:2121114/7‑96 (MQ=25)
                                                                                ggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                     >  2:1042106/7‑100 (MQ=18)
                                                                                ggataaggcgGTTACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTG‑CA‑ATCT‑CCGTGTAGGCATGATAAGCGTAGCGAATCAGGccat                       >  1:4190648/7‑95 (MQ=17)
                                                                                  aTGCGGCGTTAACGCCTTATCCTGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑TA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:802814/100‑1 (MQ=255)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:923938/100‑1 (MQ=255)
                                                                                   taaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                  >  1:3555881/4‑100 (MQ=18)
                                                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                  >  2:2022169/1‑100 (MQ=255)
                                                                                    gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                 <  1:1834323/100‑1 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  2:4160287/1‑100 (MQ=255)
                                                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               <  2:3556852/100‑1 (MQ=255)
                                                                                                ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  2:3395845/1‑100 (MQ=255)
                                                                                                  ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  1:209480/100‑1 (MQ=255)
                                                                                                   |                                                                                                       
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084

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