Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,017 C→T 31.3% intergenic (‑227/+81) 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,706,0170CT31.3% 441.5 / 224.1 298intergenic (‑227/+81)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (101/86);  new base T (83/2);  total (210/88)
Fisher's exact test for biased strand distribution p-value = 2.43e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 300 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGT  >  NC_000913/3705918‑3706114
                                                                                                    |                                                                                                          
  gcCACCGGTTATACAGATTGTTGACAAAGTGCCCTTTGTTTATGCCGGATGTGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                                                            <  2:3355360/98‑1 (MQ=255)
  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                                                                          >  2:2606558/1‑100 (MQ=255)
   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                                                         >  2:2841549/1‑100 (MQ=255)
   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                                                         >  2:3938411/1‑100 (MQ=255)
       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑a                                                                                                   <  1:169121/99‑1 (MQ=255)
         ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGaa                                                                                               >  2:2171039/1‑100 (MQ=255)
            ataCAGATGGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtt‑‑‑t                                                                                         <  1:341159/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑g                                                                                    <  2:2995149/100‑1 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTa                                                                                  <  1:3150194/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAgg                                                                                <  1:247495/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAgg                                                                                >  2:340911/1‑100 (MQ=255)
                    tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGc                                                                               >  2:4266395/1‑100 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAGAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGcc                                                                              >  1:3088648/1‑100 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGcc                                                                              <  2:715622/100‑1 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTGACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTg                                                                            <  1:2210155/100‑1 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTg                                                                            >  1:1683117/1‑100 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTg                                                                            <  2:1470226/100‑1 (MQ=255)
                       aaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTg                                                                            >  2:1834330/3‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGa                                                                           >  1:3744466/1‑100 (MQ=255)
                            aagcgcacTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAg                                                                       >  2:2574175/8‑100 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGc                                                                      <  1:1976622/100‑1 (MQ=255)
                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                 <  1:2371245/100‑1 (MQ=255)
                                    tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAgcgc                                                               >  2:209628/10‑100 (MQ=18)
                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                              <  2:1735579/100‑1 (MQ=255)
                                     ctttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAgcgcagcgca                                                              >  2:2700622/9‑100 (MQ=21)
                                      tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCAt                                                             >  2:1117994/8‑100 (MQ=21)
                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                             >  2:3772392/1‑100 (MQ=255)
                                      cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAgcgc                                                               >  2:945490/8‑98 (MQ=17)
                                       ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                           >  2:3553759/7‑100 (MQ=18)
                                       ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATc                                                            >  2:4343869/7‑100 (MQ=18)
                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                            >  1:141880/1‑100 (MQ=255)
                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                            <  2:3531107/100‑1 (MQ=255)
                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                           >  1:2461817/1‑100 (MQ=255)
                                        tgtaggCCGGATAAGGCGTTCACGCCGCATCAGGCCTACAAAGACCTGCAAATTCAATATATTGC‑‑AGGATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                           >  2:2981499/6‑100 (MQ=17)
                                        tgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                            >  2:2264423/6‑99 (MQ=17)
                                         gtaggCCGGATAAGGCGTTCACGCTGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                          >  2:645012/5‑100 (MQ=18)
                                         gtaggCCGGATAAGGCGCTTACGCCGCATCCGGCCTACAAGAACGTGCAAATTCAACAAATTGTG‑‑ATGCATC‑‑‑CT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                           >  2:515598/5‑100 (MQ=18)
                                          taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                         >  1:1483394/4‑100 (MQ=17)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                         >  1:3253062/4‑100 (MQ=18)
                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                         <  1:3801484/100‑1 (MQ=255)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGC‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                        >  1:4033404/3‑100 (MQ=17)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                        >  1:2863329/3‑100 (MQ=17)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                        >  2:380383/1‑100 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                       >  2:685377/1‑100 (MQ=255)
                                            ggctggataaggcgTTCACGCCGCATCCAGCCTACAAAATTGTGCAAAATCAATGGATTGC‑‑ACAGCT‑‑‑AGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  2:2329378/11‑100 (MQ=25)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCACTATATTGC‑‑AGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  1:732366/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  1:966751/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AATTCT‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                                                       >  2:868667/2‑99 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                        >  2:1218092/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGT‑‑GGGCTC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  2:2075891/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACT‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                                                       >  2:2470596/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAACCGCGCAAATTCAATATATTGCGGAGAAAA‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                      >  2:3562998/2‑98 (MQ=21)
                                             gtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAAC‑‑‑GAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                      >  2:3266383/10‑99 (MQ=21)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                      >  1:3624096/1‑100 (MQ=255)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                      >  1:4109146/1‑98 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                      >  2:3930721/1‑98 (MQ=18)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGC‑‑‑CAC‑‑GTAGTCCTGATAAGCGTAGCGCATCAGGCag                                                      >  1:1008699/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACa                                                       >  2:1010024/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                      >  2:168247/1‑99 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑‑AGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                      >  2:2227271/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                       >  1:3517417/1‑99 (MQ=17)
                                              ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGC‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                     >  2:1252887/9‑98 (MQ=17)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                                                     <  1:1117991/92‑3 (MQ=25)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑AAGATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                                                     >  2:4129596/9‑98 (MQ=17)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑‑AGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                     >  1:1203728/9‑98 (MQ=18)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGC‑‑AGAGAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                     >  2:4052065/9‑98 (MQ=14)
                                              ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                    >  1:1279166/1‑100 (MQ=21)
                                              ccGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                     >  2:4174774/1‑100 (MQ=18)
                                               cggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                    >  2:2769333/8‑97 (MQ=18)
                                               cggataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                    >  2:2418941/8‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGT‑‑GGGCTC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                   >  1:3320088/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCACTTAAATTCAATATATTGC‑‑ACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattc                                                   >  1:2987483/7‑99 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAAC‑‑‑GAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                                                   >  1:1517992/7‑96 (MQ=21)
                                                 gataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCATGATAAGCGGAGCGCATCAGGCAttttt                                                  >  1:1375457/6‑100 (MQ=18)
                                                   taaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                  >  1:1881978/4‑98 (MQ=18)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                  >  1:2104398/1‑98 (MQ=255)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                >  2:2022169/1‑100 (MQ=255)
                                                    gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                               <  1:1834323/100‑1 (MQ=255)
                                                           tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGCGCGCATCAGGCATTttttccttttc                                          >  1:1315856/1‑96 (MQ=255)
                                                           tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttttttttttttt                                        >  1:3867332/1‑90 (MQ=255)
                                                             aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                      <  1:3283167/100‑1 (MQ=255)
                                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                     >  1:2268796/1‑100 (MQ=255)
                                                                    aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATAAGGCATTTTTGCCTTTGCCATCAGTCtt                                >  1:571339/1‑99 (MQ=255)
                                                                      cctgCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                             >  2:1970570/4‑100 (MQ=255)
                                                                      ccGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                             <  1:846601/100‑1 (MQ=255)
                                                                      ccGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                             >  2:47104/1‑100 (MQ=255)
                                                                       cGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTCTTTGCTTTTAGTCGTGTac                            >  1:4230648/1‑99 (MQ=18)
                                                                       cGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                            >  1:2680948/1‑100 (MQ=255)
                                                                        ggCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                           >  1:531758/1‑100 (MQ=255)
                                                                            tACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                       <  1:2293729/100‑1 (MQ=255)
                                                                            tACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                       >  1:3912872/1‑100 (MQ=255)
                                                                              cAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                     >  1:1793546/1‑100 (MQ=255)
                                                                                  aaTCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGTCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                  <  2:442052/100‑1 (MQ=255)
                                                                                  aaTCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                  <  2:1883434/100‑1 (MQ=255)
                                                                                    tCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                <  2:2710519/100‑1 (MQ=255)
                                                                                      gTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt              <  2:3121656/100‑1 (MQ=255)
                                                                                          aaaTTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt          <  2:1347561/100‑1 (MQ=255)
                                                                                           aaTTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa         <  1:3185528/100‑1 (MQ=255)
                                                                                               cAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga     <  2:1372329/100‑1 (MQ=255)
                                                                                                  tACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTACTTATTCGATCGAGACGt  >  1:4062190/1‑100 (MQ=255)
                                                                                                    |                                                                                                          
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGT  >  NC_000913/3705918‑3706114

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: