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.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑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)
                  cACGCAGAGGCACAGCGCCACCGGGTATACATATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                           <  2:505950/100‑1 (MQ=255)
                   aCGCAGCGCCACAGCGCCACCGGGTCTACAGATTGTTGACAAAGCGCGCATCCGTTTCACCCGATGCGACGTTAACGCCTTATCCCGCCTCCAACATc                                                                                            >  1:1735371/1‑98 (MQ=16)
                   aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                          <  2:1707329/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                         >  1:2574668/1‑100 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                        <  2:2958851/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                        <  1:2846191/100‑1 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGTCCTACAAAATCGTGCa                                                                                       >  1:3459258/1‑100 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                       <  2:2011895/100‑1 (MQ=255)
                         agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                                    <  2:104045/100‑1 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                               <  2:2370798/100‑1 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                               <  2:401821/100‑1 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                               <  1:2910280/100‑1 (MQ=255)
                                  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                                           >  1:1187763/1‑100 (MQ=255)
                                   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                          >  1:112978/1‑100 (MQ=255)
                                    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                                         >  2:3558869/1‑100 (MQ=255)
                                       cGGGTATACAGATTGTTGACAAAGTGCTCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑Ag                                                                   <  1:998186/100‑1 (MQ=255)
                                       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑Ag                                                                   <  1:2811991/100‑1 (MQ=255)
                                           tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtt                                                              >  2:3153657/1‑100 (MQ=255)
                                            ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGTGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtt‑‑‑t                                                          <  2:138914/100‑1 (MQ=255)
                                            ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtt‑‑‑t                                                          <  1:3570604/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑g                                                     <  2:646015/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑g                                                     <  2:1187769/100‑1 (MQ=255)
                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑Gt                                                    >  2:3621596/1‑100 (MQ=255)
                                                 gcttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTa                                                   >  2:2531750/3‑100 (MQ=255)
                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTa                                                   <  1:3621568/100‑1 (MQ=255)
                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTa                                                   <  1:1220134/100‑1 (MQ=255)
                                                  aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAg                                                  <  1:2373010/100‑1 (MQ=255)
                                                   gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAgg                                                 >  1:2204182/7‑100 (MQ=18)
                                                        gACAAAGTGCGCTTTGTTTATGCCAAATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTg                                             <  1:3108082/99‑1 (MQ=255)
                                                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGAt                                           <  1:272825/100‑1 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATa                                          <  2:1156422/100‑1 (MQ=255)
                                                           aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATaa                                         <  1:848009/100‑1 (MQ=255)
                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGc                                       >  2:2984976/1‑100 (MQ=255)
                                                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                     >  1:3081429/1‑100 (MQ=255)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                   >  2:3616599/1‑100 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                  <  1:1773930/100‑1 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                 >  1:372536/1‑100 (MQ=255)
                                                                    tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                                >  2:1896797/1‑100 (MQ=255)
                                                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                               <  2:1435350/100‑1 (MQ=255)
                                                                     ttGTGTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAGAGTGCTAATTAAATACATTGC‑‑TGTATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                     <  1:440247/94‑1 (MQ=17)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                             <  1:1275269/100‑1 (MQ=255)
                                                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            <  1:1201295/100‑1 (MQ=255)
                                                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            >  2:500551/1‑100 (MQ=255)
                                                                        tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:2562065/6‑100 (MQ=18)
                                                                        tgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCATACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            >  1:1597669/6‑100 (MQ=18)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           <  1:1171033/100‑1 (MQ=255)
                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATGACAAAAACCTGCAAATTCAATATATTGC‑‑AGGATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  2:1528971/5‑100 (MQ=17)
                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAACCTGCAAATTCAATATATTGC‑‑AGGATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:494504/5‑100 (MQ=18)
                                                                         gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCATACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           >  1:2095719/5‑100 (MQ=17)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                          >  1:3339126/4‑100 (MQ=21)
                                                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2753248/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:3307256/3‑100 (MQ=17)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:3500583/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAAAA‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                        >  1:1938436/3‑99 (MQ=21)
                                                                            ggtcggataaggcgTTCACGCCGCATCCGGCCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAACG‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:2762599/11‑100 (MQ=18)
                                                                            ggtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAACG‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:1978380/11‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2468600/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                         >  2:3188461/2‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑AAGATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2696508/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAGTTCGTGCAAATTCAATAAATTGC‑‑AATAT‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga                       >  2:1880054/2‑98 (MQ=21)
                                                                             gtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAACG‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  2:1288420/10‑99 (MQ=21)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  2:1675406/1‑99 (MQ=21)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGAC‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag                       >  2:2581391/1‑99 (MQ=18)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:3334806/9‑98 (MQ=18)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:746243/9‑98 (MQ=17)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:2291139/9‑98 (MQ=18)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                      >  1:2325546/9‑98 (MQ=25)
                                                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:1962587/1‑100 (MQ=255)
                                                                              ccGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGC‑‑AGAAAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  1:1059579/1‑98 (MQ=17)
                                                                              ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  2:87263/1‑98 (MQ=18)
                                                                               cggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                     >  2:2603850/8‑97 (MQ=17)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1769186/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtat                    >  1:2271363/7‑96 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑‑AGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:1724518/7‑100 (MQ=18)
                                                                                ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    <  1:1389415/100‑1 (MQ=255)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAAAGAATCAGGCAgtttt                   >  2:1362151/6‑100 (MQ=17)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:1847113/100‑1 (MQ=255)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  2:774892/100‑1 (MQ=255)
                                                                                   tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 <  2:1774810/100‑1 (MQ=255)
                                                                                    gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                <  1:2951913/100‑1 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  2:3567587/1‑100 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:2109391/100‑1 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:2026018/100‑1 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:2218269/100‑1 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc          <  2:3461796/100‑1 (MQ=255)
                                                                                            aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATAGGGCATTTTTGCCTTTGCCAt        >  2:1604981/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:2246440/100‑1 (MQ=255)
                                                                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa      <  1:1441139/100‑1 (MQ=255)
                                                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  2:3389178/100‑1 (MQ=255)
                                                                                                   tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  >  1:776253/1‑99 (MQ=255)
                                                                                                   |                                                                                                             
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑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: