Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 32.4% 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,9850TC32.4% 370.8 / 225.7 241intergenic (‑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 (69/94);  new base C (73/5);  total (142/99)
Fisher's exact test for biased strand distribution p-value = 9.70e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.66e-01

Only 100 of 245 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705886‑3706083
                                                                                                   |                                                                                                             
tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                               >  1:1022670/1‑100 (MQ=255)
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                              >  2:23952/1‑100 (MQ=255)
 gctgcCGGATGACGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                              <  1:1093019/100‑1 (MQ=255)
  ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                             >  2:3712696/1‑100 (MQ=255)
       ggATGGCGCTACACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTCTACGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                        <  1:2461596/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                       <  1:727678/100‑2 (MQ=37)
         gtGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                      >  2:2910299/2‑100 (MQ=255)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                      <  2:382467/100‑1 (MQ=255)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                      <  1:3000303/100‑3 (MQ=25)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                      <  2:1302352/97‑3 (MQ=25)
                  cACGCAGAGGCACAGCGCCACCGGGTATACATATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCg                                                                                           <  2:505950/100‑1 (MQ=255)
                   aCGCAGCGCCACAGCGCCACCGGGTCTACAGATTGTTGACAAAGCGCGCATCCGTTTCACCCGATGCGACGTTAACGCCTTATCCCGCCT‑CC‑AACATc                                                                                            >  1:1735371/1‑98 (MQ=16)
                   aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGt                                                                                          <  2:1707329/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTg                                                                                         >  1:2574668/1‑100 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGc                                                                                        <  2:2958851/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGc                                                                                        <  1:2846191/100‑1 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGTCCT‑AC‑AAAATCGTGCa                                                                                       >  1:3459258/1‑100 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCa                                                                                       <  2:2011895/100‑1 (MQ=255)
                         agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAAt                                                                                    <  2:104045/100‑1 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAAt                                                                               <  2:2370798/100‑1 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAAt                                                                               <  2:401821/100‑1 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAAt                                                                               <  1:2910280/100‑1 (MQ=255)
                                  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACAt                                                                           >  1:1187763/1‑100 (MQ=255)
                                   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACAtt                                                                          >  1:112978/1‑100 (MQ=255)
                                    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTg                                                                         >  2:3558869/1‑100 (MQ=255)
                                       cGGGTATACAGATTGTTGACAAAGTGCTCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑g                                                                   <  1:998186/100‑1 (MQ=255)
                                       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑g                                                                   <  1:2811991/100‑1 (MQ=255)
                                           tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAt‑t                                                              >  2:3153657/1‑100 (MQ=255)
                                            ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGTGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAt‑t‑‑‑t                                                          <  2:138914/100‑1 (MQ=255)
                                            ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAt‑t‑‑‑t                                                          <  1:3570604/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑g                                                     <  2:646015/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑g                                                     <  2:1187769/100‑1 (MQ=255)
                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑Gt                                                    >  2:3621596/1‑100 (MQ=255)
                                                 gcttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTa                                                   >  2:2531750/3‑100 (MQ=255)
                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTa                                                   <  1:3621568/100‑1 (MQ=255)
                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTa                                                   <  1:1220134/100‑1 (MQ=255)
                                                  aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAg                                                  <  1:2373010/100‑1 (MQ=255)
                                                   gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAgg                                                 >  1:2204182/7‑100 (MQ=18)
                                                        gACAAAGTGCGCTTTGTTTATGCCAAATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTg                                             <  1:3108082/99‑1 (MQ=255)
                                                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGAt                                           <  1:272825/100‑1 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                          <  2:1156422/100‑1 (MQ=255)
                                                           aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATaa                                         <  1:848009/100‑1 (MQ=255)
                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGc                                       >  2:2984976/1‑100 (MQ=255)
                                                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                     >  1:3081429/1‑100 (MQ=255)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                   >  2:3616599/1‑100 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                  <  1:1773930/100‑1 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                 >  1:372536/1‑100 (MQ=255)
                                                                    tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                                >  2:1896797/1‑100 (MQ=255)
                                                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                               <  2:1435350/100‑1 (MQ=255)
                                                                     ttGTGTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAAGAGTGCTAATTAAATACATTGC‑‑T‑GTAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                     <  1:440247/94‑1 (MQ=17)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                             <  1:1275269/100‑1 (MQ=255)
                                                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            <  1:1201295/100‑1 (MQ=255)
                                                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            >  2:500551/1‑100 (MQ=255)
                                                                        tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:2562065/6‑100 (MQ=18)
                                                                        tgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            >  1:1597669/6‑100 (MQ=18)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           <  1:1171033/100‑1 (MQ=255)
                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GGAT‑A‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  2:1528971/5‑100 (MQ=17)
                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑AC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GGAT‑A‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:494504/5‑100 (MQ=18)
                                                                         gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           >  1:2095719/5‑100 (MQ=17)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GTTT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                          >  1:3339126/4‑100 (MQ=21)
                                                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2753248/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TGAG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:3307256/3‑100 (MQ=17)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑AAT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:3500583/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑‑AAACCGCGCAAATTCAATATATTGCGGA‑GAAA‑A‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                        >  1:1938436/3‑99 (MQ=21)
                                                                            ggtcggataaggcgTTCACGCCGCATCCGGCCC‑AC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAAC‑G‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:2762599/11‑100 (MQ=18)
                                                                            ggtcggataaggcgTTCACGCCGCATCCGACCC‑AC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAAC‑G‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:1978380/11‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GAGT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2468600/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑‑‑AAT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                         >  2:3188461/2‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCATGAAAATTCAATACATTGC‑‑A‑AGAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2696508/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑AC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑‑ATA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga                       >  2:1880054/2‑98 (MQ=21)
                                                                             gtcggataaggcgTTCACGCCGCATCCGACCC‑AC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAAC‑G‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  2:1288420/10‑99 (MQ=21)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  2:1675406/1‑99 (MQ=21)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAACCTTGCTAATTCAATATATTGC‑‑A‑GGGA‑C‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag                       >  2:2581391/1‑99 (MQ=18)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GAAA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:3334806/9‑98 (MQ=18)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GAGT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:746243/9‑98 (MQ=17)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GGAC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:2291139/9‑98 (MQ=18)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GTTT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                      >  1:2325546/9‑98 (MQ=25)
                                                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:1962587/1‑100 (MQ=255)
                                                                              ccGGATAAGGCGTTTACGCCGCATCTGGCCT‑AC‑AAATTCTTCCAAATTCAATATATTGC‑‑A‑GAAA‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  1:1059579/1‑98 (MQ=17)
                                                                              ccGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑AAT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  2:87263/1‑98 (MQ=18)
                                                                               cggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GAGT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                     >  2:2603850/8‑97 (MQ=17)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GAAA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1769186/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑AC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TGAG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtat                    >  1:2271363/7‑96 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCT‑AC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GGAG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:1724518/7‑100 (MQ=18)
                                                                                ggATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    <  1:1389415/100‑1 (MQ=255)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCT‑AC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TGAG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAAAGAATCAGGCAgtttt                   >  2:1362151/6‑100 (MQ=17)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:1847113/100‑1 (MQ=255)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  2:774892/100‑1 (MQ=255)
                                                                                   tGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 <  2:1774810/100‑1 (MQ=255)
                                                                                    gcggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                <  1:2951913/100‑1 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  2:3567587/1‑100 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:2109391/100‑1 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:2026018/100‑1 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:2218269/100‑1 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc          <  2:3461796/100‑1 (MQ=255)
                                                                                            aaCGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATAGGGCATTTTTGCCTTTGCCAt        >  2:1604981/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:2246440/100‑1 (MQ=255)
                                                                                              cGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa      <  1:1441139/100‑1 (MQ=255)
                                                                                               gCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  2:3389178/100‑1 (MQ=255)
                                                                                                   tATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  >  1:776253/1‑99 (MQ=255)
                                                                                                   |                                                                                                             
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705886‑3706083

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